V V Kafarov Experimentează metode de optimizare în chimie tehnologii contabilizarea e indemnizatiile d p * p ( , ) Probabilitatea unui anumit eveniment U este egală cu unu: P(C)= , ( , ) iar probabilitatea unui eveniment imposibil este AND-zero: Р(V) = ( , ) Prin urmare, J /(x)dx ( , ) —oo Din aceasta, se derivă o altă proprietate importantă a densității distribuției: oo p(x)dx= , -( întrucât lovirea unei variabile aleatoare în intervalul -°° O- ) = pentru continuu B[X] = J(x - mx) f(x)dx ( , ) -oo Alte denumiri pentru dispersia Dx, stx , o Rădăcina pătrată a celui de-al doilea moment central se numește abatere standard (sau standard)” °x = VD[X] = V G ( - ) Al treilea moment central, împărțit la x , se numește coeficient de asimetrie: •> i = • ( , ) În ceea ce privește momentele inițiale, J se exprimă astfel: R-z \u d tz - "it - t \\ ( , ) Al patrulea moment central este calculat prin formula P- \u d " - mjm + m m - Zt * ( , ) Valoare \u d ( H / ° x) - ( , ) se numește coeficient de curtoză Pe fig prezintă exemple de densități de distribuție cu coeficienți de asimetrie și curtoză non-zero Pentru comparație, linia întreruptă arată o curbă cu aceeași așteptare matematică tx și varianță σ , dar cu valori zero ale coeficienților de curtoză și asimetrie Momentele există dacă integralele sau seriile corespunzătoare pentru mărimi discrete converg Pentru variabile aleatoare ale căror valori sunt limitate, momentele există întotdeauna Dacă o variabilă aleatoare X are un prim și un al doilea moment ii, este posibil să se construiască o variabilă aleatoare normalizată: Xo \u d (X - mx) /? g ( - ) Pentru o variabilă aleatoare normalizată M[XO]= , D[XOJ= ( , ) Multe tabele de distribuții sunt construite special pentru variabile aleatoare normalizate Există următoarele relații între funcțiile de distribuție, respectiv - Orez Densitatea de distribuție cu coeficienți de asimetrie și curtoză non-nuli corespunzătoare mărimii normalizate Xo și nenormalizatei X\ / x—mx \ f(x)=-f (xt) = -f [ - , ( , ) fi(^o) = °x/(JC) = ax/("x + ox^o)- (I- ) F (x) = L (x ) = L ( X ~ OTx ) , ( , ) \ °x / L(*o) = F(x) = F(mx + °x*o)- ( - ) Momentele sunt caracteristici generale (integrale) de distribuție Al doilea grup de parametri caracterizează valorile individuale ale funcției de distribuție Acestea includ cuantile Cuantila xp a distribuției unei variabile aleatoare X cu funcția de distribuție F(x) este soluția ecuației F(xp) - p, ( , ) adică xp este o astfel de valoare a unei variabile aleatoare încât P(X b Deoarece aria mărginită de curba de distribuție este egală cu unitatea: c(b - a) = , atunci c*=\I(b-a), iar densitatea distribuției f(x) are forma f (x) = pentru a b ( - ) Funcția de distribuție este exprimată prin aria curbei de distribuție situată în stânga lui x Prin urmare, pentru x b ( - ) Graficul funcției F(x) a unei variabile aleatoare de pe segmentul [a, /f este egal cu prezentată în fig „Așteptarea ca X să aibă o distribuție uniformă pe a + b ( , ) Datorită simetriei distribuției uniforme, mediana lui X este, de asemenea, egală cu (a + b)/? Folosind formula ( ), determinăm varianța variabilei aleatoare X: (b - a) ( , ) Coeficientul de asimetrie Ui=Sch/stx este egal cu zero (distribuția este simetrică) Pentru a determina coeficientul de curtoză, găsim al patrulea moment central: b = - I x b - a ,] \ A (b - ") de aici g = - s \u d - Probabilitatea ca o variabilă aleatoare X distribuită uniform să lovească segmentul [a, p], care este o parte a segmentului [a, ] (Fig ), este determinată de raportul dintre lungimea segmentului [a, p] la lungimea întregului segment [a, b]; - a P (a Pentru o variabilă aleatoare normalizată, luând în considerare ( ) și ( ), avem P(X! , ) Să găsim probabilitatea evenimentului opus folosind formula ( ): P(Ni xt F(xr, y)^F(xlt y)-, ( - ) cu y >Yi F (x, yt)^F (x, yj Peste tot pe -°° funcția de distribuție este egală cu zero: F(x, -oo) - F(-oo, y) = F(-oo, -oo) = ( , ) Cu unul dintre argumente egal cu +°°, funcția de distribuție a sistemului se transformă într-o funcție de distribuție a unei variabile aleatoare corespunzătoare celuilalt argument: F (x, w) = Ft (x), ( , ) F (+ o °, Y) \u d (Y), unde F}(x), F (y) sunt funcțiile de distribuție ale variabilelor aleatoare X și respectiv Y Dacă ambele argumente sunt +°°, atunci funcția de distribuție a sistemului este egală cu unu: F(+°° + x, c^ sunt abaterile standard ale valorilor lui X și Y Acest indicator se numește coeficient de corelație al valorilor lui X și Y Pentru variabile aleatoare independente, coeficientul de corelație este zero, dar invers nu este adevărat - coeficientul de corelație (și covarianța) poate fi egal cu zero, iar variabilele aleatoare sunt dependente: legătura, fără a afecta variațiile, se manifestă în momente de ordin superior În cazul general, o afirmație mai slabă este adevărată: variabilele aleatoare pentru care covarianța (și, prin urmare, coeficientul de corelație) este egală cu zero, se numesc necorelate Astfel, variabilele aleatoare independente sunt întotdeauna necorelate, dar independența lor nu rezultă din necorelarea variabilelor în cazul general Și numai în cazul unei distribuții normale, egalitatea coeficientului de corelație la zero indică fără ambiguitate absența unei legături între variabilele aleatoare Densitatea distribuției normale a două variabile aleatoare este exprimată prin formula ( ) unde r este coeficientul de corelație Să presupunem că variabilele aleatoare X, Y, respectând legea distribuției (I ), sunt necorelate, adică r = Obținem (x - tx)g (y-ty) f(x, y)= -exp ts Oy a k exp (x - mx) (U - tu) g = fi (*) fi (Y) ( ) adică, densitatea de distribuție a sistemului este egală cu produsul densităților de distribuție ale cantităților individuale incluse în sistem, ceea ce înseamnă că variabilele aleatoare X și Y sunt independente Astfel, pentru un sistem de variabile aleatoare supus legii normale, coeficientul de corelație zero indică nu numai necorelare, ci și independența variabilelor aleatoare, astfel încât importanța rolului coeficientului de corelație ca indicator al conexiunii în acest caz crește semnificativ Observați următoarele proprietăți ale coeficientului de corelație Coeficientul de corelație nu se schimbă de la adăugarea oricăror termeni non-aleatori la X și Y, de la înmulțirea X și Y cu numere pozitive Dacă una dintre valori, fără modificarea celeilalte, este înmulțită cu - , atunci și coeficientul de corelație va fi înmulțit cu - Prin urmare, coeficientul de corelație rxy nu se va modifica dacă trecem de la variabilele aleatoare originale la cele normalizate: X - mx Y - Shu Pe baza ( ) și ( ) avem Z)(X+K) = D(X) + Z> (Y)+ yx Vd(X)D(Y) ( ) În mod similar, pentru varianța diferenței a două variabile aleatoare, putem scrie D (X ~Y) - D (X)D (Y) ~ ryx VD(X)D(Y) ( ) Expresiile ( ) și ( ) pentru variabile aleatoare normalizate, ținând cont de faptul că D(X )=D(Y ) = , iau forma D (Xo + Uo) \u d + guh, D (X - Uo) \u d - guh ( ) Deoarece dispersia este o cantitate nenegativă, avem + yx > ; - guh ; guh și, în sfârșit ( ) Valorile extreme ale coeficientului de corelație yx =± corespund dependenței funcționale liniare Y \u d bl + bіX, iar semnul coeficientului b, corespunde semnului coeficientului de corelație În cazul general, când valorile Xi sunt legate printr-o dependență stocastică arbitrară, coeficientul de corelație poate avea o valoare în intervalul - >> există o corelație pozitivă între mărimile X și Y, iar pentru rxy !(/) = - " " - exp K - Gg Kgl Din ( ) avem f (Y *) \u d - ■ - exp □u I -g V k / x - mx y - my \ ' (i-^I ~ ~r y I • ( ) G ° y -:—- y - tu - g - (x - tx) ( - r ) o y ohv ( ) Evident, ( ) este o expresie pentru densitatea legii distribuției normale cu așteptări matematice la meu\x \u d meu + r ~ (x-mx) ( ) și abaterea standard °g|x = =y( - f ) • ( - ) Mărimea mYIX se numește așteptarea matematică condiționată a mărimii Y pentru un X dat Dependența liniară ( ) se numește regresia lui Y pe X În mod similar, linia dreaptă A, txIu~ tx + g—(U—tu) (I P ) °y există o regresie a lui X pe Y Dreptele de regresie coincid doar dacă există o dependență funcțională liniară a lui Y de X Cu X și Y independente, dreptele de regresie sunt paralele cu axele de coordonate În acest capitol au fost luate în considerare principalele caracteristici ale variabilelor aleatoare Informațiile complete despre variabilele aleatoare sunt conținute în legile distribuției Sunt luate în considerare legile uniforme și normale ale distribuției probabilităților În multe probleme aplicate, nu este nevoie să folosiți legile de distribuție În schimb, puteți folosi numere caracteristicile unei variabile aleatoare, exprimând într-o formă comprimată cele mai semnificative trăsături ale distribuției În teoria probabilităților și statistica matematică, se utilizează un număr mare de caracteristici numerice Acest capitol introduce conceptul de momente de distribuție, notează proprietățile momentelor cele mai frecvent utilizate - așteptarea și varianța matematică Se introduce conceptul de relație stocastică între variabile aleatoare și coeficientul de corelație, care caracterizează apropierea unei relații liniare între variabile aleatoare Studiul dependenței dintre variabile aleatoare este o problemă aplicată importantă Exerciții O sută de tije polimerice noi sunt testate aleatoriu pentru rezistență Condiția de neadecvare a întregului lot este prezența a cel puțin unei tije insuficient de puternice dintre cele cinci testate Care este probabilitatea ca un anumit lot să fie respins dacă conține % tije insuficient de puternice? O variabilă aleatoare X este supusă unei legi de distribuție a cărei densitate de probabilitate are forma /(*)= ax pentru Determinați așteptările matematice, varianța, abaterea standard și asimetria distribuției Funcția de distribuție a unei variabile aleatoare continue X este dată de expresia FW= O pentru x Aflați: a) coeficientul a; b) densitatea distribuţiei f(x) \ c) probabilitatea ca valoarea X să cadă în intervalul ( , ; , ) Erorile de măsurare sunt o variabilă aleatoare supusă legii normale cu parametrii m - , ; ' = , Găsiți probabilitatea ca această valoare să se încadreze în intervalul (- , ; + , ) Ca urmare a verificării acurateței dispozitivului, s-a constatat că % dintre erori nu depășesc ± ° С Determinați eroarea pătratică medie a dispozitivului dacă se știe că dispozitivul nu are erori sistematice, iar erorile aleatoare sunt distribuite conform legii normale CAPITOLUL II DETERMINAREA PARAMETRILOR FUNCȚIEI DE DISTRIBUȚIE Populația generală și eșantionarea aleatorie În practică, cercetătorul are întotdeauna doar un număr limitat de valori ale unei variabile aleatorii, care este un anumit eșantion din populația generală Populația generală este înțeleasă ca toate valorile admisibile ale unei variabile aleatorii Când se analizează orice variabilă tehnologică aleatoare care se schimbă continuu timp (de exemplu, temperatura, presiunea etc n ), valorile observate ale unei variabile aleatorii sunt înțelese ca valori ale tehnologice parametru la momente discrete separate de un astfel de interval la care valorile învecinate pot fi considerate obținute din experimente independente Eșantionul se numește reprezentativ (reprezentant) dacă oferă o idee suficientă a caracteristicilor populației generale Dacă nu se știe nimic despre populație, selecția aleatorie este singura garanție a reprezentativității În foarte multe studii, este necesară selecția aleatoare sau amestecarea aleatorie (randomizarea) a datelor Tabelele cu numere aleatorii pot fi folosite pentru a simula selecția aleatoare Să presupunem că trebuie să selectați articole dintr-o populație care conține de articole Pentru a face acest lucru, trebuie să numerotați elementele populației generale de la la Apoi, începând din orice loc din tabele, scrieți ultimele două cifre din zece numere consecutive De exemplu, începând de la primul număr, numerele (dacă numerele se repetă, acestea trebuie omise) Numerele rezultate arată ce elemente trebuie selectate Secvența selectată nu poate fi modificată Încălcarea aleatoriei, de regulă, duce la o denaturare a rezultatelor Similar cu selecția, elementele sunt randomizate În acest caz, trebuie să scrieți numere aleatorii până când acestea acoperă toate elementele date Din natura aleatorie a eșantioanelor, rezultă imediat că orice judecată despre populația generală din eșantion este în sine aleatorie Să presupunem că în urma experimentului se obține un eșantion x}, x , , xn de valori ale unei variabile aleatoare X Fie x un punct pe axa reală x; notăm cu nx numărul de valori ale eșantionului situat în stânga lui x pe aceeași axă Raportul nx/n este frecvența valorilor variabilei aleatoare X obținute în eșantion care sunt mai mici decât x Această frecvență este o funcție a lui x Se notează cu F„(x): Fn(x) = nx/n (BUNĂ) Funcția de distribuție Fn(xy obținută din eșantion se numește funcție empirică sau de distribuție a eșantionului (în contrast cu distribuția populației generale, sau distribuția teoretică) Pentru fiecare eșantion, funcția de distribuție empirică va fi diferită, dar toate empirice funcțiile de distribuție ale aceleiași variabile aleatoare vor avea ceva general, care este informații despre funcția de distribuție a acestei variabile aleatoare Se poate demonstra (teorema lui Glivenko) că cu probabilitatea ca n -> diferența maximă dintre funcțiile de distribuție ale variabilelor aleatoare Fn(x) și F(x) tinde spre : P(sup |F(x)-Fn(x)| -> ) = , (II ) L—>OO —OO A \u d , i - G, P (II ) fn(x) = pentru X > xn Pe fig prezintă graficul funcţiei Fn(x) Toate elementele eșantionului se dovedesc a fi puncte de discontinuitate ale acestei funcții În punctul de discontinuitate x-xk, funcția Fn(x) sare de la valoarea (k- )/n (în intervalul x (pag ) În cazul general, este necesară estimarea simultană a mai multor parametri ai unei distribuții unidimensionale sau multivariate Dacă a și x sunt înțeleși ca vectori, atunci formularea principiului probabilității maxime va fi păstrată; este necesar să se găsească un astfel de set de valori admisibile ale parametrilor a^ , a* , aț , care inversează funcția de adevăr asemanare cu maxima Condițiile necesare pentru un extremum sunt date de sistemul de ecuații dTskh, ah, ah) — - = aproximativ / = , , (I ) da j în timp ce caracterul nenegativ al matricei /&L\ , i, j= , , (pag ) \ da da este o condiție suficientă pentru ca extremul local să fie maximul funcției de probabilitate Să găsim, prin metoda maximei probabilități, o estimare a parametrului A a distribuției exponențiale cu densitate / (x) = Ae x, , din ( ) avem (xz - m) = , de unde găsim estimarea - ("" x pentru parametrul m\ Diferențiând funcția de probabilitate în raport cu o , obținem p o" l (хі - m) = , (pag ) P Deoarece /( o ) = , avem P - P i = , de unde găsim estimarea R: S = P (pag ) Estimarea așteptării și dispersiei matematice Metoda probabilității maxime conduce întotdeauna la estimări consistente, deși uneori părtinitoare, care au cea mai mică variație posibilă, pe măsură ce dimensiunea eșantionului crește fără limită Pentru o variabilă aleatoare distribuită normal, se obțin estimări de următoarea formă: medie aritmetică ~x pentru așteptarea matematică tx P x— Y, Xi/P = și varianța eșantionului n pentru varianța D [X] Ultima estimare este oarecum părtinitoare: M G s?l = D[X] (pag ) l p - Pentru a obține o estimare imparțială, trebuie înmulțit cu n/(u—): l p ( , ) Scăderea numitorului din ( ) cu unu este direct legată de faptul că valoarea lui x, în raport cu care se iau abaterile, depinde ea însăși de elementele eșantionului Fiecare valoare care depinde de elementele eșantionului și este inclusă în formula varianței eșantionului se numește relație Se poate dovedi că numitorul varianței eșantionului este întotdeauna egal cu diferența dintre dimensiunea eșantionului n și numărul de legături / impus acestui eșantion Ea diferență f = n - l ( , ) se numește numărul de grade de libertate ale eșantionului În calculele practice, pentru varianța N, formula este adesea convenabilă care decurge cu ușurință din ( ): ( , ) \u d n \ [ [( - xjx + x ) + (\u e? - x x + x ) + - + (x p~ xpx+ x )] = Ț [( + ••• +^)- - x (Xi + x + - xn) + (x + x N + x )] = Avantajul formulei ( ) este că nu conține operații de scădere a numerelor apropiate, ca în formula ( ), ceea ce duce la o pierdere a preciziei În formula ( ), această operație se aplică o singură dată Media și varianța eșantionului conform datelor grupate din tabel se calculează prin formule kk ( , ) ( ) Valoarea #/ se numește corecția lui Sheppard, este legată de schimbarea varianței în timpul grupării Clasificarea erorilor de măsurare Fiecare rezultat al măsurării este o variabilă aleatorie Abaterea rezultatului real de la cel adevărat se numește eroare de observație Eroarea de observare este, de asemenea, o variabilă aleatorie - este rezultatul acțiunii doar a unor factori aleatori (neluați în considerare) Dacă notăm rezultatul adevărat cu a, eroarea cu Ai, rezultatul măsurării cu X, atunci X - a = bx ( , ) Există trei tipuri de erori Erorile grave apar din cauza încălcării condițiilor de bază de măsurare Rezultatul, care conține o eroare grosolană, diferă brusc ca magnitudine de restul măsurătorilor Unele criterii pentru excluderea erorilor grave se bazează pe aceasta Erorile sistematice sunt constante în întreaga serie de măsurători sau se modifică după o anumită lege Dezvăluirea lor necesită studii speciale, dar de îndată ce sunt descoperite erori sistematice, acestea pot fi eliminate cu ușurință prin introducerea de corecții adecvate în rezultatele măsurătorilor Erori aleatorii - erori de măsurare rămase după eliminarea tuturor erorilor grosolane și sistematice identificate Cu această definiție, factorii aleatori care generează o eroare aleatoare nu includ factorii cu o acțiune constantă (erori sistematice) și factorii cu o singură acțiune, dar foarte puternică (erori grosolane) Erorile aleatorii sunt cauzate de un număr mare de astfel de factori, ale căror efecte sunt atât de mici încât nu pot fi izolate separat (cu nivelul actual de tehnologie de măsurare) În acest caz, distribuția erorilor aleatoare este simetrică față de zero: erorile care sunt opuse ca semn, dar egale ca valoare absolută apar la fel de des Din simetria distribuției erorilor rezultă că adevăratul rezultat al observației este așteptarea matematică a variabilei aleatoare corespunzătoare Deoarece de la ( ) X \u d a + AX și în absența erorilor grosolane și sistematice M [DH] = , ( , ) Acea L [X] = a (II ) În cele ce urmează, vor fi luate în considerare numai erorile aleatorii de măsurare Legea adunării erorilor Pentru variabile aleatoare independente, varianțele au proprietatea de aditivitate, dar nu și erorile pătratice medii Dacă Xp X , , Xn sunt variabile aleatoare independente; a}, a , , an sunt variabile non-aleatoare și Z - aiXl -f-arXr + • • • + anXn, (II ) G atunci varianța eșantionului a lui Z este determinată după cum urmează: „ , II „ SZ - " SX + a SX, + -T ansX ( , ) Dacă punem a, = a = = an = /n, atunci Z= În acest caz n Sy SȚ - n ns (II ) P si - î= unde sx = - n Dacă valorile Xv X , Xn sunt interpretate ca n observații independente ale aceleiași variabile aleatoare X, atunci sț = = = d = Atunci obținem s^ =s> x/n ( , ) Din expresia ( ) rezultă o concluzie practică: atunci când se evaluează acuratețea celor două metode, trebuie luată în considerare durata analizei Folosind o metodă expresă mai puțin precisă, se poate face un număr semnificativ mai mare de experimente în același timp și se poate obține o precizie mai mare decât o oferă metoda exactă consumatoare de timp Erori de măsurători indirecte Măsurătorile sunt împărțite în linii drepte și indirect În primul caz, mărimea de determinată este măsurată direct; în măsurătorile indirecte, este dată de o anumită funcție a mărimilor măsurate direct Fie variabila aleatoare z să depindă de observațiile xt, xn conform legii cunoscute: z = /(*i, xr, , xn) Valoarea adevărată a lui z poate să nu coincidă cu așteptarea matematică a lui A/ , dar este determinată de aceeași lege: az^f(™, t , t \ Valoarea lui a se numește media măsurării indirecte Varianta măsurătorii indirecte se determină în același mod ca și varianța obișnuită, doar abaterile se iau de la media măsurării indirecte ar Se poate constata prin cunoașterea varianțelor observațiilor individuale și a formei funcției / În practică, se determină variațiile eșantionului , iar din acestea varianța eșantionului a măsurătorii indirecte $ , care servește ca estimare a generalului varianta a? Pentru a găsi $, extindem funcția z-=f(xv x , xn) într-o serie Taylor în punctul (mx Numărul de grade de libertate pentru variația totală a reproductibilității, determinat de formulele ( ) și ( ), este mult mai mare decât pentru fiecare varianță particulară separat Prin urmare, varianța de reproductibilitate totală estimează mult mai precis varianța populației st resp Când se calculează variația de reproductibilitate pentru măsurătorile curente, sunt combinate numai acele eșantioane care pot fi considerate eșantioane din populații cu variații egale Mai mult, fiecare dintre valorile, •••, poate fi considerată o estimare a aceleiași variații generale Exemplul Rezultatele determinării concentrației (%) de P O în sistem (NH^sHPCh - K CO - H O prin metoda colorimetrică sunt prezentate în tabel Determinați eroarea metodei colorimetrice din măsurătorile curente Numărul eșantionului , , , , , , , , , , , , , , , , , - , , - , , , bUu , , , , , , , , , , , , , , , , tі Soluţie Conform tabelului, calculăm s? conform formulei sj GL* - Vai ~ - ti Pentru a calcula varianța totală folosind formula ( ), avem nevoie de termeni ai formei Si (mi - /mi Ui \u=l Ui— - Pentru primul test, avem (mj- )^ = , - ' ' = , În mod similar (t - ) « -= , ; (m - ) s| = , ; (m - ) = , ; (m -l)sj= , ; (m, - ) sg = , ; (m, - ) s^ = , ; (me-l)s| = , Numărul de grade de libertate a variației de reproductibilitate totală este /play= A= - = - = i=i "= Și variația de reproductibilitate S \u d - \u d \u d , / U PLAY/WOS PR Eroarea metodei colorimetrice, determinată din măsurătorile curente, este egală cu saocnp = SBocnp = V > = , Intervalele de încredere și probabilitatea de încredere Parametrii eșantionului sunt variabile aleatoare, abaterile lor față de cei generali (erori) vor fi și ele aleatorii Evaluarea acestor abateri este probabilistică în natură - se poate indica doar probabilitatea unei anumite erori Pentru a face acest lucru, statisticile matematice folosesc intervale de încredere și probabilități de încredere Să se obțină o estimare nepărtinitoare a* din experiență pentru parametrul general a Este necesar să se evalueze posibila eroare în acest caz Să atribuim o probabilitate suficient de mare — astfel încât un eveniment cu o probabilitate p poate fi considerat practic sigur și să găsim o astfel de valoare е =/( р ) = , pentru care Р( I а*-а I , parametrii necunoscuți sunt înlocuiți în expresia estimărilor lor; ) de la variabila aleatoare a* trec la o altă variabilă aleatoare, a cărei lege de distribuție nu depinde de parametrul estimat a, ci depinde doar de mărimea eșantionului n și de tipul legii de distribuție a mărimii X Astfel de mărimi au fost studiate în cele mai multe detalii pentru distribuția normală a variabilei aleatoare X Ca limite de încredere a' și a" iau de obicei cuantile simetrice * * „( - )/ sau luând în considerare ( ) ap/ Ѳр, se respinge ipoteza H și se acceptă ipoteza alternativă H și invers, dacă Ѳ aț și a* aț sau a* și Ѳ , , adică semnificative (nealeatoare) sunt considerate Ѳ*, care au luat valorile * și Ѳ*>> ^ , • Cu un test de semnificație unilaterală, una dintre aceste inegalități (de exemplu, * va fi semnificativă Probabilitatea ultimei inegalități este , și, prin urmare, nivelul de semnificație va fi , Astfel, dacă pentru testul de semnificație cu o singură coadă sunt utilizate aceleași numere critice ca și pentru cel cu două cozi, aceste valori vor corespunde la jumătatea nivelului de semnificație De obicei, același nivel de semnificație este luat pentru un test cu o coadă ca și pentru un test cu două cozi În aceste condiții, ambele criterii oferă aceeași eroare de primul fel Pentru a face acest lucru, un test cu o singură coadă trebuie să fie derivat dintr-un test cu două cozi corespunzător unui nivel de semnificație de două ori mai mare decât cel acceptat Pentru a menține nivelul de semnificație p = , pentru un test cu o singură coadă, p = , pentru un test cu două cozi, care oferă valori critice de - , și - Dintre acestea critice va exista o singură valoare pentru un test unilateral, de exemplu o,e • Nivelul de semnificație pentru un test unilateral este , Același nivel de semnificație pentru testul cu două cozi corespunde valorii critice Ѳ , Dar , ceea ce înseamnă că, cu un criteriu unilateral, un număr mai mare de ipoteze va fi respins și, în consecință, eroarea de al doilea fel va fi mai mică Estimarea așteptării matematice a unei variabile aleatoare distribuite normal În absența erorilor grosolane și sistematice, așteptarea matematică a unei variabile aleatoare coincide cu rezultatul adevărat al observațiilor Prin urmare, estimarea așteptărilor matematice este de mare importanță în prelucrarea observațiilor Cel mai simplu mod de a estima așteptările matematice este cu o varianță cunoscută a populației generale (vezi capitolul P ) Varianța generală și în formula '( ) pentru intervalul de încredere, parametrul general Folosind distribuția Student, mai întâi determinăm regiunea critică pentru un criteriu cu două fețe La p - , p - , și cuantila fj-p/ - - , la ffl (Tabelul din apendice) Valorile critice ale ipotezei nule conform ( ) vor fi >izh+ /] r/ sx/V~n- Doar prima inegalitate are sens fizic: x Xp- ( , ) Cuantilele Z p pentru diferite puf sunt date în tabelul aplicații Deoarece varianța eșantionului z prin elementele eșantionului este determinată de formula 'de la ( ) avem X =/"x/'G (II ) Înlocuind ( ) în ( ), obținem vulpe " Xp/ , * , , , Extragând rădăcina pătrată din toate părțile inegalităților, obținem o estimare pentru eroarea de reproductibilitate , * o* * , Datorită numărului mic de grade de libertate, limitele de încredere s-au dovedit a fi puternic asimetrice Pe măsură ce numărul de grade de libertate crește, asimetria curbelor de distribuție scade, iar asimetria limitelor de încredere scade în mod corespunzător Se poate demonstra că pentru n > eșantionul standard s este distribuit aproximativ normal cu așteptarea matematică ms = a și eroarea standard ', = 'x// - (P- ) Un standard general necunoscut în expresia (A ) cu n> este înlocuit cu unul selectiv; ,s^&xiyTf ■ ( , ) În conformitate cu ( ), limitele de încredere pentru standardul general sunt determinate de inegalitatea ss sx~^ ^ u'~PI 'x * s* + ui-PIi • (P- ) Exemplul Estimați eroarea de reproductibilitate ox pentru un eșantion de de observații cu un eșantion standard sx - , Probabilitatea de încredere p este luată egală cu , Soluţie Construim un interval de încredere pentru eroarea de reproductibilitate folosind distribuția X Conform tabelului aplicații, cu numărul de grade de libertate / - și probabilitatea de încredere p - , , găsim Xo os - , și Xo - , Conform formulei (A ), determinăm estimarea încrederii pe două părți pentru varianța de reproductibilitate: , • , • , , ' , , Limitele de încredere pentru eroarea de reproductibilitate sunt determinate de inegalitatea , * Ox* , Să determinăm limitele de încredere pentru ax folosind distribuția normală Prin formula ( ) definim „ca: yîf-}l • Conform tabelului aplicații pentru probabilitatea de încredere p - , găsim „ , - , În conformitate cu ( ), limitele de încredere pentru eroarea de reproductibilitate sunt determinate de inegalitatea , - , • , >/ unde Γ(/) este funcția gamma F-distribuția depinde numai de numărul de grade de libertate și f Pe fig prezintă curbele de densitate de probabilitate ale distribuției F pentru unele valori ale / și / Curbele sunt asimetrice În tabel Anexa prezintă cuantilele Fr p pentru nivelul de semnificație p= , și numerele de grade de libertate și / Pentru a determina cuantilele Fp pentru valorile p, se folosește relația Fp(J '= (A, L)' (II ) Sub ipoteza nulă = oȘ și oȘ / o? = și, prin urmare, distribuția F poate fi utilizată direct pentru a estima raportul variațiilor eșantionului $ /s| Cu un nivel de încredere de - p, estimarea de încredere bifață a lui F are forma fz) F Ff—p/zifit fz)- Luând în considerare ( , ) -P/ (A, A) (IL ) Sub ipoteza nulă F=^/^ , prin urmare, cu probabilitatea - p, inegalitatea cu două laturi trebuie satisfăcută sî - - - st , adică dacă o varianță mai mare a eșantionului evident nu poate corespunde uneia generale mai mici În acest caz, diferența dintre dispersiile conform (IL ) trebuie considerată semnificativă dacă fi-P (A, Al ții ) Valoarea lui F\ p(fi, fi) pentru p = , poate fi determinată din tabel aplicații Testul de semnificație cu două cozi (IL ) este utilizat pentru ipoteza alternativă > a ■ O sț Ft pf (f , fi) pentru râuri mici Diferența dintre variații ar trebui considerată semnificativă dacă ''î/ > ^ -p/ ^ > ^)- (I- ) Exemplul La evaluarea acurateței determinării conținutului de Р О asimilabil într-un îngrășământ complex prin metoda acidului sulfuric, dispersia de reproductibilitate a fost R = , ; f = Este necesar să se compare această metodă de analiză a P O digerabil cu o metodă mai precisă a citratului bazată pe rezultatele a patru determinări paralele de P O : , ; , ; , ; Soluţie Dispersia de reproductibilitate a metodei citratului cu numărul de grade de libertate fi-З În funcție de condițiile problemei, pentru a evalua semnificația diferenței dintre varianțele y și r, , poate fi utilizat testul de semnificație unilaterală ( ) Raportul de dispersie F- , / , - , trebuie comparat cu cel tabelar pentru nivelul de semnificație p = , și numerele de grade de libertate / i = și / - Fi - p (J \, fi) - Astfel, raportul de dispersie a probei este mai mic decât cel tabelar, iar datele experimentale nu ne permit să considerăm că acuratețea metodelor este semnificativ diferită Testul lui Fisher poate fi folosit pentru a compara variațiile dacă una dintre variații este generală Numărul de grade de libertate ale dispersiei generale este considerat egal cu Comparația mai multor varianțe Când se determină estimarea varianței pentru măsurătorile curente folosind formula /a? + / ' + ' ' ' + In^n ' $ — , — Y + / + • ■ ■ + fn f a fost acceptată ipoteza nulă a egalităţii varianţelor generale corespunzătoare Această ipoteză poate fi testată pentru mostre de dimensiuni diferite utilizând testul Bartlet Bartlet a arătat că în ipoteza nulă, raportul B/C, unde B \u d , (/lg ^ ~ fi lg ^ ), \ = / ( , ) distribuit aproximativ ca X cu u - grade de libertate dacă toate f > Ipoteza egalității varianțelor generale este acceptată dacă D/S xf /), criteriul Bartlett este calculat complet Exemplul în producerea fosforului prin sublimare din fosfaţi, gradul de reducere a fosforului a fost măsurat la patru temperaturi diferite Tabelul prezintă rezultatele unei analize statistice a omogenității dispersiunilor de reproductibilitate a rezultatelor la diferite temperaturi Temperatura /' Igs? Mg ? ' g , , , , , b , , , , , Tz , , , , , , , , , , , , , Determinați dacă acuratețea analizei se modifică cu temperatura Soluţie Conform tabelului, variația de reproductibilitate este egală cu \u d , / \u d , , lg * \u d , , /lg = • , = , , |// = / = , B \u d , ( , - , ) \u d , , Conform tabelului Găsim aplicații la trei grade de libertate și un nivel de semnificație de P = , , Xo = , Valoarea LC este , și, prin urmare, la nivelul de semnificație p = , , putem accepta ipoteza egalității varianțelor generale Valoarea lui C nu a putut fi calculată Astfel, criteriul Bartlett ne permite să considerăm că acuratețea analizei nu depinde de temperatură Varianțele eșantionului sunt omogene, prin urmare, ca estimare a varianței de reproductibilitate, putem lua varianța medie ponderată $ cu numărul de grade de libertate / egal cu Dacă variațiile eșantionului sunt obținute din eșantioane de dimensiuni egale mx = m = = mn = m, se folosește un test Cochran mai convenabil și mai precis pentru a le compara Cochran a examinat distribuția varianței maxime a eșantionului la suma tuturor variațiilor: /-g, ( ) Distribuția variabilei aleatoare G depinde numai de numărul de varianțe însumabile n și de numărul de grade de libertate f cu care este definită fiecare varianță: f = m - În tabel Anexa prezintă cuantilele G[p pentru nivelurile de semnificație p = , Dacă valoarea criteriului Cochran găsită din varianțele eșantionului este mai mică decât cea tabelară G °Х , ° у „/ , \ a = a- a-= - -= a -| i ( ) x U p \ Pі Pі ) Compuneți o variabilă aleatorie normalizată - t- (x - u) - (tx - ti) -i- = ——, ( ) ar R K M + /" care are o distribuție normală standard Dacă standardul general st este înlocuit cu cel eșantion, obținem o valoare care are o distribuție Student (ІІ ) s К М + Ms cu numărul de grade de libertate f egal cu / = u, + « Omogenitatea variațiilor eșantionului sx și y* poate fi verificată folosind testul Fisher Cu o probabilitate de încredere B = -p, avem o estimare cu două părți pentru diferența mx-tu: x Y 'i-zg^y P + n ~-y- t^ps I/" /A! + ! /n ( ) În ipoteza nulă, mx = mu, iar inegalitățile ( L ) și (II ) oferă un criteriu pentru testarea acestei ipoteze Ipoteza nulă este respinsă într-un test cu două cozi dacă | x - y / > t cu un criteriu unilateral, dacă I x - y I\u e ii-p ^ Ui / fh + / pg ( , ) Aceste criterii nu pot fi utilizate dacă Criteriile de mai sus nu pot fi utilizate dacă variațiile generale t ( , ) ( , ) ( , ) Criteriul formulat este bifat; devine unilateral când p/ este înlocuit cu p Exemplul S-a studiat procesul de polimerizare radicalică a sărurilor pe bază de -vinilpiridină și două halogenuri de alchil diferite: iodură de butii - C H I și bromură de etil - CsHsBr Comparați reactivitatea halogenurilor de alchil dacă randamentul mediu de polimer (n - ) în timpul sintezei cu C H I a fost , % (x), iar randamentul mediu (n - ) cu CrH Br este de , % (y) Eroarea de reproductibilitate a procesului de polimerizare este de - , % Soluţie Ca ipoteză nulă, se consideră ipoteza egalității reactivității alchililor halogenuri Numărul de grade de libertate al dispersiei de reproductibilitate este /- + - - Deoarece datele medii sugerează că reactivitatea bromurii de etil este mai mare, testul cu o coadă ( , ) poate fi utilizat pentru a evalua semnificația diferența La / - și p - , t , - , De aceea G-ps/ /^ - /p, = , - , V / + / = , ; y - x - , > , Prin urmare, la un nivel de semnificație de %, ipoteza nulă este respinsă și diferența de reactivitate a halogenurei de alchil ar trebui considerată semnificativă Dacă nu se fac ipoteze că reactivitatea bromurii de etil este mai mare, trebuie utilizat un test cu două cozi ( ) pentru a testa ipoteza nulă La / - si p - fo - De aceea , Astfel, chiar și cu un test cu două cozi, ipoteza nulă este respinsă Adesea, în practică, un eșantion de observații este alcătuit din mai multe subgrupuri obținute într-o ordine sau alta (de exemplu, din diferite părți ale populației generale) Pentru a combina astfel de subgrupuri într-un singur eșantion, este necesar să vă asigurați că mediile pe subgrupuri sunt omogene Pentru a face acest lucru, verificați semnificația diferenței dintre mediile subgrupurilor și media generală a întregului eșantion conform testului Student Să fie k subgrupuri de valori ale unei variabile aleatoare X de volum /i,, m , ,mk Pe baza acestor date, mediile subgrupurilor xp sunt media întregului eșantion x și pătratul mediu abaterea totală a eșantionului s: rnj p xh * - = - = unde n = /u, + /u + + mk; Denota W=(jy-*)/"• ( ) Se poate dovedi că valoarea (P І) p - t ] - mjy^ are o distribuție Student cu f=n- grade de libertate Folosind acest criteriu, se testează ipoteza nulă H° egală cu valorile așteptărilor matematice din subgrupul y-a m și așteptările matematice ale întregului eșantion m: = t Ipoteza nulă este respinsă dacă inegalitățile (II ) pentru un test cu două cozi sau inegalitățile ( ) și ( ) pentru un test cu o singură coadă nu sunt satisfăcute De obicei, ca medie x) luați în considerare cea mai mare (cea mai mică) medie dintre subgrupurile medii Compararea mai multor medii Când comparați mai multe mijloace, puteți utiliza testul f, comparând în perechi Cu toate acestea, pentru a fi utilizată în compararea informațiilor complete despre toate mijloacele, o astfel de comparație este efectuată folosind testul de rang multiplu al lui Duncan Să se obțină k valori medii pentru k eșantioane de diferite dimensiuni: Xk xg, xg, , X], pі X Пі Dispersiile generale sunt egale între ele, adică Atunci când se aplică criteriul Duncan, trebuie: ) să se ierarhească la valorile medii, ordonându-le în ordine crescătoare; ) determinați eroarea de reproductibilitate sx cu numărul corespunzător de grade de libertate fx; ) determinați eroarea pentru fiecare medie: = j/" s xlnj , j = , , , k; ) scrieți din tabelul Duncan (tabelul din apendice) (k-\) valorile de rang cu nivelul de semnificație selectat, numărul nD=fx și p = , , L; ) înmulțiți aceste valori ale rangului cu sx și deci determinați (k-l) rangurile cel mai puțin semnificative; ) se verifică semnificația diferenței dintre medii, plecând de la cele extreme din seria de clasare; comparați diferența dintre valorile maxime și minime ale mediei cu cel mai mic rang semnificativ la p = k, apoi găsiți diferența dintre media maximă și a doua medie din seria de clasare și comparați-o cu cel mai mic rang semnificativ la p = k- etc Continuați această comparație pentru al doilea cu valoarea mediei, care este comparată cu cea mai mică, și așa mai departe, până când diferențele dintre toate perechile k(k- / sunt examinate pentru semnificație Exemplul S-a studiat procesul de polimerizare radicalică a sărurilor pe bază de -viilpiridiu şi diferite halogenuri de alchil: CH I, C H I, C H I, SHCheI, C H B , C H Br, C H C , C H C Randamentul mediu de polimer pe opt replici pentru fiecare halogenură de alchil este prezentat mai jos: N" p / p Halogen alchili CH I C H I C HgI C H I C H Br C H Br C H CI CsHtCI In medie randament, % , , , , , , , , Eroarea de reproductibilitate în măsurarea randamentului polimerului sx este , Numărul de grade de libertate /x - În conformitate cu regula de aplicare a criteriului Duncan, aranjam rezultatele medii în ordine crescătoare: , , , , , , , , ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) Eroarea medie este sx - , / - , Să scriem rangurile semnificative din tabelul Duncan (Tabelul din apendice) pentru m - și un nivel de semnificație de , : R Ranguri , , , , , , , Cele mai puțin semnificative ranguri (LSR) înmulțite cu eroarea mediei sunt: p NZR • sx , , , , , , , ( ) - ( ) = , - , = , > , - diferență semnificativă ( ) - ( ) = , - , == , > , - diferență semnificativă ( ) - ( ) = , - , = , > , - diferență semnificativă ( ) - ( ) = , - , = , > , - diferență semnificativă ( ) - ( ) = , - , = , , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) - ( ) = , , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) - ( ) = , , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) ( ) = , , - diferență semnificativă ( ) - ( ) = , > , - diferență semnificativă ( ) - ( ) = , , — diferență semnificativă ( ) - ( ) - , > , - diferență semnificativă ( ) - ( ) = , go v Prin urmare, la un nivel de semnificație de p- , , valoarea gradului de extracție, egală cu %, trebuie considerată eronată, trebuie exclusă din eșantion și recalculată x și sx Comparație între distribuția eșantionului și distribuția populației Testarea ipotezelor privind parametrii de distribuție ai populației generale a fost efectuată sub ipoteza unei distribuții normale a variabilei aleatoare observate Ipoteza despre normalitatea distribuției studiate în statistica matematică se numește ipoteza principală Această ipoteză este testată pe un eșantion folosind teste de bunătate Criteriile de bunătate sunt utilizate pentru a testa ipoteza despre forma așteptată a legii distribuției Criteriile de bunăstare a potrivirii fac posibilă determinarea probabilității ca, conform unei legi de distribuție ipotetică, abaterea observată în eșantionul luat în considerare să fie cauzată de cauze aleatorii, și nu de o eroare în ipoteză Dacă această probabilitate este mare, atunci abaterea de la legea distribuției ipotetice ar trebui recunoscută ca aleatorie și ipoteza despre legea distribuției propusă nu trebuie respinsă Natura probabilistă a criteriilor nu permite acceptarea sau respingerea fără ambiguitate a ipotezei testate Criteriul ne permite să afirmăm că ipoteza nu contrazice datele experimentale dacă probabilitatea abaterii observate de la legea ipotetică este mare, sau că ipoteza nu este de acord cu datele experimentale dacă această probabilitate este mică Cel mai adesea, se folosește unul dintre cele două teste de bună adaptare: testul Pearson (testul X ) și testul Kolmogorov Pentru a aplica criteriul X (chi-pătrat) întregul interval de măsură- Variabila aleatoare din eșantionul de mărime n este împărțită în k intervale Numărul de intervale k este de obicei luat în funcție de dimensiunea eșantionului în intervalul de la la Numărul de intervale poate fi determinat printr-o formulă semi-empirică (IL ) Numărul de elemente eșantion care au intrat în intervalul i-lea va fi notat cu n Histograma construită (vezi Cap II ) a distribuției eșantionului sau considerațiile generale despre mecanismul de apariție a unei variabile aleatoare servesc ca bază pentru alegerea tipului de lege de distribuție Parametrii acestei legi pot fi determinați fie din considerente teoretice, fie prin găsirea estimărilor acestora dintr-un eșantion Pe baza legii de distribuție adoptată se calculează probabilitățile p, lovirea unei variabile aleatoare X în intervalul i-lea Valoarea care caracterizează abaterea distribuției probei de la cea așteptată este determinată de formulă * XI (nt-npt) ( ) la unde k este numărul de intervale; l este dimensiunea eșantionului Suma ( ) are o distribuție aproximativ X cu f = k-c- grade de libertate, unde c este numărul de parametri ai legii distribuției ipotetice determinate din eșantion Pentru o distribuție normală, c = dacă atât x, cât și sx sunt determinate din eșantionul dat Ipoteza despre tipul de lege de distribuție adoptat este acceptată la un nivel de semnificație dat p dacă ^ /), se ajunge la concluzia că ipoteza nu este în concordanţă cu distribuţia eşantionării Când se utilizează criteriul X , este de dorit ca dimensiunea eșantionului să fie suficient de mare: n - - , iar numărul de elemente n;> - Dacă oricare dintre u; * este frecvența relativă determinată prin formula (II ) Datele (cantitățile) necesare pentru calcul sunt date în tabel I хі Р* Xj Xj - X (Xj - xР (X, - x)z (X/ - X/ P'ț'Xj - XJ p*fxz - X) p^x, - X/ - , - , - , , - , - , , - , - , - , , - , - , , - , - , - , , - , - , , - , - , - , , - , - , , , - , - , , - , - , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Sume: - Ca urmare, obținem x - , microni, -Șț - , microni, P * - , microni, q * - microni Să definim coeficienții de asimetrie și curtoză conform formulelor ( ) și ( ): ? = — = - , , * = —— = — , și variațiile acestora - conform formulelor (A P ) și ( ): D (îî)= ( °-~ -= , ; /o( |') ® , ; ' ( + ) ( + ) " ѵ ' = - ( - ) ( ) = V I ' ( -g I) ( + ) ( + ) / D(i ) = , Astfel, avem z y "£> ( О) = , și I = , Să transformăm eșantionul într-o serie variațională Xi distribuția produsului /ko^ este aproape de o distribuție limitativă, pentru care au fost întocmite tabele Valorile critice pentru valoarea /ko sunt determinate din aceste tabele În tabel prezintă quantilele (n(o ), ^p Dacă valoarea calculată /ko este mai mică decât valoarea tabelată (/ko ), ,,, atunci ipoteza despre coincidența legii de distribuție teoretică F(x) Tabelul Cuantile de distribuție n'- p p , , , , , , , , o s , , , , , , , , , , , cu Fn(x) selectiv nu este respins Când ly >(nuF^-p, ipoteza este respinsă Nivelul de semnificație p este de obicei ales egal cu , Testul m este mai complet decât testul Pearson, folosește informațiile conținute în datele eșantionului În gruparea datelor care se produce la aplicarea testului Pearson, există o anumită grupare în sine duce la o oarecare pierdere a informațiilor conținute în eșantion În plus, distribuția u converge către legea limitei mult mai rapid decât X , mai ales în regiunea valorilor mari de u , care sunt esențiale doar pentru o estimare probabilistică Criteriul Wilcoxon Testul Wilcoxon este utilizat pentru a testa ipoteza că două eșantioane aparțin aceleiași populații generale Să existe mostre de variabile aleatoare X și Y de tip dimensiune Să transformăm mostrele în serii variaționale: X Ui și nX) numărul total de inversiuni și este distribuit aproximativ normal cu așteptările matematice ti - tp / ( ) și dispersie o = -^-(/u + n- ) ( ) La nivelul de semnificație p = , , conform ( ), valorile critice pentru ipoteza nulă vor fi h si ■ ti + , s„ ( ) Exemplul Una dintre cele mai comune metode de identificare a modelelor matematice ale proceselor de inginerie chimică este obținerea și prelucrarea curbelor de răspuns la aplicarea perturbațiilor de impuls pentru a obține așa-numitele curbe C A fost studiat modul hidrodinamic de funcționare al unui extractor industrial În acest scop, am folosit curba C a distribuției timpului de rezidență al particulelor în fază solidă în reactor, obținută de angajații NIUIF prin metoda pulsului la combinatul chimic Voskresensky din magazinul pentru producția de acid fosforic extractiv de către metoda dihidratului Obiectul studiului a fost un extractor cu secțiuni, la intrarea căruia a fost introdus un indicator în faza solidă (concentrat de apatită) Aurul radioactiv (g— μm) a fost folosit ca indicator Amestecarea se realizează cu mixere cu două niveluri (cu șase lame pe fiecare nivel) Din ultima secțiune, o parte din pulpă revine la prima Volumul primelor patru secțiuni = Vi = De la = Vi = m Volumul ultimelor patru secțiuni De la \u d Yv \u d Și? = De la - mz Curba C experimentală a fost luată la o sarcină pe celuloză Gn - , m /h și o valoare de reciclare R - , Timp de ședere în - tronsoane ti = ta = tz = Ta = I / bp - / , - , ore Timp de ședere în - tronsoane tz = te = r? \u d acelea - / , - , ore Curba C experimentală (Fig ), luată la ieșirea pulpei din extractor, a fost construită conform următoarelor date experimentale z, h С, mcr/mz , , , , , , , , , , , G, h C, mcr/mz , , , , , , , , , , , , , , Reshen și e Curba C experimentală a fost comparată cu soluția sisteme de ecuații diferențiale pentru un model de celulă cu reciclare Organizarea fluxului în reactor poate fi reprezentată sub forma unei diagrame bloc (Fig ) Sistemul de ecuații diferențiale pentru concentrația indicatorului în faza solidă are forma: = (si + Dc - ( + R)cJ - dt tj ■ - = ( + R) ci -( + R)c ] - dt tj de, „ ” = (( + R) c - ( + R) c ] » dt ~~~ = [( + R) C - ( + R) c*J • dt ti (! -Ș- \u d [(l + *) c - (l + R) c ] - de - = [( + R) c - ( + R) c,J - dt t ~~ = [( + #)ce - ( + #)c ] - • dt t - \u d (( + R) c-, - ( + R) c l - dt t Orez Curba C experimentală Orez Schema structurală a reactorului unde C este concentrația indicatorului, microni/m , - R este valoarea reciclării Condiții inițiale pentru sistem g md = — = —— -= , mcr/m ; t = vi s n- = - Sistemul de ecuații ( ) a fost rezolvat prin metoda Runge-Kutta pe calculatorul Minsk- Rezultatele obţinute sunt prezentate sub forma unui grafic (vezi Fig , curba ) G, h cu calcul, micro/mz , , , , , , , , , , , , , , Z, H Calculat, mcr/m , , , , , , , , , , , , , , Comparația a două curbe C: experimentală (y) și calculată (x) - a fost efectuată conform criteriului Wilcoxon A fost testată ipoteza H° de apartenență a două eșantioane (y) și (x) la aceeași populație generală Pentru a face acest lucru, elementele ambelor mostre au fost aranjate într-o secvență ascendentă comună: X , , X , Y , X , Y , X , Y , X , Z , X , x Y Y X Y X Y y Y Y X , , , , , , , , , , , Y Y X Y X Y X Y X Y x , , , , , , , , , , , Y x Y X Y x Y X Y X y , u i , I , , , , , , , X Y X Y Y X Y x Y XX , , , , , , , , Numărul total de inversiuni a fost determinat: „= - + + + + + + + - + + + + + + + + + + + + + + + + = Ipoteza I' nu este respinsă dacă numărul de inversiuni nu se încadrează în regiunea critică definită prin formulele ( ) Variabilă aleatorie și aproximativ distribuită normal cu așteptări matematice ty \u d tp / \u d - / \u d și standarde °și ("+"-!) *ZO -—( - - ) = , La un nivel de semnificație de p - , , valorile critice pentru ipoteza nulă, conform ( ), vor fi: „ - , - , ; „> , În consecință, numărul de inversiuni, egal cu , nu se încadrează în regiunea critică, iar diferența dintre curbele C comparate poate fi considerată nesemnificativă statistic Astfel, regimul hidrodinamic din fiecare secțiune a reactorului este aproape de amestecarea ideală Testarea ipotezei de normalitate pe un set de probe mici Să existe un număr n suficient de mare de eșantioane independente de aceeași mărime m Se impune testarea ipotezei de normalitate a populațiilor generale din care sunt prelevate probele, cu condiția ca parametrii acestor populații să poată avea valori diferite Luați în considerare abaterea relativă Hk t = „L ( ) unde хік este elementul i al eșantionului k; medie și rms abaterea tipică a probei k-a Se poate demonstra că distribuția valorii m nu depinde de parametrii populației generale m și o, ci depinde doar de dimensiunea eșantionului m Densitatea de probabilitate a valorii m este egală cu pentru |t| unde numărul de grade de libertate f \u d m - Din (IL ) pentru diferite valori ale lui m obținem: ( ) „i = , f(T) = Din (IL ) rezultă că atunci când /u = , abaterile relative din probele individuale sunt supuse unei distribuții uniforme dacă populațiile inițiale sunt normale Acesta poate fi folosit pentru testarea ipotezei de normalitate dacă numărul de probe este suficient de mare Cu m , din lipsa tabelelor necesare, se trece de la valoarea lui m la valoarea lui r): ■i—,, ■ ( ) V f + -m " Se poate dovedi că, sub populaţiile normale iniţiale, cantitatea m| are o distribuție Student cu f=m~ grade de libertate La testarea ipotezei de normalitate pe un număr mare de eșantioane mici, o valoare este selectată aleatoriu din fiecare probă O oarecare simplificare este posibilă aici - doar primele măsurători pot fi selectate, doar a doua, etc O astfel de selecție poate fi, de asemenea, considerată aleatorie Dacă numărul de elemente din eșantioane este mare, de exemplu, m > , atunci se pot face mai multe teste independente ale ipotezei, de exemplu, pe primul și ultimul element al fiecărui eșantion Atunci, dacă m = , pentru fiecare valoare selectată, conform formulei ( ), se calculează r, dacă m > , conform formulei (II ) q După trecerea la valorile ghidului pentru a testa ipoteza unei distribuții uniforme a lui r sau a unei distribuții a lui Student a lui q (și, prin urmare, normalitatea distribuției inițiale), oricare dintre testele de bunătate considerate anterior poate fi aplicat Exemplul Se cere testarea ipotezei distribuției normale a concentrației (g/l) de azotat de amoniu în aburul secundar după aparatul de reacție în producerea azotatului de amoniu conform rezultatelor a patru determinări în de probe ( tabelul de mai jos) Numărul eșantionului Хі Х х Х Numărul eșantionului Хі Х ХЗ Х , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Soluţie Luăm prima (хі) din rezultatele a patru determinări paralele ale fiecărei probe și calculăm valoarea lui m pentru fiecare dintre cele de valori folosind formula ( ) Rezultatele calculului sunt rezumate în tabel: Numărul eșantionului X sx XI - X m = ^ Numărul eșantionului X X Xi - X , , — , - , , , - , - , , , , , , , , , , , , - , , , - , - , , , , , , , - , - , , , , , , , , , , , - , - , , , - , - , , , , , , , , , , , - , - , , , - , - , , , - , - , , , , , , , - , - , , , , , , , - , - , , , , , , , , , , , , , , , - , - , , , - , - , , , , , , , - , , , , - , - , , , - , - , , , - , - , , , - , - , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , - , - , Din formula ( ) rezultă că la m = densitatea de distribuție a mărimii ? este egal cu / x) = -când |m| ) p - , (vezi Tabelul ) Valoarea calculată a lui lsh este mai mică decât valoarea tabelului Prin urmare, ipoteza distribuției normale a concentrației de azotat de amoniu în vaporii de suc nu este respinsă Capitolul descrie conceptele de bază ale statisticii matematice: populația generală și eșantionarea aleatorie, estimările și proprietățile acestora, metode de testare a ipotezelor statistice și de construire a intervalelor de încredere pentru așteptarea și varianța matematică Pentru obținerea estimărilor se folosește metoda maximă de probabilitate, care duce la obținerea unor estimări consistente, eficiente, deși uneori părtinitoare Testarea ipotezelor privind parametrii de distribuție ai populației generale se efectuează sub ipoteza unei distribuții normale a variabilei aleatoare observate Ipoteza normalității distribuției este testată folosind testele Pearson, Kolmogorov, Co , Wilcoxon și testele de bunătate a eșantionului mic Exerciții La fabricarea standardului cu un număr mare de măsurători, s-a stabilit în mod fiabil conținutul mediu al substanței, egal cu %, cu o eroare pătrată a unității care măsurare egală cu , % Determinați probabilitatea cu care se poate aștepta ca, la analize repetate, rezultatul mediu a nouă măsurători paralele să fie var rație în , h - , % Construiți un interval de încredere pentru așteptarea matematică a unei variabile aleatoare X distribuite normal cu o probabilitate de încredere de - , , dacă media eșantionului x- , se obține din patru măsurători, considerând dispersia egală cu , : a) generală; b) selectiv Estimați eroarea în determinarea densității unei substanțe folosind următoarele rezultate de măsurare: masă , g, eroare de măsurare a masei , g, volum , cm-s, eroare de volum , cmT Estimați eroarea de reproductibilitate sau alți factori privind variabilitatea mijloacelor este sarcina analizei varianței Analiza dispersiei folosește ceea ce s-a discutat în cap proprietatea de aditivitate a varianței variabilei aleatoare studiate, datorită acțiunii factorilor independenți R A Fisher în a definit pentru prima dată analiza varianței ca „separarea varianței atribuite unui grup de cauze de varianța atribuită altor grupuri” În funcție de numărul de surse de dispersie, se distinge analiza de varianță monofactorială și multifactorială Analiza varianței este deosebit de eficientă atunci când se examinează mai mulți factori În metoda clasică de cercetare, un singur factor este variat, iar restul sunt lăsați constant În același timp, pentru fiecare factor, se realizează propria sa serie de observații, care nu este utilizată în studiul altor factori În plus, cu această metodă de cercetare, nu este posibilă determinarea interacțiunii factorilor cu schimbarea lor simultană În analiza varianței, fiecare observație servește la evaluarea simultană a tuturor factorilor și a interacțiunilor acestora Analiza dispersiei constă în identificarea și evaluarea factorilor individuali care determină variabilitatea variabilei aleatoare studiate Pentru aceasta, varianța totală a eșantionului este descompusă în componente determinate de factori independenți Fiecare dintre aceste componente este o estimare a varianței populației Pentru a decide dacă influența unui anumit factor este semnificativă, este necesar să se evalueze semnificația variației eșantionului corespunzătoare în comparație cu variația de reproductibilitate datorată factorilor aleatori Semnificația estimărilor de varianță este verificată folosind testul Fisher (vezi Capitolul II, ) Dacă valoarea calculată a criteriului Fisher se dovedește a fi mai mică decât valoarea tabelată, atunci nu există niciun motiv să se considere semnificativă influența factorului luat în considerare Dacă valoarea calculată a criteriului Fisher se dovedește a fi mai mare decât valoarea tabelară, atunci factorul luat în considerare afectează variabilitatea mediilor În cele ce urmează, vom presupune că sunt îndeplinite următoarele ipoteze: ) erorile aleatorii ale observațiilor au o distribuție normală; ) factorii afectează doar modificarea valorilor medii, iar dispersia observațiilor rămâne constantă; experimentele sunt egale Cerința unei distribuții normale determină alegerea factorilor principali în studiul procesului prin metoda analizei varianței Dacă este necesar să se obțină o distribuție normală a valorii de ieșire, este de dorit să se atribuie aleatoriu doar acei factori a căror influență asupra valorii de ieșire este foarte mică O excepție poate fi făcută numai pentru acei factori care prin ei înșiși (din alt motiv) dau o distribuție normală a rezultatelor Factorii luați în considerare în analiza varianței sunt de două feluri: ) cu niveluri aleatorii și ) cu niveluri fixe În primul caz, se presupune că alegerea nivelurilor se face dintr-un set infinit de niveluri posibile și este însoțită de randomizare În același timp, rezultatele experimentului au mai multe valoare, întrucât concluziile experimentului pot fi extinse la întreaga populație Dacă toate nivelurile sunt alese aleatoriu, modelul matematic al experimentului se numește model cu niveluri aleatorii de factori (model aleatoriu) Când toate nivelurile sunt fixe, modelul se numește model de nivel de factor fix Când unii dintre factori sunt luați în considerare la niveluri fixe, iar nivelurile celorlalți sunt alese aleatoriu, modelul se numește model de tip mixt Uneori nu există nicio diferență în criteriile utilizate pentru diferite modele, iar singura diferență este generalitatea concluziilor, în alte cazuri există o diferență în criterii Analiza varianței poate fi aplicată sub diferite forme, în funcție de structura procesului studiat; alegerea formei adecvate este de obicei una dintre principalele dificultăţi în aplicarea practică a analizei Analiza unidirecțională a varianței Luați în considerare acțiunea unui singur factor A (cantitativ sau calitativ), care ia k valori diferite (niveluri ale factorilor) La nivelul z-lea, se fac m observații, ale căror rezultate pot fi scrise după cum urmează: Sus Ugі • • Ukg > Uіg Ugg • • • Ukg > Utg Ugpg ■ • ■ Ukp • Vom presupune că rezultatul oricărei observații poate fi reprezentat ca model unghi/ = n + dz+ (HI- ) unde u este efectul total în toate experimentele; rf, este efectul factorului A la nivelul i-lea (/ = , , k); — eroare de măsurare la nivelul z Să presupunem, de asemenea, că observațiile la un nivel fix al factorului sunt distribuite în mod normal în jurul mediei p + dj cu o varianță totală de o Numărul total de experimente este N: N = zzj + n + • • • + nk- (III ) Se testează ipoteza nulă de egalitate a valorilor medii la diferite niveluri ale factorului A: t = tg = = tk = t Cele mai simple calcule se obțin cu un număr egal de experimente la fiecare nivel al factorului A: = n = = pc = n (Tabelul ) În acest caz, numărul total de observații N este egal cu kn Notați cu y valoarea medie a observațiilor la nivelul i-lea P Tabelul Date inițiale pentru analiza unidirecțională a varianței cu un număr egal de repetări ale experimentelor Numărul de observație Nivelurile factorului A a \ ak Тп У Укі Ui! U Uk n ^ІL U p Ukp p p l Rezultate A'^yv A = £yv Ak = \UC /" Ui ■ și valoarea medie totală pentru întregul eșantion de N observații k p k == (III ) Pentru a efectua analiza varianței, variația totală a eșantionului s se descompun în componente care ar caracteriza contribuția factorului A și a factorului de șansă Factorul de aleatorie este ușor de evaluat datorită prezenței experimentelor repetate la fiecare nivel Definim varianța eșantionului la fiecare nivel: Dacă nu există încredere în acuratețea egală a experimentelor, omogenitatea dispersiilor s і , poate fi verificată prin testul Cochran (vezi cap , ) Dacă nu există diferențe semnificative între varianțe, varianța eșantionului este utilizată pentru a estima varianța generală c , care caracterizează factorul de aleatorie Numărul de grade de libertate ale varianței este k(n - ) = N - k O estimare aproximativă a varianței factorului A poate fi obținută după cum urmează: (III ) O estimare mai precisă pentru st poate fi obținută luând în considerare abaterile mediilor y la niveluri individuale de la media totală a întregului eșantion y Într-adevăr, De aici (W ) (W ) Varianța factorului A pentru modelul cu niveluri fixe (al ) nu este asociată cu nicio variabilă aleatoare, acesta este un nume condiționat pentru așteptarea matematică a pătratului mediu al abaterilor datorită influenței factorului A Această notație este convenabilă , deoarece definește dispersia cauzată de influența factorului A în mod similar cu factorul aleator al indicatorului de influență, care vă permite să comparați direct factorul A cu efectul aleatoriei De asemenea, introducem următoarea notație: k £ = y-, (Yi - Y ) « + • (III ) i=i Această dispersie are k - grade de libertate Dacă varianța este semnificativ diferită de ipoteza nulă nc = t = = tk = t este respinsă și influența factorului A este considerată semnificativă Ipoteza nulă este testată folosind criteriul Fisher Deoarece inegalitatea >cgH este o alternativă la = = , testul cu o singură coadă al lui Fisher este folosit pentru a testa ipoteza Influenţa factorului A este considerată semnificativă dacă s Și h=k(nl) = Nk (III ) ®ОШ Analiza varianței poate fi efectuată conform următorului algoritm: ) se calculează totalurile pentru coloane P A = ^Uiii •••> k' (III ) /=i ) sumele pătratelor tuturor observațiilor k pі = , , diferența dintre haloalchili ar trebui recunoscută ca fiind semnificativă După ce am stabilit cu ajutorul analizei de dispersie faptul că valorile medii ale randamentelor polimerului în ansamblu diferă semnificativ unele de altele, să trecem la o comparație a influenței haloalchililor individuali Să efectuăm această comparație conform testului lui Duncan (vezi Capitolul II, ) cu o probabilitate de încredere de - , Eroarea normalizată a mediei este s-~y , / = , y g Aranjați valorile medii în ordinea crescătoare a valorilor lor și scrieți din tabel aplicații ranguri semnificative pentru / - și p - , : s(C H, ) Yb (C H Br) s(C H, ) y (C,H Br) F (CH ) , , , , , R • • • Rangi, dl ■ ■ domnule Sy , , , , !/i - Kz = , > , - diferența este semnificativă Vx — Ul= , > , — diferența este semnificativă Wow — !/»= , , - diferența este semnificativă Uz - Ub~- , > , - diferenta este semnificativa Uz - Ug ~ , , - diferența este semnificativă Y r - Ya = , ,, b , ,bm Să presupunem că pentru fiecare combinație a nivelurilor factorilor A și B se efectuează n observații paralele (Tabelul ) Tabelul Date pentru analiza bidirecțională a varianței cu repetări B A Q а аі ак Ai ^ 'L'ng ■■■' J'l\P Y ' ^ • ••> ,G Іl Vі\ ' VІ\ -'Ubp Uk\y Uk\ •••> Uk-\ p V, b > Y ■••'Y\ P ,^ ' Y •••' Ucp UI ] ' Y Î UііP Uk \ ' Uk ■■■■ Uk%p bi Yi\ ■ v / ■■■• Ui)n Yi/Y ' Y I! ■■■' Y i/n Y kj\ Uk/ •••' Ukjn Bj bm V ,Uіt ■■■'Uіtl U W ' U TP •••' Uimn U ITP\ ' U ITP •••> U imn Uktp\' Ukpі ■■■' Uktp W Rezultatele I Ai Ak Numărul total de observații este N = pct Rezultatul observației poate fi reprezentat sub următorul model: = P + + + ( , ) unde c este media generală; a, — efectul factorului A la nivelul i-lea, / = , , ,k; fij - efectul factorului B la nivelul j-lea, j = , , ,li; a, B - efectul interacțiunii factorilor Efectul de interacțiune este abaterea mediei asupra observațiilor din seria (ij)-it de la suma primilor trei termeni din model ( ), iar EiJq (?= , , ,l) ia ținând cont de variația în cadrul seriei observații (eroarea de reproductibilitate) Vom presupune, ca și înainte, că este distribuit în mod normal cu așteptări matematice zero și varianță o sh Dacă presupunem că nu există interacțiune între factori, atunci putem adopta un model liniar: Уіі = й + аі + ₽; + Chі- (III ) Acest model este de obicei utilizat în absența observațiilor paralele (Tabelul ): Tabelul Date pentru analiza transversală bidirecțională fără replicare B A Rezultate au a ak Lu Y zi Ukі v, -^ Uіі Uk, Bz bm v m v' m vkt w togi Ai ■ k Luați în considerare mai întâi modelul liniar ( ) Notăm cu y, respectiv y) mediile pe coloane și rânduri: - La Ui \u d - t v U \u d - k prin y este media tuturor rezultatelor: (III ) (W ) = Y \u d - - km Uii (III ) Împrăștierea în medii peste coloanele y, y, , yk în raport cu media totală n nu depinde de factorul B, deoarece toate nivelurile factorului B sunt mediate Această împrăștiere se datorează influenței factorului A și a factorului aleator Deoarece varianța mediei (vezi Cap II, ) este de m ori mai mică decât varianța unei singure măsurători, avem La k G" (III ) T La rândul său, împrăștierea în mediile rândurilor nu depinde de factorul A și este asociată cu influența factorului B: m ^GG^(^ ^) i;av + T- (Sh- ) /= Egalitățile (ІІІ ЗЗ) și (ІІІ ) fac posibilă estimarea influenței facurei A și B, dacă se cunoaște estimarea varianței o Pentru a estima factorul aleatoriu în absența observațiilor paralele, procedăm după cum urmează Să găsim varianța observațiilor în /-a coloană m si = ^~F (III ) /= Această dispersie se datorează influenței factorului B și știri factorul a se intampla- Egalitatea va deveni mai precisă dacă în loc de s? variația medie ponderată pentru toate coloanele: utilizare ° ~ -uu-b (III ) Să notăm estimarea obținută (IIL ) pentru varianța o cu v UI Numărul de grade de libertate este k = TZT ~ ~yU ™l + , (ІН ) T /= Valorile lui z și i pot fi considerate ca variații ale eșantionului cu (k- ) și respectiv (m - ) grade de libertate Testez ipoteze tip glonț despre nesemnificația influenței factorilor A și B conform criteriului Fisher Dacă relaţia de dispersie F= - L-r (Sh ) DESPRE* influența factorului B este considerată semnificativă Când se testează ipotezele nule, se folosește testul Fisher unilateral, deoarece alternativa la egalitate = o w este inegalitatea analiza varianței în termenii unui model liniar ( ) este convenabil să se utilizeze următorul algoritm de calcul Găsiți: ) totaluri pentru golbtsam Dі= UC, Rezultatele calculului sunt prezentate convenabil sub forma unui tabel Tabelul Analiza bidirecțională a varianței (fără repetare experimente) Sursa de dispersie Grade de libertate Suma pătratelor Pătrat mediu Așteptarea pătratului mediu A k- ssA \u d ss - sst Sa k - taA + aosh În t - SSB = - SS SSB în m - l , kV + °osh Reziduu (Ar- ) (t- ) SS oct “ — SS% — SS ■+• Oct $ — — — - br (AD(ml) °SH kt - valoare totală ZZtot = - După ce a stabilit semnificația influenței acestui factor folosind analiza varianței, ei află apoi, folosind testul Student sau testul de rang Duncan, care valorile medii ale lui y sunt diferite Modelul liniar (III ) este valabil dacă nu există o interacțiune între factorii A și B În caz contrar, această interacțiune ca factor are propria sa varianță cjB Interacțiunea AB servește ca măsură a cât de mult depinde influența factorului A de nivelul factorului B și invers, cât de mult depinde influența factorului B de nivelul A În algoritmul de mai sus, dacă există o interacțiune între factorii ajB , ca parte integrantă, este inclusă în varianța în prezența observațiilor paralele Să fie efectuate n experimente paralele pentru fiecare combinație de niveluri ale factorilor A și B Deci, în tabel în celula formată din intersecția coloanei i-a și a rândului y, există o serie întreagă de observații Vіl> Vy - ■■■' Ucn- Să păstrăm notația din spatele rezultatului mediu din celulă Varianta eșantion de rezultate în fiecare celulă P = (PI ) are n- grade de libertate, celulele sunt omogene, pot fi ponderate variația medie Dacă variațiile eșantionului peste toate sunt mediate și rezultatul s (III ) ca o estimare a varianței de reproductibilitate a o Numărul de grade de libertate w este egal cu mk(n - ) O formulă mai convenabilă pentru calcularea varianței de reproductibilitate — dosh (III ) mk(n - ) k t p unde y!} este suma observațiilor din celula ij-a Când se efectuează ANOVA pentru un model neliniar, este convenabil să se utilizeze următorul algoritm de calcul Conform tabelului găsiți: ) sumele observațiilor din fiecare celulă Uiii - ^^uni' (III ) u= \u d , m\ i \u d , , , k; ) pătratul sumelor observațiilor din fiecare celulă Uii Wu Ni (III ) ) totalurile coloanei t p At= uni-, = și= (II ) ) totaluri de rând (W ) ) suma tuturor observațiilor (total general) p k k t (III ) ) suma pătratelor tuturor observatii k t p SS Wiii' (III ) ) suma pătratelor totalurilor coloanei împărțită la numărul de observații din coloană, k (III ) ) suma totalurilor de rânduri pătrate împărțită la numărul de observații pe rând, t SS = (III ) ) pătratul totalului mare împărțit la numărul tuturor observațiilor (termen de corecție), / k t p k V ss = mkn mkn ) suma pătratelor pentru o coloană = SS - SS ; ) suma pătratelor pentru un rând SSfi = SS - SS ; (III ) (III ) ) suma pătratelor pentru varianța reproductibilității k t a ssoH = sS — ) suma totală a pătratelor, egală cu diferența dintre suma pătratelor tuturor observațiilor și termenul de corecție, tot = SS - ss ; ) suma reziduală a abaterilor pătrate pentru efectul de interacțiune A B (III ) (III ) SSAB - SSo ui~SS^~SSB~SSoiu' (III ) N ) dispersie s\ - /(* - Oi (III ) ) discurs s^b sb = SSB/(m~ ); (III ) ) Varianta ^^ LV (III ) ~ (k - ) (m - ) ) dispersia reproductibilității SS osh (III ) °w ~ mk(n- ) Testarea ipotezei despre semnificația interacțiunii factorilor A și B se efectuează conform criteriului F în același mod pentru modelele cu niveluri aleatoare și fixe Cu toate acestea, testarea ipotezelor despre semnificația factorilor A și B este efectuată diferit pentru diferite modele În tabel Figura prezintă o analiză bidirecțională a varianței cu execuții repetate pentru un model cu niveluri aleatorii Tabelul Analiza bidirecțională a varianței pentru modelul cu niveluri aleatorii (cu experimente repetate) Sursa de dispersie Grade de libertate Suma pătratelor Pătrat mediu Așteptarea pătratului mediu A k- SSA \u d SS - SS , \u d A * - | Fri°A + P°AW + Posh V t - SSg - SS - SS SSB SR - o t - PaAV ^ ° osh AB (k- ) (t- ) S$AB = ® ' - SS ~ SS + s SSAB AB Fi-p (ti• fi)> (III ) unde p este nivelul de semnificație; f\*=k- ; fi = (kl)(m - ) În mod similar, influența factorului B este considerată semnificativă dacă $v/ SAB > Fi-p (fi' ft)' (III a) unde / i \u d m - , / \u d (/ c - ) (t - ) Dacă inegalitățile (III ) și (ІІІ a) nu sunt satisfăcute, influența factorilor A și B ar trebui considerată nesemnificativă Pentru un model matematic cu niveluri fixe, termenii corespunzători interacțiunii dispar din sumele abaterilor pătrate SSA și SSB Ca rezultat, pentru a evalua semnificația factorului A, o relație de dispersie a formei (PH ) al cărui numitor este estimarea variaţiei de reproductibilitate Raportul de dispersie rezultat este comparat cu tabelul f p(f, ff) F \u d sb / osh> (III a) care se compară cu tabelul Fx ț(f\, fi) pentru numerele de grade de libertate /i = m - și / = mk (l - ) Dacă rapoartele de dispersie (PI ) și (III a) sunt mai mari decât cele tabelare SV soiu >Fi-p(fi' fi) și (III ) SO U> Fl-P (fl•fi)' influența factorilor A și B trebuie considerată semnificativă Dacă inegalitățile (PI b) nu sunt îndeplinite, influența factorilor A și B este nesemnificativă Pentru a verifica semnificația efectului de interacțiune, o relație de dispersie a formei și comparați-l cu tabelul (f, fi) la nivelul de semnificație p și numărul de grade de libertate f = (k- )(t - ) și fi^mk (l- ) Dacă raportul de dispersie rezultat este mai mare decât tabelul %w/%w > F(fi, fi), efectul interacțiunii factorilor trebuie considerat semnificativ În caz contrar, dacă Fx p(f, fi), influența efectului este reciprocă acțiunile trebuie considerate nesemnificative Exemplul S-a studiat influenţa asupra procesului de sinteză organică a doi factori; A este tipul de solvent la nivelurile ai, az, az, al și B este tipul de haloalchil la nivelurile li, bz, bz, bt Rezultatele (randamentul polimerului în procente) sunt prezentate în tabel: B A au " az al b, , , , , , , , , b , , , , , , , , bz , , , , , , , , Z» , , , , , , , , Pentru fiecare combinație de tip de solvent și haloalchil, s-au făcut două experimente paralele Este necesar să se evalueze semnificația influenței tipului de solvent și haloalchil asupra procesului de sinteză Soluţie Modelul matematic al experimentului este un model cu niveluri fixe Nivelurile factorilor A și B nu au fost alese aleatoriu, deoarece este necesar să se stabilească efectul asupra procesului de sinteză al acestor patru tipuri de solvenți și haloalchili Calculul se efectuează în conformitate cu algoritmul de mai sus, conform formulelor ( ) - ( ): Determinați sumele observațiilor din fiecare celulă (tabel) B A Rezultate ai ag az iadul h , , , , , bb , , , , , bz , , n o , Ll , , , , Rezultate , , , , , Să punem la pătrat sumele rezultate Prezentăm rezultatele lui yfj sub forma unui tabel: B A ai ag az iadul b' , , , , /> , , , , bz , , , , Ba , , , Calculați totalurile pentru coloane De exemplu = , + , + , + , = , Să calculăm totalurile pe linii De exemplu, B \u d , + , + , + , \u d , Să definim totalul general - suma tuturor observațiilor: ' uo "= n * = in / \u d і - Să determinăm suma pătratelor tuturor observațiilor; SSi= î /^ = , Să determinăm suma pătratelor din coloană, totalurile coloanei împărțite la numărul de observații |r \u d - - ( , + , + , + , ) SS =- - -= , Determinați suma pătratelor totalurilor rândului împărțită la numărul de observații din rând, ss = \u d yi-^ ( , + , l + , + , ) \u d , = , Determinați pătratul totalului mare împărțit la numărul tuturor observațiilor ss = , , — = -^- = ' - Să determinăm sumele abaterilor pătrate pentru factorii Liy: SSX = SS - SS = , - , = , , SSfl = SS - SS = , - , = , Determinați suma pătratelor pentru varianța de reproductibilitate: k t este suficient de mare Deci, există de pătrate latine X , de pătrate latine X Planificarea experimentului conform schemei pătratului latin este utilizată atunci când se studiază influența a trei factori A, B și C asupra procesului În acest caz, factorii A și B pot fi asociați cu studiul în sine și eterogenitatea materialul este considerat factor C Toți cei trei factori din pătratul latin au același număr de niveluri (a/r b-, ct) Deci, în plan (Tabelul ), fiecare factor se modifică la două niveluri pătrat latin suprapus X Matricea de planificare - tabelul corespunzător , planul experimental, care include trei coloane și patru linii, este prezentat în tabel Pătratul latin face parte din plan - conform schemei pătratului latin, în planificare este introdus al treilea factor C Cu toate acestea, întregul plan (Tabelul ) este de obicei numit pătrat latin În pătratul latin, fiecare element se repetă o singură dată în fiecare rând și în fiecare coloană, deci, indiferent de proprietățile perturbatoare ale elementului pătratului, acestea vor afecta în egală măsură calculul mediilor pe coloană și pe rând Date în tabel Planul este o replică pe jumătate a PFE (Tabelul ) Experimentele incluse în semi-replică sunt marcate cu asteriscuri Rezultatul observației obținute prin experimentul factorial complet poate fi reprezentat ca următorul model: Ui ;, • (III ) În tabel prezintă planul experimentului după schema pătratului latin X Tabelul Pătrat latin X A B Rezultate b\ bb bz a\ C S'i C U Cs Ui A, U С ^ СЗ З'Б С Ув Аъ az SZ U C Uv C ue Az Rezultă „ Vz Din punct de vedere structural, pătratul latin X poate fi considerat ca o replică de / a experimentului factorial complet În cazul general, pătratul latin și X și poate fi considerat ca / și o replică a PFE u Când se efectuează analiza varianței pătratului latin fără experimente repetate, este convenabil să se utilizeze următorul algoritm de calcul Pentru a face acest lucru, determinați: ) totaluri pentru rândurile A, -, coloanele Bj și literele latine C De exemplu, pentru cele date în tabel pătrate latine X rânduri totale: Ai \u d Ui + Ui + Uz > A \u d + uv, A \u d y ; + y + y , totalurile coloanei: - Ui + Ui + Ui • - U + Uv + U "\u e B \u d Uz + Uv + U" =i /=i - ' ) suma pătratelor totalurilor rândului împărțită la numărul de observații din rând, SS = ( , ) ) suma pătratelor totalurilor coloanei împărțită la numărul de observații din coloană, ss = ( , ) ) suma pătratelor totalurilor pentru literele latine, împărțită la numărul de observații corespunzător fiecărei litere, P SS = ( , ) ) pătratul totalului general, (termen corectiv), /p\ q=\ împărțit la numărul tuturor observațiilor ss = n ( , ) P I \ S / V -d V / l osh ( ) V / V , - diferența este semnificativă ya - Vі = , - , = , > , - diferența este semnificativă y - ~y = , - , = , , - diferență semnificativă ~Yi = , - , = , C d Y bx C d Y? ay bz cs dz Uz b Cі d Uz />, C dș U bz C d\ U /> Cz d} U Tabelul A B C= C= C= C= C= D = = = = = £= £= £= £= £= F= £= F= £=• £= C \u d C \u d C \u d C \u d C \u d D = = = = = £ = £ = £ = £ = £ = £= F = £= £= F= C= C= C= C= C= = = = = = £ = £ = £ = £ = £ = £ = £= £= £= £= C= C= C= C= C= = = = = = £ = £ = £ = £ = £ = £ = £ = £= F= F= C \u d C \u d C \u d C \u d C \u d = = = = = £= £= £= £= £= F= £ = F= £= £= În pătratul greco-latin, există GR diferite combinații de niveluri de factori în loc de n combinații ale experimentului complet cu patru factori Prin urmare, pătratul greco-latin este o replică /rfi dintr-un experiment factorial complet (FFE) Deci, dat în tabel Pătrat greco-latin X este o replică / a PFE (N= ), pătratul greco-latin X - / este o replică a PFE (N= ), X este o replică / de PFE (N= ) Analiza varianței pătratului greco-latin se realizează în același mod ca și analiza pătratului latin obișnuit, ținând cont de al patrulea factor D (litera greacă) Suma pătratelor pentru o literă grecească are n- grade de libertate Numărul de grade de libertate al sumei reziduale, definit, ca mai înainte, ca diferență între suma totală a pătratelor și sumele pătratelor tuturor factorilor, este egal cu (n- )(n- ) Dacă suprapunem trei pătrate latine ortogonale unul peste altul, obținem un pătrat latin de ordinul al treilea, n pătrate ortogonale - un pătrat latin de ordinul al n-lea Pătratele rezultate se mai numesc și pătrate hiper-greco-latine Cu n niveluri, n + factori pot fi introduși în plan Numărul de grade de libertate al sumei reziduale va fi atunci egal cu zero Astfel de planuri se numesc saturate Să construim un plan saturat pentru n = Pentru a face acest lucru, suprapunem cele patru pătrate latine ortogonale obținute X [vezi (III ) - (III )], constituind o serie completă de pătrate latine ortogonale X (Tabelul ) Pătratul latin original (III ) corespunde nivelurilor factorului C, al doilea pătrat (III ) corespunde nivelurilor factorului /), etc Nivelurile factorilor sunt indicate prin numere Designul experimental corespunzător pentru șase factori este prezentat în tabel Planul rezultat este saturat, deoarece numărul de grade de libertate al sumei reziduale, determinat de formula / = (u- ) (u k + ), unde k este numărul de factori studiati, este egal la zero Planul este o replică / a PFE Astfel de planuri sunt de obicei utilizate în etapele incipiente ale cercetării procesului, când Tabelul Proiectare experimentală și = , N= Numărul experimentului А В С D Е F Numărul experimentului А В С D Е F - Q Orez cub latin de ordinul întâi Este necesar să se efectueze o enumerare complexă a factorilor calitativi pentru a evidenția combinațiile promițătoare pentru cercetări ulterioare și pentru a le elimina pe cele inacceptabile Folosirea pătratelor greco-latine și hiper-greco-latine ca planuri experimentale economisește atât numărul de observații, cât și duce la o simplificare a calculelor Principala ipoteză care stă la baza aplicării pătratului greco-latin și pătrate de ordine superioară, este ipoteza absenței interacțiunilor între factori Verificați adecvarea model liniar acceptat, ca în aplicarea pătratelor latine, posibilă numai în prezenţa unor experimente paralele cuburi latine Un experiment factorial complet pentru trei factori „ fn> ) corespunde unui aranjament cubic de n elemente, inclusiv „ poziții Trei muchii ale cubului corespund factorilor A, B și C cu nivelurile , , , , n- (Fig ) Dacă introducem al patrulea factor D în plan și plasăm nivelurile acestui factor ( , , , , n - ) în punctele aranjamentului cubic corespunzător experimentelor, atunci obținem un cub latin de marime si prima comanda Un cub latin de mărimea n de ordinul întâi este un tabel cubic de n elemente situate în poziții rfl, în care fiecare element intră rfi ori și apare în fiecare dintre cele de planuri paralele cu planurile de coordonate xyoxi, x\oxs, xroxs, la fel pentru toate elementele și egal cu n număr de ori Într-adevăr, nivelurile factorului suplimentar D (elementele cubului latin) apar în plan același și egal cu r de ori și apar în fiecare dintre cele planuri de coordonate (adică, cu nivelurile celor trei factori A, B, C) același și egal cu n număr de ori (Tabelul ) X X cub latin Tabelul A B Matricea de planificare corespunzătoare pentru un cub latin cu dimensiunile n = , r = este dată în tabel Proiectarea unui experiment pe cubul latin de ordinul întâi vă permite să includeți patru factori (A, B, C și D) în considerarea dvs Diferența față de pătratul greco-latin, care face posibilă și izu- Tabelul Plan experimental l = , N = Numărul experienței А В С D У Numărul experimentului А В С D У Yi Y U wib Uz U U U U U U U O yi Y U U U U -a U yi I Y U U Uіz U U Influența a patru factori este că în cubul latin trei factori (A, B și C) sunt considerați principali și un factor (D) constituie grupul eliminator, iar în pătratul greco-latin sunt considerați doi factori A și B cele principale, iar C viD constituie un grup cu dublă eliminare Numărul de experimente într-un cub este de n ori mai mare decât într-un pătrat greco-latin Cubul latin fără experimente repetate se aplică sub ipoteza unui model de proces liniar: ВіІцІ - Iх + аі + ?Z + 'tq + Z + zijqi • (III ) unde și este media generală; az este efectul factorului A la nivelul z-lea, / = , , , n - ; J j■ - efectul factorului B la nivelul y-lea, y - , , , , n - ; la? - efectul factorului C la nivelul q-lea, q = , , , , n - ; , efectul factorului D la nivelul l-lea; este o eroare experimentală aleatorie Analiza statistică a cubului latin de ordinul întâi fără experimente repetate este efectuată în mod convenabil conform următorului algoritm Defini; ) totaluri pentru toți factorii de la fiecare nivel; Ai(i = , , , , n - ), = , n - ), Cq(q = , , n- ), Dz(/ = , , n- ) În ceea ce privește, de exemplu, planul prezentat în tabel , avem D) \u d Ui "Ui + Uz 'I / io " Un " / + /u " Hugo " / > A \u d Ui + Yb - Uts "Uіz + / - Uіb - Uc + Ug; i + / \u e Ag \u d Uі - Uv "Uz - Un - Up + Uiv " Ug "Uiv " Ug?> Vo \u d Ui - Ui "Ui - / io " / іz “ Uiv - Uiv - Un + Ug, Ві \u d Y "Y " Y I "Y and + UN " Y n "Y go ~ Ugz " Ugv\u e Vya \u d Uz - Uv - Uz - / + / ib - Uiv - / ~ / - Uc\u e \u d Ui - Ui Uz - Ui - Ub - U in - Ui - Uv + Uya Ci = i/ + i/ + U + U + U + / + U + Uiv + U T O = Ui + U° + U« + Uch + Uch + Uch + U» + U + Uiv • Di - Ug + Uv + Ui + / + £/іz + Yn + Uch + U + U i (/ = , ) - X ( = , , ), X Z(( = , , ); Tabelul Numărul experimentului Xi X xs X l • u MPa yi • , MPa Uz, % D , , , , , , , , , , , , , , NU , , , , , , , , , , Tabelul Rezultate pentru diferite niveluri de factori Răspunsuri Aditivi (хі) Cantitatea de aditiv Cantitatea de umplutură ("") (xs) Y - În continuare, au fost determinate sumele pătratelor pentru toate sursele de dispersie • ) SSXi = SS - SSe; de exemplu, pentru tine SSX = - = ; ) SSX = SSs - SSe; de exemplu, pentru tine, SSXJ = - = ; ) SSX =SS -SS ; de exemplu pentru SSX= - - ; ) SS, = SSs - SSe; de exemplu, pentru tine, SSX = - = ; ) suma totală a pătratelor este egală cu diferența dintre suma pătratelor tuturor observațiilor și termenul corectiv: gen=SSi ~sse; de exemplu pentru bsch \u d - \u d ; ) se folosește suma reziduală a pătratelor pentru a estima eroarea experimentală SS CT = bw - (SSXI + SSXj + SSXj - SSX , — diferența este semnificativă y ^ — y^ = - , - diferența este semnificativă /Y ) - p( ) , ( ) - /Y °> = , > , - diferența este semnificativă p , - diferenta semnificativa Nivele de factor X , Valori medii/) , , , , , , , , , Clasamente, g , , , , , , , , rXsD , , , , , , , , р , — diferență semnificativă e( ) p( ) o > , — diferență semnificativă — p( ) , > , — diferență semnificativă D( ) - p(°) = , , - diferență semnificativă — JD( > — , - o( > = , = , = , , - diferență semnificativă p( ) - p = , - /U ) = , , — diferență semnificativă — Z/ ’ = , = , , - diferență semnificativă O* ' - O* ' = , x: P cov = TȚ Ui (Хі ~ ~ este mai mic de % Coeficientul de corelație al eșantionului r*, precum și r este coeficientul de corelație al populației generale, în valoare absolută nu depășește unitatea: - (IV ) Cu un număr mic de experimente și o corelație relativ mare, distribuția coeficientului de corelație diferă semnificativ de cea normală (Fig , a) Pentru a construi un interval de încredere, puteți utiliza transformarea Fisher; de aici z* = al z-lea = e z - e + (IV ) (IV ) Distribuția lui z este aproape neschimbată ca formă, cu variația r* și u și cu creșterea n se apropie rapid de normal (Fig , ) cu o medie egală cu Orez Densitatea distribuției coeficientului de corelație al eșantionului t,= - Іu g (IV ) si cu standardul (IV ) valoarea necunoscutului t (IV ) U Și ' U G=Uz Atunci cu încredere probabilitatea fi este înăuntru „ r — — — cu condiția excluderii influenței lui xr asupra y În notația pentru coeficientul de corelație parțială, acest factor exclus este plasat în index după punct Când se studiază dependența lui y de trei factori xi, x și xs, coeficientul de corelație parțială dintre y și x; cu condiția ca x și xs să fie constante, pot fi calculate prin formula r * uh, x, x, GUI Z G» G o-'izG'O-'s zK (IV ) În trecerea de la coeficienții de corelație de pereche la cei privați, nu numai valoarea coeficientului de corelație, ci și semnul se poate schimba semnificativ Să ilustrăm acest lucru cu un exemplu S-a studiat viteza de coroziune (K) a probelor de oțel care conțin sulf (S), fosfor (P) și cupru (Cu) într-o soluție de acid citric Pe baza unui eșantion de de experimente s-au obținut valorile coeficienților de corelație de pereche; r^s = + , ; r'p = + , ; + , ; ^scu ~+ , ; gKCu = - , ; gr C =+ , Folosind formula (ГV ), găsim coeficienții de corelație parțială, excluzând influența unuia dintre factori; rK S-Cu ~ + , , rKCu-s= , ; gc Cu R~— , rKS-p = - , ; rs p Cu = + , ; rK Р-cu „" " , , ^Pcu-S~ , ; /KP S = + , ; Cu p = + , ; O comparație a valorilor coeficienților de pereche și de corelație parțială arată că efectul, de exemplu, al fosforului asupra vitezei de coroziune la un conținut constant de cupru este mai mare decât la unul variabil, iar efectul fosforului asupra vitezei de coroziune la un conținut constant de sulf este mai mic decât la unul variabil: * * • 'TcP S \p SCII = - , ; rK Cu SP = - , Coeficientul de corelație de pereche dintre viteza de coroziune și conținutul de fosfor la concentrații variabile de cupru și sulf este pozitiv (r*KR = + , ); coeficient de corelație parțială fa SCu = - , Astfel, analiza corelației a permis stabilirea naturii și gradului de influență a cantității de sulf, fosfor și cupru conținute în oțel asupra vitezei de coroziune a acestuia într-o soluție de acid citric În cazul general, pentru a calcula coeficienții de corelație parțială, puteți utiliza matricea de corelație eșantion; Coeficientul de corelație parțială dintre x, - și y este determinat de formula Ay (IV ) „ (xn) astfel încât probabilitatea de acest eveniment P este maxim Probabilitatea p ca variabila aleatoare Y să cadă în intervalul y! - e/ , y! + e/ în prima aproximare, este egal cu: Pi = e exp I ti/;—? (hoi ! • Să găsim probabilitatea ca sistemul de variabile aleatoare independente Ui, U , ,U„ să ia un set de valori din intervalele Ui -e/ , y,+e/ , i = , , , ,n: P / P \u d gP P 'r \u d - exp - I I * "- * (X') g me a ( { P - Evident, pentru un anumit o , probabilitatea maximă P va primi - Acesta este cazul când suma din exponentul exponentului este minimă, adică P [«/i - ? (х/)] * * = tip = Folosind metoda celor mai mici pătrate, orice date experimentale pot fi prelucrate, dar optimitatea acestei proceduri este dovedită doar pentru o distribuție normală În acest caz, se poate vorbi de statistici suficiente, adică de astfel de funcții din rezultatele observațiilor (estimări pentru parametrii populației generale), cu ajutorul cărora sunt extrase toate informațiile despre acești parametri conținute în experiment Sarcina de a determina coeficienții ecuației de regresie folosind metoda celor mai mici pătrate se reduce practic la determinarea minimului unei funcții a mai multor variabile Dacă y = /(x, & , &j, bx, bk) (IV ) este o funcţie care este diferenţiabilă astfel încât P f \u d S mai mult \u d o condiţie necesară pentru minimizarea egalităţilor df f —= , = dbp db sau P [ Ui — f(xi " • , bk) df (*i) dbk Sistemul de ecuații (IV ) conține atâtea ecuații cât numărul de coeficienți necunoscuți b , bu b , ,bk este inclus în ecuația de regresie și se numește sistem de ecuații normale în statistica matematică Funcția Φ> pentru orice bj, bv b , ,bk, prin urmare, trebuie să aibă în mod necesar cel puțin un minim Prin urmare, dacă sistemul de ecuații normale are o soluție unică, atunci acesta este minimul pentru funcția Ф Când se studiază dependența de un parametru variabil, este util să se construiască o linie de regresie empirică pentru a determina tipul de ecuație de regresie Pentru a face acest lucru, întregul interval de modificări x în câmpul de corelare (Fig ) este împărțit în k intervale egale Dx Toate punctele care se încadrează în acest interval Dx sunt raportate la mijlocul lui xy Pentru aceasta se calculează mediile parțiale pentru fiecare interval pi Șwi Ui =—— , (IV ) unde n; - numărul de puncte din intervalul Dx-; în care k %n, = n- (IV ) /=i l este dimensiunea eșantionului Apoi conectați secvențial punctele (x;, y) cu segmente de linie dreaptă Linia întreruptă rezultată se numește linia de regresie empirică y pentru x După forma dreptei de regresie empirică, puteți alege ecuația de regresie y \u d f (x) Regresia liniară de la un parametru Este necesar să se determine coeficienții ecuației de regresie liniară folosind metoda celor mai mici pătrate l Y \u d b pentru o probă de volum n Sistemul are forma o + M (IV ) ecuații normale ~ (& + b Xi) = °' — f="l p p %YiXi - + b *i = °' = = p p na+bz ^ = ^' i=i (IV ) Wad este numărul de grade de libertate ale dispersiei adecvării; (IV ) / este numărul de coeficienți din ecuația de regresie; Svosp ~ suma pătratelor asociată cu variația de reproductibilitate i = lU= - m' Ui=—'U Uii- (IV ) și=\ SBocnp ~ ^play//play> (IV ) Losp este numărul de grade de libertate ale dispersiei de reproductibilitate; /play - (ті О "(IV ) r-M n ti ®play ■ “, (IV ) (ti - ) >= S SM, - suma reziduală a pătratelor; "ti l "oem \u d ^ (Uіi-Ui) ', = ȘI= P f rest ~ t = (IV ) ventilator este numărul de grade de libertate ale dispersiei reziduale s C m SSqct SOCT~f ma opresc ic-^)a Dacă F = s / s r\d' play (IV ) (IV ) se dovedește a fi mai mică decât valoarea tabelară F p(f\, fit) pentru nivelul de semnificație p și numărul de grade de libertate /i și / = y e ref, ecuația este adecvată experimentului Pentru același număr de experimente m, = ms = = m, = = mn = m, calculele sunt simplificate: s ad - t (Ui - Ui) - • (IV ) n - p t c VOSpr \u d (Uіi-Uі) și - - (iv ) n (t ~ ) Dacă experimentele au fost efectuate fără paralele și pentru a obține dispersia de reproductibilitate s-au efectuat o serie separată și tone de experimente, atunci Și " L ( ' = De îndată ce Na+ rem încetează să fie semnificativ mai mic decât s k CT, creșterea gradului k ar trebui oprită Semnificația diferenței dintre s și j +) este verificată prin testul lui Fisher: F = &sk+i • Dacă raportul E obținut este mai mic decât F tabelar (fi, fi) pentru nivelul de semnificație selectat și numerele de grade de libertate /і =fk și fi =fk+i, creșterea gradului de k trebuie oprită Polinoame Cebyshev Ecuația de regresie exprimată în termeni de polinoame Chebyshev are forma y = bllP (x) + blP (x)+ +bhPk(x), (IV ) unde P (x),>Pi(xfi ,Pk(x) sunt polinoame Chebyshev ortogonale pe mulțimea punctelor xt, X , ,xn Aceasta înseamnă că pentru toate j relațiile = (IV ) i=i unde Pk+\(x) depinde doar de mărimea eșantionului n Cunoscând polinoamele Cebyșev Pk+i(x), cu fiecare creștere a gradului ecuației de regresie, este necesar să se calculeze doar coeficientul bk+I Polinoamele Chebyshev sunt determinate de formule P (x)= , (IV ) l - (IV- ) fe (P k ) Pk+ (fi = (x) Є (x) - Pk-i (fi- (IV ) De exemplu, P (x) = x - (l - ) x + ~ ” + (IV ) O r, , (n + ) n + n + (l -f- ) (l + ) (l + ) P (x) \u d Xi ' - X -L- - X - v - - (IV ) Ptw= + '|,+ u x -a±daі+d±d> (l + ) (l + ) (l - ) (l - ) (IV ) Determinând coeficienții b , b\, ,bk ai ecuației de regresie (IV ) folosind metoda celor mai mici pătrate, obținem P I= calculați prin formulele (IV ) nu depind de care va fi ordinea ecuației de regresie care se determină La găsirea ecuației de regresie prin metoda rafinărilor succesive, se folosesc toate bj găsite anterior Creșterea ordinii ecuației de regresie cu duce la definirea unui singur coeficient În același timp, formulele pentru calcularea varianței reziduale pentru ecuația de regresie de ordinul k sunt obținute în mod convenabil: rest = ~ ' (, V' ) unde sumele abaterilor pătrate , ' $' sunt determinate de formula recurentă SSk = SSA , - bk V (Xi) (IV ) = Este necesar doar să precalculați SS : ss" \u d -b ° po (xi)j \u d (yi - \u d y\ ~ i= i= i= P Cu valori echidistante ale argumentului X \u d X + L; X =X + L; ; xn \u d Xj + (n - ) A, unde h este pasul de interpolare, calculele coeficienților sunt simplificate Să facem o schimbare de variabile: z = ^-^+ (IV ) h Apoi fiecare valoare x, va fi înlocuită cu numărul ei, adică z, = i Să determinăm coeficienții ecuației de regresie a formei Y = „aP =i (IV ) yt рі w d= p w i=i (IV ) pupa= Уі Pk (i) t=l ■ î p k(') i=l (IV ) Sumele din numitor pot fi determinate prin formula prescurtată: >=i (A!) n(n - )(n - ) (n - k ) [( * - ) P} ft ( k + ) (IV ) unde ( k - D-'-' este produsul tuturor numerelor impare de la la k~ inclusiv În special, = , (IV ) l (l - ) (l - ) l l (l - ) (l - ) (l - ) ha (L - ) (L - ) (L - ) (L - ) (IV ) (IV ) (IV ) Aceste sume sunt, de asemenea, utilizate pentru a calcula sumele SSk necesare pentru a determina varianța reziduală: n o(z) + „i Pi (z)\ P (r)= ; Pr(z) = z~& Conform formulei (IV ) ov \u d / C \u d , ; l U Pi(i) conform formulei prescurtate (IV ) este egal cu = lum - L = până în ' r=i iar coeficientul ai conform (ГV ) este egal cu a, = / = , Ecuația de regresie de ordinul întâi este y = , + , (r - ) Să determinăm varianța reziduală ^Rem Pentru a face acest lucru, folosind formulele (IV ) și (IV ), calculăm sumele: ( Ulі ,=! „ oo = z-, Y = = , * Și unsprezece SS = SSe— oh? V R? ( = - , - = , - = Unde Definim acum ecuația de regresie de ordinul doi: Y \u d "" ( ) + a A (r) + " ^ (r) Prin formula (IV ) Pi(z) este egal cu P (r) = r - r + Prin formula (IV ) găsim unsprezece ^ și se determină „g - / = , Ca rezultat, obținem o ecuație de regresie de ordinul doi: L ^i-p(/i / ), prin urmare, ecuația de ordinul doi este o rafinare esențială a ecuației de ordinul întâi Să definim ecuația de regresie de ordinul trei: L Y \u d "L (r) - (r) - " ^ (r) + "spz (r) ■ Înlocuind n - în formula (IV ), obținem Ps (r) = r - r + , g- , Prin formula (IV ) găsim unsprezece ( - ) ( - ) ( - ) = , , Unde „ \u d , / , \u d , , iar ecuația de regresie de ordinul trei este y \u d , + , (z - ) + , (z - z + ) + + , (z + z + z - , ) Să verificăm semnificația tranziției de la ecuația de regresie de ordinul doi la ecuația de ordinul trei Pentru aceasta definim SS " , - , - , , s oct = - = - = -~ -= T= , s > s oprire L oprire ' Prin urmare, trebuie să ne concentrăm pe regresia de ordinul doi Înlocuind y și r cu a și I, obținem (i - , ) = , + , (- + , și, în sfârșit i \u d , - , / + , f + Ecuația rezultată este cea finală din clasa polinoamelor Să estimăm etanșeitatea conexiunii găsite folosind relația de corelație (IV ): c \u d y e, e \u d SS / SS \u d , / \u d , , O \u d / - , \u d , Raportul de corelație Ѳ este apropiat de , prin urmare, relația găsită este apropiată de strict funcțională Regresia transcendentală Pentru dimensiunile mici ale eșantionului n, o creștere de ordinul polinomului poate duce la o creștere a varianței reziduale Pentru a reduce numărul de coeficienți determinați, se utilizează regresia transcendentală Calcularea coeficienților de regresie transcendentală poate fi foarte consumatoare de timp din cauza necesității de a rezolva un sistem de ecuații neliniare Calculul este simplificat printr-o schimbare de variabile De exemplu, dependențe de tip exponențial și putere fracțională =V* (IV ) = bou' (IV ) sunt liniarizate luând un logaritm: lg y \u d lg bn + X lg (IV ) POZITIV „O=lg” (IV ) )gy=z, lg bv = a U = gx=/ obținem ecuații liniare cu privire la noile variabile: Coeficienții ai, a, b sunt determinați prin metoda celor mai mici pătrate Pe baza a și a obținute se determină coeficienții b și by Cu toate acestea, trebuie avut în vedere faptul că coeficienții ecuațiilor de regresie (IV ) și (IV ) obținuți în acest fel sunt estimări părtinitoare pentru coeficienții generali corespunzători Estimarea etanșeității conexiunii neliniare Dacă presupunem că ecuația de regresie este găsită cu suficientă acuratețe, atunci varianța reziduală se datorează numai prezenței varianței de reproductibilitate, adică ~ stop ~ play Cu cât este mai mică proporția conexiunii dintre Y și conexiune Prin urmare forța Ser în varianța totală sj, cu atât mai mult X, deoarece proporția mai mică a aleatoriei conexiunii poate fi caracterizată prin valoarea (IV ) mai puternică în asta Unde ("- despre sL "= — (l - ) S* (IV ) ost p - i h = La Și n - eu P Ui — „=I Y = - P Conexiunea este mai puternică, cu cât valoarea este mai mică (IV ) se numește corelație Cu cât Ѳ este mai mare, cu atât conexiunea este mai puternică este o suprafață de regresie pentru k \u d și o suprafață pentru k> Această suprafață se numește suprafață de răspuns La construirea suprafeței de răspuns, sunt reprezentate grafic valorile numerice ale parametrilor (factorilor) pe axele de coordonate ale spațiului factorilor Materialul statistic inițial este prezentat în tabelul Tabelul Sursa de material statistic la „scara de viață” scala la una nouă prin glisare normalizat Trecând de la natural ku a tuturor valorilor variabilelor aleatoare conform formulelor: o Ui~ U y' = ~ - i= , , , l; /= , despre XH~XJ P~ s k (IV ) unde yfi x$ - valorile normalizate ale factorilor relevanți; y, x, - valori medii ale factorilor; sy, ^—abateri standard ale factorilor: În tabel prezintă materialul statistic original la o nouă scară: Tabelul Sursă de material statistic la scară adimensională Numărul de experiență la L \ * x? * xb x * Hk Ug ѵ * * * hіz Uz p ѵ° Xtn * X% „ Xk „ U ° „ Pe noua scară avem: x° = , Y = (IV ) Și ^= > V = E Coeficientul de corelație al eșantionului în acest caz este egal cu P ' : i g , = eu xt P - O O x, xt h L - (IV ) P oh oh /, m = , , Coeficientul de corelație al eșantionului calculat prin formula (IV ) este egal cu coeficientul de corelație dintre variabilele exprimate la scară naturală r * x , r? Xt Ecuația de regresie dintre variabilele normalizate nu are termen liber și ia forma -"' aproximativ y° \u d a x + o x + - - + i—l (= i=l = " - *і + " ( x /) ^ + "A Jj AiXki = *ЧУі> (IV ) = = i=li=l “i Y + “ T *lі* і - h “l T ( ) = S xmui = = = = Înmulțim laturile stânga și dreapta ale sistemului de ecuații (IV ) cu /(n - D Ca urmare, pentru fiecare coeficient o, obținem, conform (IV ), coeficientul de corelație al eșantionului r * Ținând cont de faptul că P obţinem un sistem de ecuaţii normale sub forma „ + + Vzd + • • • + akrxtxk = , „ * +“a+“z * + ••• + a»r'x^k = ' (IV ) „ rxLx, + a rxftx, + „ rxh*, - - ak = rft’ În sistemul de ecuații (IV ) r*/xm=^*x/- Pentru procesele multiparametrice, sistemul (IV ) este de ordin înalt și este necesar să se folosească un calculator pentru a-l rezolva După rezolvarea sistemului (IV ), se calculează coeficientul de corelație multiplă R: R - r ^ ary xy + a ^ ry , + ■ • • + (IV ) Coeficientul de corelație multiplă servește ca măsură a puterii relației pentru regresia multiplă: o , unde yo este media pe experimente paralele Numărul de grade de libertate Uvosp este egal cu Dispersia reziduală se determină prin formula (IV ): Y (y i=lt = li=li=l L R nn ^ Х ІХ І + bl Xli " • • " bk xiixki ~ Х іУі I = = = = (IV ) P nn n * xkix i + * xktxii + ■ • • + bk хіі = ^іхкіУі-i=li=li=li = l Sub formă de matrice, sistemul de ecuații normale poate fi scris după cum urmează: Xt XB = Xt Y (IV ) Într-adevăr, înmulțind matricea Xm și X, avem xthx = p p p хоіХ І • • • xoixkl = = = nnn Х ІХ І X}ț • • xiixki i=li=li=l (IV ) nnn xkixoi xkîxii * X^ i=li=li=f XTX este matricea momentului Înmulțind matricea Xm X cu matricea coloanei B, obținem matricea coloanei: HTXV = p p p " •*' " ^ -^ki t==l Z= z- nnn Xji - bk xnxt i=ii=ii=i (IV ) nnn bt xktxoi - * xkixii - • • • ~*fe *йі i=li=l = Înmulțind matricea Xm cu vectorul de observație Y, obținem Xm Y = P hoiu I= P chiUi i=i (IV ) P хкіУі i= Din ecuația (IV ), matricea coloanei coeficienților B se determină după cum urmează: B = (Xr X)' Xm Y, (IV ) unde (X X) este matricea inversă matricei (ХТХ); (x ^) " \u d coo coi • • • coL C ° c" C k (IV ) -cho cki ••• ckk Elementele matricei inverse sunt determinate de relație sii - - ' det (Xt X) (IV ) unde det (XTX) este determinantul matricei X X iar a este complementul algebric al elementului £ xn,x,, din matricea XTX Pentru existența matricei inverse (XTX) trebuie să fie nedegenerată În acest sens, atunci când se utilizează metoda de calcul luată în considerare, este necesar ca variabilele Xp x^ ,x,t să fie liniar independente Apoi, în matricea variabilelor independente, elementele unei coloane nu vor fi o combinație liniară a elementelor corespunzătoare ale altor coloane Pentru a determina varianța reziduală, se determină o matrice de coloană (IV ) Numărătorul varianței reziduale se obține prin înmulțirea matricei: (IV ) Notăm cu B vectorul coloană al coeficienților adevărați de regresie, în timp ce așteptarea matematică a lui B este egală cu M(B) = R T atunci : prin cel mai bun polinom este statistic nesemnificativ diferită de eroarea de reproductibilitate Pentru a verifica acuratețea aproximării, calculăm erorile rădăcină pătratică medie / pe mulțimea Mv; [Ui - ZJ (OF , /= , , , lts (IV ) Găsiți cel mai mic dintre toate cele și notați Dacă se dovedește că este statistic nesemnificativ diferit de eroarea de reproductibilitate, atunci procedura de aproximare este finalizată Dacă este semnificativ mai mare decât eroarea de reproductibilitate, este necesară a treia etapă de aproximare, la care polinoamele z}, z , , zm, vor fi folosite ca argumente fictive În absența informațiilor despre eroarea de reproductibilitate, valoarea poate fi utilizată pentru a încheia procesul de aproximare (IV ) Dacă E este relativ mic, de exemplu, mai mic de , - , , atunci aproximarea ulterioară este ineficientă De remarcat că în multe cazuri valoarea lui crește odată cu numărul de etape de aproximare Acest fenomen se numește instabilitate MGUA Instabilitatea GMDH se manifestă și printr-o modificare a valorii de *, care mai întâi scade odată cu creșterea numărului de etape și apoi începe să crească rapid Valoarea s* este mai puțin sensibilă la creșterea numărului de etape decât , deoarece punctele aparținând setului de testare Mn participă la formarea polinomului de aproximare Se recomandă finalizarea procedurii de aproximare succesivă la etapa corespunzătoare valorii minime a lui s* Multe aspecte legate de aplicarea și justificarea metodei de analiză de grup a argumentelor trebuie încă dezvoltate în detaliu Un avantaj important al acestei metode este posibilitatea în fiecare etapă de aproximare de a face față unei matrice informaționale de ordin care nu depășește X , care în cazul general îmbunătățește condiționalitatea acesteia Metoda componentelor principale Metoda componentelor principale a fost propusă de Pearson ( ) și dezvoltată mai târziu în detaliu de Hotelling ( ) Metoda face posibilă trecerea de la factorii măsurați direct Xj( =l, , k) la combinațiile lor liniare necorelate: k vi= ^UjiXi, i, j = , , , k (IV ) ale căror varianțe scad, adică Coeficienții combinațiilor liniare, încărcările variabilei i-a în componenta j-a și /(- sunt elemente ale vectorilor proprii ai matricei de covarianță, iar variațiile componentelor sunt egale cu valorile proprii ale acestei matrice Fie X o matrice de observații centrate cu zero înseamnă vector: x Xfo * • »• xh% / * xJ Xlo Buze x%n•••xkn Apoi matricea N N xm x L \u d - * - - - \u d (n- ) - n - Х ІХ І N хііхкі (IV ) NN xMx*t ••• 'ki N N * ×*' matricea de dispersie-covarianță a observațiilor; elementele sale diagonale sunt varianțe, elementele în afara diagonalei sunt covarianțe Dacă trecem de la variabile centrate la cele normalizate (vezi IV ), obținem o matrice de corelație; Gxx • ■ ■ rxtxh rx X ■ ■ ■ GC X r" r* ••• » xkxl xkx' (IV ) Pentru orice matrice simetrică L (sau/î), există o matrice ortogonală L astfel încât /*i UT LU -* = I° \ (IV ) unde Ăr (r - , , , /c) sunt valorile proprii ale matricei Z, iar U este o matrice ortogonală în care a r-a coloană este al-lea vector propriu corespunzător celei-a valori proprii Suspin matricea / sunt definite ca rădăcinile polinomului caracteristic (A -Xt) = , (IV ) unde E este matricea de identitate; I L-XE i \u d (- ) * X * + (- ••• + m (IV ) Primul, al doilea și ultimii termeni ai polinomului -caracteristic al matricei L sunt definiți simplu: /l = tr(L) este urma matricei L (suma elementelor diagonale ale acesteia); Pk = \ L I Este mult mai dificil să se determine membrii rămași ai polinomului caracteristic Se propun următoarele formule pentru coeficienții polinomului caracteristic: Lt \u d L Pi \u d tr (D - Li - pkE - L = Ytf % == ^* (IV ) ^k- - LBk- Bk-i ~ fhi Pk-jE Pk-l — k t* (*-fe- Lk - Pk = - tr (Lk) Bk = Lk - PkE unde E este matricea de identitate Din punct de vedere geometric, găsirea componentelor principale se reduce la o tranziție la un nou sistem de coordonate ortogonale Se determină prima axă de coordonate astfel încât combinația liniară corespunzătoare să extragă cea mai mare varianță posibilă, apoi se găsește o axă ortogonală cu aceasta, care face același lucru cu varianța rămasă etc Din noile coordonate, puteți reveni la cele vechi : k Xi= Y UfiVf, i, i= , , , k, (IV ) /= unde Vj este j-a componentă principală; uJt este masa componentei j din variabila a i-a Proporția de varianță, exprimată în %, explicată prin componenta a r-a este definită după cum urmează: X = -*■ % (IV ) ^ Într-o serie de probleme, se pot trage concluzii interesante dintr-un număr mic de componente În aceste cazuri, metoda componentelor principale face posibilă reducerea dimensiunii problemei datorită faptului că combinațiile liniare cu variații mici sunt eliminate și sunt luate în considerare numai combinațiile liniare cu variații mari De exemplu, un proces este caracterizat nu de o valoare de ieșire, ci de un întreg set de indicatori corelați Este de dorit să știm care combinațiile de indicatori explică o mare parte a varianței Acest lucru va determina ce să studieze în continuare Componentele principale pot fi considerate ca un nou set de măsurători care rezultă dintr-o combinație liniară a măsurătorilor originale Folosind metoda componentelor principale, se poate încerca să reducă un număr mare de caracteristici de proces înalt corelate la una sau două caracteristici sau componente generalizate, ceea ce va facilita evaluarea calității procesului la fiecare punct de observație Toate acestea vor avea sens numai dacă combinațiilor liniare rezultate li se poate da o interpretare fizică rezonabilă Metoda componentelor principale este, de asemenea, utilizată ca o modalitate de ortogonalizare a matricei de variabile independente Dacă V este matricea valorilor componentelor principale obținute prin transformare ortogonală: V^XU, (IV ) atunci noile variabile μ sunt independente și matricea lor de covarianță este diagonală: Vt V=UT Xt XU=(n—l)UT LU (IV ) Ținând cont de (IV ), avem IrV = (n - ) X (IV ) Dacă Y este un vector de valori ale mărimii de ieșire și toate condițiile preliminare ale analizei de regresie sunt îndeplinite, este posibil să se obțină estimări ale coeficienților ecuației de regresie liniară Y-VB, (IV ) unde B este vectorul coeficienților În conformitate cu (IV ) obținem V = (n - ) - * £ - * VG, (IV ) sau pentru fiecare coeficient - " bj \u d (n - l HI; ^uUi (IV ) i=i evaluarea semnificației coeficienților de regresie pentru componentele principale se efectuează conform testului T Student destul de corect, deoarece toate covarianțele dintre coeficienți sunt egale cu zero: = II II s^j ^(u - ) X; $play Se poate arăta că coeficienții de regresie pentru variabilele independente x și pentru componentele principale sunt legați după cum urmează: (IV ) * bj= Ujtatl , (IV ) unde a, - coeficient de regresie pentru variabila x,, / - , , , k-, bj - coeficient de regresie pentru componenta Vj; w„ sunt coordonatele J-ro ale vectorului propriu; /, este un factor de scară Componentele principale nu sunt invariante la schimbările în scara acelor scale pe care sunt numărate variabilele De obicei lucrează cu variabile normalizate, reciproca abaterii standard /, - = este luată ca factor de scară, adică toate calculele sunt efectuate pentru matricea de corelație Metoda componentelor principale a fost aplicată cu succes pentru a construi un model matematic al procesului industrial de flotare a minereului de potasiu în condițiile uzinei Beloruskali În timpul experimentului pasiv la uzină, au fost înregistrate valorile a de parametri de proces de intrare și doi de ieșire Au fost construite ecuațiile de regresie pentru componentele principale pentru indicatorii de proces yt și y , s-a estimat semnificația coeficienților conform formulei (IV ) Conform ecuațiilor de regresie și încărcărilor componente, s-a evaluat gradul de influență al variabilelor asupra indicatorilor de proces y, u , s-au tras concluzii despre relația procesului și s-au dat recomandări pentru îmbunătățirea acestuia Capitolul este dedicat metodelor de analiză de corelare și regresie, care sunt utilizate pe scară largă în prelucrarea rezultatelor observațiilor Modelul matematic al procesului este reprezentat de un polinom ai cărui coeficienți sunt determinați prin metoda celor mai mici pătrate Metoda celor mai mici pătrate este fundamentată ca un caz special al metodei probabilității maxime cu o distribuție normală a variabilei aleatoare observate Modelul polinom este foarte convenabil, deoarece vă permite să îmbunătățiți aproximarea prin creșterea ordinii polinomului, duce la un sistem liniar de ecuații normale la determinarea coeficienților ecuației de regresie prin metoda celor mai mici pătrate Sunt luați în considerare diferiți algoritmi de analiză de regresie pentru procesarea unui experiment pasiv Corelația dintre coeficienții ecuației de regresie obținute prin prelucrarea experimentului pasiv complică analiza statistică și interpretarea rezultatelor Metodele de experiment activ evidențiate în capitolul următor fac posibilă depășirea acestor neajunsuri ale analizei regresiei clasice Exerciții Tabelul prezintă caracteristicile membranelor omogene schimbătoare de anioni obținute în diferite condiții de nitrare a copolimerului: хі — temperatură, °С; X - durata, h; xs este gradul de nitrare, %; yi - rezistenţa electrică volumetrică specifică, Ohm ■ cm: y > - rezistenţa la tracţiune, MPa Numărul experimentului Xi X xa l ■ uch ■ Numărul experimentului Xi X xs J'i • Y ■ Determinați coeficienții de pereche și corelația parțială între parametrii procesului Determinați dependența conținutului de Fe, % (y), în cristalele de sulfat de cupru CuSO - H O de conținutul de FeSO , g/l (x), în lichidul mamă Fiecare experiență se repetă de două ori X Y X Y , , , , , , , , , , , , a) Evaluează omogenitatea dispersiilor; b) determinarea varianţei de reproductibilitate; c) alegeți tipul de dependență funcțională y \u d f (x)\ d) determinați ecuația de regresie; e) efectuează analiza de regresie și corelație a rezultatelor Determinați dependența rezistenței hidraulice a stratului de umplutură Dr (/i = m, umplutura inelului Raschig X X , mm) de viteza de curgere fictivă uCp folosind polinoame ortogonale Chebyshev: csr , cm/s Dr- , MPa , , , , , , , , , , La obținerea fosforului prin sublimare din fosfori de calciu s-a studiat dependența gradului de reducere a fosfatului (y) de temperatura (x) pentru fosforitul Karatau Experimentele au fost repetate de - ori X, °C % X, °C y % , , , , , , , , , , , , , , , a) Verificați omogenitatea variațiilor eșantionului; b) se calculează coeficientul de corelație dintre temperatură și gradul de reducere; c) determinați coeficienții ecuației de regresie; d) evaluați semnificația coeficienților și adecvarea ecuației de regresie la experiment ) S-a investigat dependența gradului de oxidare (y) al cromitului CrO la cromat CrO de durata calcinării (x) a încărcăturii la °C Fiecare experiment a fost repetat de două ori ordine), prin metode de algebră liniară pentru a determina ecuația de regresie și pentru a efectua analiza de regresie și corelație a rezultatelor Partea Metode de proiectare a experimentelor CAPITOLUL V METODE DE PLANIFICAREA EXPERIMENTELOR EXTREME Un mare număr de probleme experimentale din chimie și tehnologia chimică sunt formulate ca probleme extreme: determinarea condițiilor optime de proces, a compoziției optime a compoziției etc Datorită dispunerii optime a punctelor în spațiul factorilor și transformării liniare a coordonatelor, este posibil să se depășească deficiențele analizei clasice de regresie, în special, corelația dintre coeficienții ecuației de regresie Alegerea planului experimental este determinată de formularea problemei de cercetare și de caracteristicile obiectului Procesul de cercetare este de obicei împărțit în etape separate Informațiile obținute după fiecare etapă determină strategia ulterioară a experimentului Astfel, apare posibilitatea controlului optim al experimentului Proiectarea unui experiment vă permite să variați toți factorii simultan și să obțineți estimări cantitative ale principalelor efecte și efecte de interacțiune Efectele de interes pentru cercetător sunt determinate cu o eroare mai mică decât în cazul metodelor tradiționale de cercetare În cele din urmă, utilizarea metodelor de planificare crește semnificativ eficiența experimentului Experiment factorial complet La planificarea conform schemei unui experiment factorial complet (FFE), toate combinațiile posibile de factori sunt implementate la toate nivelurile alese pentru studiu Numărul necesar de experimente N în timpul PFE este determinat de formulă W = l* unde n este numărul de niveluri; k este numărul de factori Dacă experimentele sunt efectuate numai la două niveluri, cu două valori ale factorilor și, în același timp, toate combinațiile posibile de k factori sunt efectuate în timpul experimentului, atunci configurarea experimentelor conform unui astfel de plan este numit un experiment factorial complet de tip k Nivelurile factorilor reprezintă limitele zonei de studiu pentru un parametru tehnologic dat De exemplu, se studiază efectul a trei factori asupra randamentului produsului (y, %): temperatura (zi) în intervalul – °C, presiunea (z ) – – Pa și timpul de rezidență (z ) ) – de minute Nivelul superior din punct de vedere al temperaturii z,max este egal cu °C, nivelul inferior zj”’” este egal cu °C Apoi, pentru z, avem max g tip max -tip O „ Г і —гі \ \u d \u d ; D = J -= iad În general, pentru orice factor zt max g type Aproximativ / „G / / = k-, (V ) max tip / — / ' L • (V ) Punctul cu coordonatele (z?, ) se numește, uneori, centrul planului se numește nivelul principal; Az, este intervalul de variație de-a lungul axei zz Din variabilele z,, zk se trece la cele noi - x,, xk prin următoarea transformare liniară; zj - g°, ; /=b , , k (V ) Pentru variabilele xp xk nivelul superior este egal cu + , nivelul inferior este , coordonatele centrului planului sunt egale cu zero și coincid cu originea coordonatelor În exemplul luat în considerare, k = Numărul de combinații posibile N a trei factori la două niveluri este V = k = = Planul de realizare a experimentelor (matricea de planificare) este scris sub forma unui tabel Prezentat în tabel planul la o scară adimensională poate fi interpretat geometric ca opt vârfuri de cub (Fig ) Să introducem în PFE o coloană a așa-numitei variabile dummy ho = (Tabelul ) Date în tabel matricea de planificare are următoarele proprietăți; N х“ *; / * O; N ^ = ^; / = o, i, *, i= unde k este numărul de factori independenți; N este numărul de experimente din matricea de planificare Prima proprietate (V ), dispariția produselor interioare a tuturor vectorilor coloană, se numește proprietatea de ortogonalitate a matricei de planificare Această proprietate reduce dramatic dificultățile asociate cu calcularea coeficienților ecuației de regresie, deoarece matricea coeficienților ecuațiilor normale (XTX) devine diagonală și ei Tabelul Experimentul factorial complet T? Numărul de experiență Factori la scară naturală Factori la scară adimensională Z Z XI Х Y - - - ' + - - - + - + + - - - + + - + - + + + + + elementele sunt egale cu numărul de experimente din matricea de planificare diagonală niya N Coeficienții ecuației de regresie cu cele mai mici pătrate sunt determinați după cum urmează; = (XrX)~lXTY &o bt bt Ț ^ • Matricea momentelor (X AG), corespunzătoare tabelului arată ca — ~ Х І *ІІ *• « *•* Х,і*зі х І i=l B = - Luând în considerare proprietățile (V ), obținem (V ) Matricea inversă matricei momentelor se dovedește a fi egală cu Prin urmare, orice coeficient al ecuației de regresie bj este determinat de produsul scalar al coloanei y și al coloanei corespunzătoare Xj, împărțit la numărul de experimente din matricea de planificare N; N Folosind planul prezentat în tabel , calculați mai întâi coeficienții ecuației de regresie liniară y = bd + bxxx + bgxi + ft x (V ) De exemplu, pentru a determina coeficientul, obțineți suma produselor: la x, este necesar La - " - - + + - - + X = + - - + + - - + - + xcui = x\>y, bl= În mod similar, obținem /u = , ; /> = - , ; />h = + , Dacă introducem în considerare o ecuație de regresie mai completă cu coeficienți de interacțiune Y \u d b + b X + b x + b x + bl xlx + bi xlx + b x x + b ( z+++z (V ) apoi pentru a determina coeficienții bm, bіz, bіz (efectele interacțiunii perechi) și b, z (efectul interacțiunii triple), este necesară extinderea matricei (Tabelul ) după cum urmează (Tabelul ) Tabelul Matrice de planificare extinsă pentru experimentul factorial complet E Numărul experimentului X XI X Xs X X XіXs X XE X X X Y + - - - + + + - + + - - - - + + + - + - - + - + + + + - + - - - + - - + + - - + + + - + - + - - + - + + - - + - + + + + + + + + Efectele de interacțiune sunt definite în mod similar efectelor liniare Deci, pentru a determina coeficientul este necesar: (*iXih Uі \u d - i \u d i “ i(ХіХі)іУі Restul coeficienților sunt definiți în felul următor: £>іz — , * £> з , , £> — , Dacă se realizează experimente paralele suplimentare, este posibil să se determine ssr, să se verifice semnificația coeficienților de regresie și, în prezența gradelor de libertate, adecvarea ecuației Datorită faptului că matricea de covarianță (ХТХ)~' pentru un experiment planificat este o matrice diagonală G /LG ~ (xthx) - = W (V ) coeficienții ecuației de regresie nu sunt corelați între ei Semnificația coeficienților ecuației de regresie poate fi verificată pentru fiecare coeficient separat folosind testul t al lui Student Excluderea din ecuația de regresie (VII) a unui coeficient nesemnificativ nu va afecta coeficienții rămași În acest caz, coeficienții eșantionului £> se dovedesc a fi așa-numitele estimări nemixte pentru coeficienții teoretici corespunzători ; £>;-> îb (V ) adică valorile coeficienților ecuației de regresie caracterizează contribuția factorului corespunzător la valoarea lui y Elementele diagonale ale matricei de covarianță sunt egale între ele, deci toți coeficienții ecuațiilor (V ) și (V ) sunt determinați cu aceeași precizie: sbj ~ *repro/ (V ) De exemplu, trei experimente paralele suplimentare au fost plasate în centrul planului și s-au obținut următoarele valori ale lui y: ; S"! / = ^-= , ; Sft și ) sunt nesemnificative și ar trebui excluse din ecuație După eliminarea coeficienților nesemnificativi, ecuația de regresie are forma y \u d , "T , * i" T , xd - , ■ Să verificăm caracterul adecvat al ecuației obținute după criteriul Fisher: g / l r - °odihnește/ joacă ' l » s „OST N—l - = ; secnp - O' ®' / - numărul de coeficienți semnificativi din ecuația de regresie, egal cu Atunci F- / , = , A)- , , F „*■ ? + ?із> bz -*■ ?з + ?і > (V ) unde P sunt așteptările matematice pentru coeficienții corespunzători Acești coeficienți generali nu pot fi estimați separat conform unui plan care include doar patru experimente (Tabelul ), deoarece cu În acest caz, coloanele pentru termenii liniari și produsele de pereche sunt aceleași Dacă, de exemplu, pe lângă coloanele date în tabel , calculați o altă coloană pentru produsul x, x , rezultă că elementele acestei coloane sunt exact egale cu elementele coloanei x Astfel, reducerea numărului de experimente duce la estimări mixte pentru coeficienți Pentru a determina ce coeficienți generali sunt amestecați, este convenabil să folosiți această tehnică: punând x în loc de x}x (Tabelul ), obținem relația X = X X , (V > numit raport generator Să înmulțim ambele părți ale raportului generator cu x : în timp ce în stânga obținem o singură coloană: / = x x x (V ) Produsul (V ) se numește contrast definitoriu, cu ajutorul acestuia este convenabil să se determine ce coloane au aceleași elemente Înmulțind pe rând contrastul determinant cu x , x^, obținem x, = x? x x = x x ; x = XjX ; Xs - XiX (V: ) Relaţiile obţinute (V ) corespund sistemului de estimări mixte (V ) Atunci când utilizați TFE, este necesar să aveți o idee clară despre așa-numita rezoluție a replicii fracționale, adică să determinați în prealabil care coeficienți sunt estimări nemixte pentru coeficienții generali corespunzători Apoi, în funcție de sarcină, se selectează o replică fracționată, cu ajutorul căreia se poate extrage maximum de informații din experiment De exemplu, într-o problemă cu patru factori, k = , se poate lua ca relație generatoare х =хіх хз (V ) și oricare dintre produsele pereche de factori, de exemplu x =x x (V ) Matricea de planificare cu raportul generator (V ) este dată în Tabel Tabelul Semi-replică din PFE cu raport de generare (V ) Numărul experimentului XO Xi X Xs X Numărul experimentului XO Xi X Xs X + + + + + + + + - - + - - + + + - - - - + - + + - + - + - + + + - + - + + - - + Folosind contrastul definitoriu /=x x x x , obținem următorul sistem de estimări comune pentru coeficienții ecuației de regresie: * = * * xit Xi = x xaxt, * = X * xi' Xi = * X, x , * X = * * * * = *a* * X = X, X , -*■ ₽ + ₽ > ^ "* ? + Ș " ₽ * ? > ->■ tU + ? ^ „*■ ₽X " ? i ^ „*■ ? + ? i ^ „*■ ? + ? - (V ) În problemele reale, interacțiunile triple sunt egale cu zero mult mai des decât cele duble Dacă estimările pentru efectele liniare sunt de cel mai mare interes, trebuie luată relația generatoare x =x]x x Cu raportul de regenerare (V ), matricea de planificare are forma (Tabelul ) Tabelul Numărul experimentului XO X X xs X Numărul experimentului XO X X xe X + + + + + + + + - + + - - + + + - „ - + + - + + - + - + - - + + - + - + + - - - Contrastul definitoriu este exprimat prin raportul I = x, x x Rezultă următorul sistem de rating: Xj \u d xs x , Xj = Xj X , * = * X, x x , X = Xj ha, * * = * * * * * = * * * * * = * * * bl „► ? + ? bl -*■ p + ? > ba “*■ ?з * ?i > bl -*■ ₽ + Ріа ^ ? + ? bia -*■ ? z ?іz ^ „*• ? + ? (V ) Prin urmare, este logic să folosiți o replică fracțională cu un raport generator x = x}x dacă coeficienții , p și sunt de cel mai mare interes În semi-replica din PFE cu raportul de generare x = x x x x , toate efectele liniare și efectele interacțiunii perechilor sunt amestecate numai cu efecte de ordin triplu și superior: * - * * * * * "bd - * % " V" * \u d * * * * , &i - *? -? S Xj = XjX X XS, „*■?“+ P ”“ x = Xj Xj X X", ^ "*■ ? + ? • X = Xj Xj Xz X , bl ₽ + Ș S- x = * x x X , b, "*■ ?S + ? J > * * = * x x", bc "*• ? + ? " > XjXs = X X X , bts -*■ ? + ₽ S> X X = X X XS, bts : „*■ ₽n+ Șj S> Xj x = X X X , bi„ „*■ + ? > x = * x x , bas ? + ₽ > Xj X = Xj x X , bc + ? S" Xj X = Xj X X bc ? S + ? > * = * Xj x , bsi „»■ ? + ₽ J X XS= XjXjX , dar „*■ ? b+ • (V ) Neglijând efectele de interacțiune mai mari decât ordinul doi, putem presupune practic că pentru k > semireplicile PFE oferă estimări nemixte pentru efectele de interacțiune liniare și de pereche De asemenea, în practica experimentală sunt utilizate replici UD, replici '/v etc O replică fracţională, în care p efectele liniare sunt echivalate cu efectele de interacţiune, este notă Pentru un sfert de replica, de exemplu, în planificarea pentru k= tip , pot fi setate relații de generare: X =XjXjX > X = XjXj Contrastele definitorii pentru această replică vor fi rapoartele / \u d x xx x , / \u d x xgx (V ) Înmulțind contrastele definitorii între ele, obținem așa-numitul contrast definitoriu generalizant, care, ținând cont de relațiile (V ), caracterizează complet rezoluția replicilor de un grad ridicat de fracționare: * = Xj x x x Xj xg x = x x x (V ) Rezultă următorul sistem de evaluări comune: X, = Xj Xj X X = Xj x x = x x x , Xj \u d Xj X X \u d Xj X - Xj X X Xb, Xj = Xj X X = Xj X = Xj X X X , X = Xj Xj X = Xj Xj X Xj = X X , X = Xj Xj X = Xj Xj X Xb = X X , X = Xj Xj X X X = Xj Xj = x X , - * Zo + Zіzz + Zi "+ Zz -*■ + + Zai + "• bt -+■ " + + is + "> b + + " + ₽ B> bi -► Z + + "+ " > bt -► "+ + + Z > Xj X = X X = Xi " Xs = Xj X Xs> j -► j + ₽ + ? b + ?!"• XjX = X X = X X X = XjXjiXb, b -> ? +? S + ? b + ? l> (V ) Planul experimental corespunzător este prezentat în tabel Tabel Sferturi de replica de la PFE cu rapoarte generatoare (V ) Numărul experienței Xt X Xs X X Numărul experienței Xt X Xs x XB + + + + + + + - - + - - + + + - - - - + - + + - - - + - + - + - + - - + - - + - Rezoluția acestui sfert-replică este scăzută - toate efectele liniare sunt amestecate cu efectele interacțiunii perechilor DFE poate fi suplimentat la un experiment factorial complet prin implementarea replicilor fracționale lipsă În exemplul luat în considerare, pentru restul de trei sferturi de replici, relațiile generatoare vor fi: X =X X X , Xb= - X [X , x = - x x x , xx = Xi Xi, (V ) X =-XiX X > Xb=-XjX În acest caz, contrastele definitorii generalizatoare au forma / \u d Xj x x x - - Xj x x \u d - x x x > I = - Xi Xi x x = Xj Xi x = - x x Xi, (V ) I \u d - Xj Xi X X \u d - Xi Xi Xi \u d X X Xi Ca urmare a implementării acestor replici trimestriale complementare, se obțin estimări nemixte pentru toți coeficienții teoretici Numărul de experimente dintr-o replică fracțională trebuie să satisfacă inegalitatea k + /N, obținem /= Valoarea poate fi luată ca măsură a informațiilor conținute unde p este pătratul razei sferei în spațiul ^-dimensional, reciproca lui s în ecuația de regresie Conform (V ), cantitatea de informație scade proporțional cu pătratul razei sferei p și este aceeași pentru toate punctele echidistante (Fig ) Planificarea cu această proprietate se numește planificare rotativă Să prezentăm în termeni generali schema analizelor de dispersie și regresie a unui experiment planificat, când fiecare experiment din matricea de planificare a fost repetat de m ori (Tabelul ) Orez Proprietatea de rotație a unui plan liniar Tabelul Matricea de planificare și rezultatele măsurătorilor Numărul experimentului XO Xi X xk Y Y + + - + J'n J'IS, > Wit Jl ? + - - + U ■■ ■' U?t si + + + + N + - + - ■Ev ' -Ev ' ■ ■'yNm yN $ În fiecare linie a matricei de planificare, valoarea medie a valorii măsurate este determinată pentru m experimente paralele: T Ui~-= — - , i= , , , N- (V ) m si varianta T ^{Woo-Woo s? = ; i = , N (V ) 'm - Omogenitatea variațiilor eșantionului este verificată prin testul Cochran Pentru aceasta, raportul dintre variația maximă și suma tuturor variațiilor este compilat: c* °max Raportul rezultat este comparat cu cel tabelar: G, - ^ unde / = , ; / \u d m - l;f \u d N Dacă dispersiile sunt omogene Apoi, ca estimare a varianței de reproductibilitate, putem lua varianța medie cu numărul de grade de libertate / aOspR = N (m - / Coeficienții ecuației de regresie sunt determinați prin formula Ș HL Ui b}z=^N - ' (V- ) Având în vedere că varianța volumului t obținut din probă este de t ori mai mică decât varianța unei singure măsurări i s~ ~ seocnp !m' (V ) în exemplul luat în considerare (Tabelul ), varianța coeficienților zr se determină după cum urmează: ^=" play/'^- (V ) Semnificația coeficienților este verificată prin testul t Student Sub ipoteza nulă Н°:Ș; = , raportul dintre valoarea absolută a coeficientului ecuației de regresie și eroarea acesteia are o distribuție Student Pentru toți coeficienții ecuației de regresie, este compilat un raport z SUNT , sbj (V ) care se compară cu tabelul (f) pentru nivelul de semnificație p = , și numărul de grade de libertate f=N(ml) Daca atunci se acceptă ipoteza egalității la zero a coeficientului general de regresie = , iar coeficientul eșantionului corespunzător bj este eliminat din ecuația de regresie ca fiind nesemnificativ În același timp, datorită ortogonalității matricei de planificare, coeficienții rămași nu trebuie să fie recalculați Adecvarea ecuației de regresie la experiment se verifică în același mod ca și la prelucrarea unui experiment pasiv, după criteriul Fisher În matricea de planificare (Tabelul ), fiecare experiment a fost repetat de m ori Pentru a verifica caracterul adecvat, este compilat un raport de dispersie g aad / play > unde $ D este varianța de adecvare, determinată de formulă / este numărul de coeficienți semnificativi din ecuația de regresie Dacă raportul de dispersie rezultat este mai mic decât cel tabelar: E df-> df — gradf = -i + j + + -i-k, (V ) dxi dxg dxk —> —► —> unde i, j, k sunt ortele axelor de coordonate Se presupune că funcția f este continuă, diferențiabilă, cu o singură valoare și nu are puncte singulare Box și Wilson au propus o metodă în trepte pentru deplasarea pe suprafața de răspuns În vecinătatea punctului L, se realizează un experiment pentru o descriere locală a suprafeței de răspuns printr-o ecuație de regresie liniară: L y \u d bd + b x + b xi + + bxk (V ) Apoi, ne deplasăm de-a lungul suprafeței de răspuns în direcția gradientului aproximării liniare: DU h du ~ \u d Oi\u e a dx dx LA d-± = bh dhk (V ) La stabilirea experimentelor, dimensiunea pasului ar trebui să fie proporțională cu produsul coeficientului bj cu intervalul de variație bj&Zj Dacă o aproximare liniară nu este suficientă, atunci pune O nouă serie de experimente centrată pe punctul care corespunde celei mai mari valori a lui y și se găsește o nouă direcție pentru deplasarea de-a lungul suprafeței de răspuns Acest proces treptat continuă până la atingerea unei regiuni apropiate de extremum, sau „regiune aproape staționară” Direcția gradientului depinde de intervalul ales de variație a factorilor independenți Când intervalul de variație pentru un factor /-lea se modifică de n ori, dimensiunea pasului pentru acest factor se modifică de n ori, deoarece coeficientul de regresie bj se modifică de n ori și intervalul de variație se modifică și el de n ori Doar semnele componentelor gradientului rămân invariante la modificarea intervalului Alegerea cu succes a intervalului de variație este în mare măsură asociată cu disponibilitatea informațiilor a priori despre sensibilitatea parametrică a procesului Intervalul de variație ar trebui să fie suficient de mare, astfel încât intervalul de modificare a valorii de ieșire să fie de câteva ori (de cel puțin - ori) mai mare decât eroarea de reproductibilitate În același timp, pentru majoritatea proceselor, aproximarea liniară a suprafeței de răspuns este adecvată experimentului doar pentru intervale mici de variație Dacă nu se impun restricții asupra intervalelor de variație, acestea tind să fie alese astfel încât să se obțină o ecuație de regresie simetrică față de coeficienții termenilor liniari Prelucrarea rezultatelor unui experiment asociat cu o ascensiune abruptă trebuie să fie însoțită de o analiză statistică aprofundată a rezultatelor obținute Exemplul Compoziţia optimă a sticlei fotocromatice a fost determinată în sistemul LizO - AI O - SiOa Densitatea optică în starea iradiată a fost considerată ca parametri de optimizare (y) A fost necesar să se determine compoziția sticlei și condițiile de topire a acesteia, oferind densitatea optică maximă Au fost aleși ca factori independenți următorii factori: zi este concentrația inițială de clor, atom g/ g sticlă; zi este concentrația inițială de brom, g atom/ g sticlă; rz, co-purtător de Ag; CI; Z - temperatura de gătire, °C; zs este timpul de expunere, h; ze este conținutul de AI O , mol acțiuni; z? - raportul LizO / SіOg Condițiile experimentale sunt date în tabel G Z Z Z Z Z Z Nivel de bază zf> Interval variabil , , , , , , iaDg; , , , , , , + , , , , , - , , , , , , Soluţie Pentru a determina coeficienții ecuației de regresie liniară Y - + b x b x + b x + b % + bjXs + beXf + b x folosit '/te din PFE ' cu rapoarte generatoare X =X X X , x =x x > Xv=x x , Xi—X X Fiecare experiență din matricea de planificare (Tabelul ) a fost repetată de două ori Valorile medii ale densității optice y au fost obținute din două măsurători Să verificăm omogenitatea varianțelor s?, i = , , , , prin testul Cochran Suma varianțelor criteriul lui Cochran s? = > °- ~ Z Y Numărul experienței z> Z Y ZJ , , Lg- , , bj , , , , — bj tzzj , , , , — Pas , , , , — , , , , , , , , , , , , , , , Valorile luate ca pași sunt de , ori mai mari decât produsele bj Azj Cel mai bun rezultat a fost obținut în al -lea experiment O creștere suplimentară a concentrației de clor și a raportului Ag: CI înrăutățește proprietățile fotocromice ale sticlei În acest sens, au fost implementate experimentele ratate și Au fost obținute următoarele valori ale densității optice a sticlei: y = , , y = , Astfel, compozitia sticlei obtinuta in al -lea experiment este recomandata ca fiind cea optima Descrierea zonei apropiate de extremum Compozit Planuri Box-Wilson Regiunea apropiată de extremum se numește de asemenea o regiune aproape staționară Aceasta este o zonă cu semnificative neliniaritatea funcției de răspuns, pentru o descriere adecvată a cărei trebuie folosite polinoame neliniare În prezent cel mai larg pentru a descrie regiunea aproape de extremum, se folosesc polinoame de ordinul doi Acest lucru se datorează faptului că, în primul rând, există planuri de ordinul doi bine dezvoltate și, în al doilea rând, pentru că suprafețele de ordinul doi sunt ușor de sistematizat și studiat pentru un extremum Și în sfârșit, o creștere a ordinii polinomului de aproximare duce la o creștere semnificativă a numărului de experimente De obicei, un experiment implementat pentru a determina condițiile optime de proces poate fi descris în mod adecvat printr-un polinom de ordinul doi În acest caz, numărul de experimente N din plan nu trebuie să fie mai mic decât numărul de coeficienți determinați în ecuația de regresie de ordinul doi pentru k factori: L Y \u d * + + b X + + bk xk -f- bis Xr + + + bk— , k xb-i xh + *n + • • • + bkk x% (V ) Coeficienții ecuației de regresie (V ) servesc ca estimări pentru coeficienții corespunzători ai ecuației de regresie teoretică: tu — ₽ + ₽ X + • • • + ?*Xh + ? X X + + kXfc + + Pux?+ - + hfcxg (V ) (V ) (V ) Numărul de coeficienți / într-un polinom de ordinul doi poate fi determinat prin formula / \u d A + i + A + C \u d L + + - \u d (* + ) (* + ) k ! (A - ) unde C* este numărul de combinații de k factori cu doi, egal cu numărul de efecte de interacțiune perechi din ecuația (V ) Dacă regiunea aproape staționară poate fi descrisă în mod adecvat printr-o ecuație teoretică de regresie de ordinul doi (V ), atunci efectele de interacțiune a factorilor determinate experimental și efectele pătratice devin semnificative Acest lucru face posibilă stabilirea faptului de a fi într-o regiune aproape staționară Proximitatea unei regiuni aproape staționare poate fi stabilită dacă, pe lângă proiectarea factorială, în centrul planului se efectuează experimente k sau k ~ p (xi = ; xr = ; ; xk = ) iar media y se calculează Media y este o estimare pentru termenul liber al ecuației de regresie teoretică (V ) în timp ce coeficientul bo, calculat în experimentul factorial conform formulei N Wu Lee Wu * = ' N este o estimare comună pentru termenul liber și suma pătraticii tipic: h * ** ?«+ ₽// ,= (V ) Prin urmare, diferența k Wo-* t=i (V ) poate servi într-o oarecare măsură ca măsură a curburii unei suprafețe Pentru a descrie suprafața de răspuns prin polinoame de ordinul doi, factorii independenți din planuri trebuie să ia cel puțin trei valori diferite Un proiect pe trei niveluri, în care sunt implementate toate combinațiile posibile de k factori la trei niveluri, este un experiment factorial complet de * În tabel prezintă matricea de planificare pentru experimentul factorial complet Z Tabelul Numărul experienței x Y Numărul experienței x Y U, + L + U? - v? - Uz + - Uy + Y* - - Y, + + U Un experiment factorial complet * necesită prea multe experimente, depășind cu mult numărul de coeficienți / deja determinat pentru k > Numărul de experimente în PFE Zk și numărul de coeficienți / din ecuația de regresie de ordinul doi sunt prezentate mai jos: La * Puteți reduce numărul de experimente dacă folosiți așa-numitele planuri compoziționale sau secvențiale propuse de Box și Wilson Miezul unor astfel de planuri este PFE * cu k Posibilitatea de a folosi semireplicile ca nucleu al planului pentru k > se datorează faptului că deja semi-replica oferă obținerea de estimări nemixte pentru efectele liniare și efectele interacțiunii perechi (vezi Cap V, ) Dacă ecuația de regresie liniară s-a dovedit a fi inadecvată, trebuie să: ) adăugați k puncte stea situate pe axele de coordonate ale spațiului factorilor Coordonatele punctului stelei: Orez Plan de compunere de ordinul doi pentru k= (*a, O, , O), (O, ±a, O O) (O, O O, ±a), unde a este distanța de la centrul planului până la punctul stea - brațul stelei; ) crește numărul de experimente în centrul planului n Luați în considerare construcția planurilor compoziționale folosind exemplul k = (Fig ) Punctele ], , , din PFE , punctele , , , sunt puncte stea cu coordonatele (± a, ) și ( , ± a), coordonatele n experimente în centrul planului sunt zero - ( , ) (Tabelul ) Tabelul Plan de compoziție de ordinul doi pentru doi factori Numărul experimentului xo * X x,x X Numărul experimentului XQ * x x,x * * + + + + + + + -a a + + - - + + + + - - + + + + + - + - + + + +a a + -a a N + + +a a Matricea informațională (matricea momentelor) a compozitului de ordinul doi are pentru k = forma ba b b b bn brr ho O O NN X»iXi i=l = N b O X O br (xthx) = O O N Х І br O O O N bc y, I=G b Y XhX i i= oh oh XliX l N xu XUX i f=l î=! хІЛі X*, І = = (V !) Unde = NN S^ = ^ = + ^, (= = NN x"ix\i = x"ix i = + a , = = N ■ £s=tl = f ^ = * + a* /ssl fs=l Numărul total de experimente din matricea planului compozițional de ordinul doi cu k factori (Tabelul ) este # \u d * + A + n pentru A (V ) Tabelul Plan de compoziție de ordinul doi pentru k factori Numărul de experiență Xo k, x Xs Xk + + - + - + - - - + + + - + + - + - + + + - + - "I + - - - - n + + + a „i+ + -a „i + L + - a N + Matricea informativă XmX corespunzătoare planului (Tabelul ) are forma (V ) Tabelul „o la „ la * * , , , , , , , , ■ * Semi-replică, x = x^x^ Unde v * + а pentru k ■ , miez Л ; NN * ІХіі = X iX f = Y r + “ pentru k ; lg * PENTRU k ; — V g + “* PENTRU k = Astfel, nalny: planuri compoziționale de ordinul doi non-ortogo- k; Oh oh XjiXui=^ și' І= ' ' ' u=hj- Alegerea valorii brațului de stea a și a numărului de experimente în centrul planului n este legată de criteriul optimității planului Planuri ortogonale de ordinul doi Planurile de compoziție se reduc cu ușurință la o alegere ortogonală a brațului stelei corespunzător a Pentru a face acest lucru, matricea (V ) a fost inversată în formă generală În acest caz, a fost suficient să inversați acea parte a acesteia care este asociată cu coloanele x și x? (Tabelul ), adică cu coeficienții b și , și determinați a din condiția ca elementul off-diagonal al matricei inverse să fie egal cu zero: a | *a - *-i (k + O no) = , miez *; a | *- a *- (k + , n) = , miez *- (V ) Valorile lui a , determinate de (V ), sunt date în tabel (vezi p ) Alegerea unui din tabel și efectuând următoarea transformare liniară a coloanelor pătratice x/ = x/~ x/ = *? (V ) lg ѣ = NN %chi = lg obținem o matrice ortogonală Astfel, proiectul ortogonal de ordinul doi pentru k = și n = are forma (Tabelul ): Tabelul Numărul experienței x A, x X,' x Numărul experienței Ho * x X, X X,' x + + + + + / + / + - + / - / + + - - + / + / + + - / + / + - - + + / + / + - - / + / + - + - + / + / + - / - / + + + / - / Datorită ortogonalității matricei de planificare, toți coeficienții de regresie sunt determinați independent unul de celălalt prin formula N hpui (V ) ^ = iar varianţele coeficienţilor sunt L L = * Joaca bj N i=i Ca rezultat al calculelor pe o matrice cu coloane transformate pentru efecte pătratice, obținem o ecuație de forma L Y = b + &i*i + b x + • • • + bhxh + b xtx H-+ + b(k-\) k *fc-i xk + bn *î) + ••• +*fth ( ~xk)- (V ) Pentru a merge la notația obișnuită, determinați b prin formula bo ~ bo bch? • • • bfyț, x (V ) și evaluați cu o varianță egală cu = sb' + (*i) sbt, +•••+( Y sl ■ (V ) oh hh Cunoscând dispersia reproductibilității, verificați semnificația coeficienților și adecvarea ecuației L ft) unde fi -N- este numărul de grade de libertate ale dispersiei de adecvare; f este numărul de grade de libertate ale dispersiei de reproductibilitate; N este numărul de experimente din matricea de planificare; / este numărul de coeficienți din ecuația de regresie de ordinul doi Semnificaţia coeficienţilor se verifică prin criteriul lui Student: O = Coeficientul este semnificativ dacă tj > tp (fi), unde fi este numărul de grade de libertate ale dispersiei de reproductibilitate Coeficienții ecuației de regresie, obținuți cu ajutorul planurilor ortogonale de ordinul doi, sunt determinați cu o precizie diferită Conform (V ) avem; sbn ~ ^play/V^FI, sbj~ sBocnp / ^ * + a , j = , , , k pentru k , nucleu *" , ^p Cup pentru k = SB Cnp HRC> , miez ^, (V ) sb = play At (*/) j= > , , k, ș sb = " Bocnp — — - pentru k> , miez * " / * ( - ) + (a - + (n + k - )py) Planurile ortogonale de ordinul doi nu au proprietatea de rotabilitate Exemplul S-a studiat procesul de descompunere a polihalitului cu acid azotic: K SO ■ MgSO ■ CaSO • Н О + HNO -► - KNOS + Mg(NO ) + H SO + CaSO Sarcina a fost de a clarifica posibilitatea fundamentală a celui mai complet transfer al componentelor utile ale polihalitului de potasiu și magneziu într-un extract de acid azotic și de a analiza sensibilitatea parametrică a procesului Soluţie S-au ales ca factori de care depinde procesul de descompunere a polihalitului de către acid azotic următorii factori: zi este temperatura procesului, °С; yr este durata interacțiunii reactivilor, min; rz este norma acidului azotic, % în greutate din stoichiometrie pentru descompunerea părții de potasiu-magneziu a polihalitului; za este concentrația de acid azotic, % în greutate Datorită faptului că scopul acestui studiu a fost acela de a analiza sensibilitatea parametrică a procesului, s-a ales ca proiect experimental un plan ortogonal de ordinul doi, care asigură că toate covarianțele dintre coeficienții din ecuația de regresie sunt egale cu zero Coordonatele centrului planului, intervalele de variație și nivelurile studiului sunt date în tabel Zjz + , , , + - , , , + - , , , , În experimente s-au folosit deșeuri de polihalit din producția de sulfat de potasiu din minereul polimineral al uzinei Kaluga Compoziția deșeurilor de polihalit, %; KSO , ; Mg , ; CaO , ; S , ; reziduu insolubil (n o ) , ; compoziție mineralogică, %; K S , ; MgS = - , ; /> = , ; M = - , ; =- , ; c \u d - , ; b = - , ; = - , ; /> = - , ; L(,= , ; sb= , ; s*= , i ui ii Pentru y : b = , ; = , ; b = , ; = , ; = , ; = - , ; = , ; = , ; J = - , ; = - , ; = , ; c \u d - , ; = - , ; = - , ; = - , sb = , ; sb= , ; sb= , / / Semnificația coeficienților ecuațiilor de regresie a fost evaluată prin criteriul /-Student în conformitate cu formula (IV ): PENTRU I/I! eu! = , ; o, o ( ) - , După ce le-a îndepărtat pe cele neînsemnate coeficienți pentru care raportul / este mai mic decât cel tabelar, ecuațiile de regresie au forma ^ \u d , + , x + , - x + , x - Z ІZx® - , x - - , ^- , x , y \u d , + , xi + , xr + , x - , x "- , x® Orez Efectul temperaturii asupra gradului de extracție a K O (/) și MgO ( ) în soluție Orez Influența duratei de interacțiune asupra gradului de extracție a K O ( ) și MgO ( ) în soluție y,% ț/,% Orez Influența normei acidului azotic asupra gradului de extracție a KrO ( ) și MgO ( ) în soluție Orez Influența concentrației de acid azotic asupra gradului de extracție KzO ( ) și MgO ( ) în soluție Ecuațiile rezultate sunt adecvate experimentului: Fy, - , ; Fy - , Valoarea tabelară a testului Fisher /''tabel - , pentru nivelul de semnificație p - , , fi - , fi - și fi - , fi- Analiza sensibilității parametrice a procesului conform ecuațiilor de regresie este prezentată în fig - Se fac calcule pentru centrul planului Gradul de extracție a KrO și MgO în soluție crește odată cu creșterea temperaturii, duratei și vitezei acidului azotic (Fig - ) Dependența gradului de extracție a MgO și KrO în soluție de concentrația acidului azotic este extremă (Fig ) Valoarea extremului (maximului) pentru gradul de extracție a KrO este de , % în aceste condiții (în centrul planului), iar MgO este de , % la o concentrație de acid azotic de , % Din datele prezentate rezultă că, în toate condițiile studiate, MgO este extras mai rapid din polihalit în soluție decât KrO Prin urmare, la stabilirea condițiilor optime pentru procesul de descompunere a polihalitului cu acid azotic, s-a ales ca indicator principal gradul de extracție a KrO Ca urmare a rezolvării problemei de optimizare prin metoda programării neliniare, s-a constatat că în intervalul studiat de modificări ale factorilor, cel mai înalt grad de extracție a KzO în soluție ( , %) se realizează în următoarele condiții: concentrație HNOa , %, norma HNO - % din stoichiometrie, durata interacțiunii - minute În aceste condiții, MgO intră aproape complet în soluție Planuri rotative de ordinul doi Box - Hunter Planurile ortogonale de ordinul doi nu au proprietatea rotativă ness Cantitatea de informații Orez Linii de informare egală pentru un plan ortogonal de ordinul doi pentru k - definit ca reciproca lui sj se dovedește a fi diferită pentru punctele echidistante Pe fig arată contururile de informații egale pentru k \u d și planul dat în gabl Suprafețele de informații egale pentru un număr mai mare de factori sunt foarte complexe Box și Hunter au sugerat ca planurile rotative de ordinul doi să fie considerate optime Planificarea personalului lui Roth va fi o astfel de planificare, în care matricea de covarianță {ХТХ) ~ este invariantă la rotația ortogonală a coordonatelor Condiția de rotabilitate pentru modelele de ordinul doi este îndeplinită dacă momentele matricea informaţională satisface relaţiile N ^i = ^ , /= , k, (V ) i=l NN = wL = jVX"' "-/=r k, u^j, (v ) О= О= unde Â și Â sunt constante arbitrare Momentele rămase ale matricei informaționale sunt egale cu zero (vezi p ) În funcție de valorile lui k și X , profilurile de informații se modifică Pe fig Figura prezintă profiluri de informații pentru k= , X - la diferite valori ale lui X Profiluri similare sunt obținute pentru alte valori ale lui k O astfel de planificare se numește planificare rotativă uniformă Determinantul matricei (V ) include ca factor valoarea [( + k) X / X ] - k X ]-/ , i e X /) > A/(d + ) (V ) Să luăm în considerare în ce cazuri condiția (V ) nu este îndeplinită Să aflăm relația dintre constantele X și Z și numărul de factori k Înmulțind relația (V ) cu k, obținem N la NN k N * ^ = *> ^ = і=і /=і ,=і і=і /=і і=і De aici (V ) (V ) unde P, este raza punctului i din spațiul ^-dimensional Folosind relația (V J, obținem și ,V N k = CP • i= >= ;= (V ) La Să exprimăm suma în termeni de rază a punctului și numărul de factori k Ținând cont că pentru orice punct i-lea Fig Profiluri de informații pentru un plan rotativ pentru k - s > II * = zUfe(fe-l)^ (V ) obținem k x* = p'-A(AJ)z>J (V ) = Înlocuind (V ) în (V ), obținem NN N * xn = ₽i—* (*—*) xui zhl • "=ii=i "=i Ținând cont de relația (V ), avem N NN d s* xl~ ₽?-*(*-DJS "i-L= i=li=l De aici N n ₽î V g \u d - -! - \u d mp (V ) (V ) (V ) N ₽î X = t=l • Nk (k - ) Apoi (V ) N b ă '■ s h ( j '+ ' (V ) Dacă toate cele N puncte ale planului rotativ sunt situate pe aceeași sferă, i e Oh oh ₽ = P ^ ••• = ₽M ȘI p = p = = Rdg , din (V ) avem X'/ X| = k!(k - ) iar determinantul matricei (X X) este zero Prin urmare, este necesar ca punctele planului să fie situate pe mai multe sfere Pentru un anumit număr de factori, raza sferei pe care se află punctele miezului planului coincide cu raza a sferei pe care se află punctele stelelor Pentru ca matricea informațională să fie nedegenerată, în planul rotativ sunt introduse puncte care se află pe o sferă cu rază zero - n puncte în centrul planului Fie că N puncte ale planului sunt situate pe s sfere cu pi puncte pe fiecare sferă Ns Ns ₽:= ""₽: ^ = ^- (V ) I= U= i=l U= Pentru Ă = avem n"P" kH k+ ' (V ) Xj- Dacă punctele nucleului planului u și punctele stea na sunt situate pe aceeași sferă Рa = Рц și există n puncte în centrul planului, atunci pl + p, + Po \u d p - "o iar din (V ) avem x kN = (W ) * + (* + )nk + ' * Astfel, prezența a n puncte în centrul planului asigură îndeplinirea condiției (V ) Valoarea brațului stelei în planuri rotative poate fi determinată din relația (V ): la k , , k i=l Coeficienții bjj sunt corelați între ei și cu termenul liber bo Prin urmare, pentru a determina coeficienții ecuației de regresie, este necesar să se rezolve sistemul de ecuații normale prin inversarea matei ritsu(xg x) В = (хт X) Хт Y • Natura specifică a matricei (Xm X) pentru planurile rotative ne permite să efectuăm procedura de inversare a acestei matrice într-o formă generală și să obținem formule pentru calcularea coeficienților ecuației de regresie și a variațiilor acestora: ba \u d ~ X * (d + ) "d- X C l i=i i=i 'tsui unde k este un număr с N х>ххПу u=^j-, u, j = f== nu - X C •• LH (A + ) s N joc • De la = - s : y N redare * = ; x^ joacă , , хііуі — (V ) L/ factori, uh) [(fe+l)X -(fe-l)]C* play ' NN LA A l L = Constant X {(th + )X th] se determină prin formula (V ) Expresiile (V ) pot fi simplificate prin combinarea constantelor: k N b = a Ui - Oi i= /=i i=i bJ = a xlui "=M' I, u= l, , k; = N k NN bjj= аъ ~a# x ai i=l /= i=li=\ sX,= aA' sbj = a^zy> ib,J = a*s': ^ = (° + "c)^- (V ) Valorile constantelor incluse în expresiile (V ) sunt date în fila Tabelul Numărul de variabile independente k Numărul de experimente Po a a a ^ a ^ i = , ; /> = - , ; b = , ; b = - , ; *s = - , ; />n = - , ; bs = , ; *zz = - , ; bu = - , ; /> = - , ; />k = , ; />i = , ; /> = , ; />i = , ; bî = , ; * = - , ; />a = , ; = , ; /> = , ; * = , ; = , ; s c/ = , ; sbjj = , Semnificația coeficienților este verificată prin criteriul Student: ; ®^d = , ; F= , ; Г > ( , ) = , Ecuația rezultată face posibilă determinarea gradului de descompunere a concentratului de flotație de fosforit Karatau la diferite temperaturi în funcție de conținutul de impurități din acid Criterii pentru optimitatea planurilor În determinarea criteriilor pentru optimitatea planurilor pentru Box și școala sa, abordarea empirico-intuitivă este caracteristică La început, ei au sugerat ca planurile ortogonale să fie considerate optime, mai târziu - rotative Un plan este ortogonal dacă are o matrice de informații diagonală corespunzătoare Estimările parametrilor obținute din planuri ortogonale sunt independente Un design este rotativ dacă matricea de covarianță corespunzătoare este invariantă la rotația coordonatelor ortogonale Îndeplinirea acestei condiții face ca orice direcție din centrul experimentului să fie echivalentă în ceea ce privește acuratețea estimării suprafeței de răspuns Proprietățile ortogonalității și rotativității planurilor sunt extrem de convenabile din punct de vedere practic, ceea ce contribuie la utilizarea pe scară largă a acestor planuri în experimente mente Modelele liniare ortogonale și l^-p au, de asemenea, proprietatea de rotație Planurile compoziționale rotative propuse de Box și Hunter nu sunt ortogonale Dacă, totuși, ortogonalitatea este aleasă ca criteriu de optimitate, atunci unele pierderi în acuratețea estimărilor parametrilor și a funcției de regresie sunt inevitabile Concomitent cu dezvoltarea ideilor lui Box, s-a dezvoltat o a doua direcție, pur teoretică, în planificarea experimentelor Cea mai mare contribuție la dezvoltarea sa a avut-o matematicianul american Kiefer Conceptul de D-optimalitate dezvoltat de Kiefer este o extensie naturală a teoriei lui Fisher a estimatorilor eficienți În teoria lui Fisher, eficiența estimărilor este dată doar de modul optim de procesare a rezultatelor experimentului La procesarea experimentelor prin metoda celor mai mici pătrate pentru o ecuație de regresie liniară, se găsesc estimări comune efective ale acestor coeficienți În acest caz, elipsoidul de împrăștiere al estimărilor are cel mai mic volum Volumul elipsoidului de împrăștiere este legat de determinant în mod informațional; matrice după cum urmează Vh= — G unde G funcția gamma (V ) În conceptul lui Kiefer, eficiența este determinată și de aranjarea optimă a punctelor în spațiul factorilor Proiectarea experimentului, în care volumul elipsoidului de împrăștiere este minimizat pe setul de planuri dintr-o zonă dată, se numește D-optimal Conform (V ), proiectul D-optim trebuie să corespundă determinantului maxim al matricei informaționale Kiefer a propus o serie de criterii pentru optimitatea planurilor Toate aceste criterii, precum și criteriul D-optimalității, sunt de fapt reduse la unele cerințe pentru tipul de covarianță și, în consecință, matricea informațională Astfel, un design se numește A-optimal dacă matricea sa de covarianță are cea mai mică urmă (suma elementelor diagonale) Proiectarea L-optimală permite reducerea la minimum a variației medii a estimărilor parametrilor Un proiect se numește E-optimal dacă valoarea caracteristică maximă a matricei de covarianță corespunzătoare a estimărilor parametrilor este minimă Aceasta înseamnă că proiectarea L-optimă minimizează axa maximă a elipsoidului de împrăștiere a estimărilor parametrilor Un plan se numește G-optim dacă oferă cea mai mică variație maximă a valorilor prezise ale lui y în domeniul de planificare dintre toate planurile și, prin urmare, asigură că nu există puncte în domeniul de planificare la care precizia estimarea suprafeței de răspuns este prea mică Box și Draper propun un alt criteriu pentru optimitatea proiectelor, care permite reducerea la minimum a părtinirii sistematice și totale care apare atunci când suprafața de răspuns este aproximată printr-un polinom de ordin inferior celui cerut pentru o descriere adecvată În prezent, teoria construirii planurilor D-optimale și (/-optimale) este cea mai dezvoltată În general, problema construirii planurilor D-optimale nu a fost rezolvată Metodele de obținere a planurilor D-optimale pentru estimarea unui parametru pot fi considerate cele mai dezvoltate În lucrările lui Kiefer, Wolfowitz, Hola și Kono au introdus conceptul de proiectare continuă și au construit modele D-optimale continue pentru regresia polinomială de ordinul întâi și al doilea sub constrângeri asupra hipercubului și ^-dimensional minge, pentru regresia trigonometrică cu diferite funcții de greutate pe interval Marinele experimentului sunt date cel mai adesea de hipercub D -Proiectele optime de ordinul întâi sub constrângeri de cub pot fi date sub forma PFE * D-Proiecte optime sunt, de asemenea, niște replici fracționale ale experimentului factorial complet, iar proiectele Nluckett-Berman pentru numărul de factori k care satisface k + sunt un multiplu de patru planurile sunt în același timp ortogonale și rotative Proiectele continue D-optimale de ordinul doi pe cuburi de dimensiuni - pentru polinoamele de ordinul doi, propuse de Kiefer și Wolfowitz, de regulă, conțin un număr foarte mare de observații; deci, de exemplu, la un astfel de plan exact ar trebui să existe mai mult de de măsurători În acest sens, cu ajutorul unui computer digital, s-au găsit planuri asimetrice de ordinul doi cu un număr destul de mic de puncte experimentale, care se apropie de D-optimal în caracteristici precum determinantul matricei informaționale, media și maximul varianta a respectiva valoare a parametrului de optimizare S-a efectuat și o evaluare comparativă din poziția de D-optimalitate a caracteristicilor unor planuri compoziționale de ordinul doi sub restricții asupra cubului pentru k - , , Alegerea unuia sau altuia plan de studiu este determinată de formularea problemei şi posibilităţile experimentului Investigarea suprafeţei de răspuns Rezolvarea problemei de optimizare O ecuație de regresie de ordinul doi care descrie în mod adecvat regiunea aproape staționară este examinată pentru a determina coordonatele optimului În plus, este de interes să se studieze proprietățile suprafeței de răspuns în vecinătatea optimului În acest caz, este util să treceți de la polinomul de ordinul doi obținut din rezultatele experimentului la ecuația standard, canonică: Y - ys - + X + + \kk X- k > (V ) unde este valoarea de ieșire în centrul suprafeței; Хі, Хі, Хк — variabile canonice care sunt funcții liniare ale factorilor x,, x , , xk, Ă , Ă , , kkk — coeficienți ai formei canonice Prima etapă a transformării canonice este transferul originii coordonatelor la un punct singular al suprafeței de răspuns - centrul suprafeței Coordonatele centrului se determină prin rezolvarea sistemului de ecuații Orez Suprafețele canonice și secțiunile lor pentru k - Când aproximăm suprafața de răspuns cu un polinom de ordinul doi, trebuie să rezolvi un sistem de k ecuații liniare Dacă determinantul acestui sistem este egal cu zero, atunci suprafața nu are centru În acest caz, puteți fie să mutați originea în punctul cu cea mai bună valoare de ieșire, fie să nu mutați deloc centrul În acest caz, pentru o suprafață necentrală, optimul se va afla la limita domeniului de definire a factorilor Dacă suprafața are un centru, atunci originea este transferată În acest caz, termenii care conțin efecte liniare dispar în ecuația de suprafață, iar termenul liber se modifică Coeficienții la puterile secunde și interacțiunile sunt invarianți la transfer A doua etapă este rotirea axelor de coordonate în noul centru în așa fel încât termenii cu efecte de interacțiune să dispară; termenul liber este invariant sub rotaţie Ca rezultat, obținem o ecuație de forma (V ) Suprafețele de ordinul doi sunt clasificate în funcție de formele lor canonice (Fig ) Toți coeficienții formei canonice au aceleași semne Suprafața este un paraboloid eliptic (Fig a) În centrul suprafeței, maximul este la Ă ,, Coeficienții au semne diferite Suprafața este un paraboloid hiperbolic, o „șa” (Fig , b) În centrul suprafeței - "minimax" Unul sau mai mulți (dar nu toți) coeficienții sunt aproape de zero În acest caz, centrul se află mult dincolo de zona de experimentare Acest tip de suprafață se numește „altitudine în creștere” („cresta”), Este posibil și un caz degenerat de planuri paralele, care nu prezintă interes practic (Fig , c) Când A \u d (Fig , d), deplasând originea în punctul (de obicei aproape de centrul planului), obținem ecuația unei parabole; y - Uv - Ac-^i + A AC- (V ) Să trecem din ecuația de regresie de ordinul doi pentru k= obținută din datele experimentale Y = b + b x + b x -|- b xtx C- btlx + b x^ , (V ) la ecuaţia canonică (V ) Să determinăm coordonatele punctului , centrul suprafeței Pentru a face acest lucru, este necesar să rezolvați un sistem de ecuații; du l = , + ft x + ft, = ; (V ) face - \u d , x + Іb x + b \u d dx Rezolvarea sistemului (V ) dă coordonatele centrului xls și x Substituindu-le în ecuația de regresie (V ), obținem valoarea mărimii de ieșire în punctul S-ys Mutăm originea coordonatelor în punctul (y , x , x ) Vechile coordonate xp x , y sunt legate de noile rapoarte x', x , y'; * = * + *i • x = * + * ' (V ) Y \u d SW + y' ■ În noul sistem de coordonate, ecuația (V ) ia forma y - ys = U (*]') + b ( x' ) +* X] x' În a doua etapă a transformării, vom scăpa de efectul de interacțiune prin rotirea axelor de coordonate Pentru a face acest lucru, este necesar să rotiți axele de coordonate cu un astfel de unghi încât ctg a = (bn - b ) / b (V ) Apoi obținem în noul sistem de coordonate X , X , y ecuația de regresie în forma canonică: y-^ = XuXf + X X (V ) Vechile coordonate xp x sunt legate de noile relații * \u d (Xi + X ) cos a - (X - x s) sin a, (V ) X - (Xj + Xls) sin a - (X + X s) COS a Pentru a determina coeficienții ecuației canonice Хп și Ă , folosim următorii doi invarianți ai ecuației (funcții ale coeficienților care nu își modifică valoarea sub nicio transformare de coordonate): A = = bn + b - const, (V ) * / * '/ * * = const Din (V ) avem bc + b - xn + x > *n b = Xn X - /Xf (V ) Deoarece în această transformare X( = ), obținem rapoartele pentru determinarea Xc și X : * + * = Xc + X • (V ) *u * - !^b\ \u d Chi X - Folosind teorema Vieta, Chi și X pot fi determinate ca rădăcini ale ecuației pătratice: K" - (ftu + b ) X + (ftn b - '/" *? ) == (V ) În problemele multidimensionale, transformarea canonică se realizează folosind metode de algebră liniară Compuneți din coeficienții ecuației de regresie de ordinul doi obținuți prin experiment, Y \u d fi + bixi + • • • + bkxk + * X + • • • + bk-i, k xk-i xk + + ^n*î + ••••+■ bkk matrice pătrată B= bc i/rir bgg Ch * k - zbik (v ) Chgkі /gk bkk unde bj = bjj Pentru a determina coeficienții Іц, , Xkk ai ecuației de regresie în forma canonică (V ), este necesar să se găsească rădăcinile polinomului caracteristic Рк (X) al matricei В- /'d(X) = |fi X£|, (V ) unde E este matricea identității sau Al X Chgb - / sbik І b b X ^(X)= (V ) Vzbfci Chgbkg bhk-X O transformare liniară ortogonală pentru coordonatele X este dată de sistemul de ecuații: xi) + m k(xr - xm) - - m k (xk - xl), Xt = /n (x - xls) + m (xr - xOS) * ••- + m k (xk - x ^), (V ) Xh = mki (xi - xis) + mk (Xi - - + mkk (xk - xks) Coeficienții mik sunt soluții ale sistemelor de ecuații omogene Pentru X, sistemul va arăta ca ( ts - X,) Шіі + / i Mj + - ѵ blk mik \u d O, : : : (V ) / bklmH + /abk mI + • • • + (bkk - Xj)mih = Deoarece soluțiile ecuațiilor sunt doar proporționale cu acele valori care sunt necesare pentru ortogonalizarea transformării liniare (V ), acestea sunt recalculate, ținând cont de condiția de ortogonalitate; "a+m? +•••+= '"= • , , L (V ) Calculând coeficienții ecuației de regresie în formă canonică, determinând astfel tipul suprafeței de răspuns Tipul suprafeței de răspuns determină strategia de căutare extremum Dacă suprafața de răspuns este un paraboloid eliptic și Xn și Ă > are valoarea maximă la limita zonei de studiu Dacă suprafața de răspuns este un paraboloid hiperbolic L Uz - Al X[ - X x (V ) și sunt determinate condiții care oferă valoarea maximă a parametrului de optimizare, sunt stabilite de valorile y>y la Xz = și se deplasează de-a lungul axei canonice X], care are un coeficient canonic pozitiv Totodată, se verifică îndeplinirea constrângerii Xj =± a, înlocuind valorile (V- ) în formule (V ) În problemele multidimensionale, căutarea condițiilor optime de proces se realizează pe un computer, utilizând de obicei una dintre metodele de programare neliniară Exemplul S-a obtinut o ecuatie pentru regresia gradului de descompunere a concentratului de flotatie al fosforitului Karatau pe temperatura si impuritatile continute in acidul fosforic: Y \u d , + , el - , x - , + , - , + , xx x Se impune determinarea condițiilor corespunzătoare gradului maxim de descompunere (UmaxL sub restricțiile impuse de o sferă cu raza egală cu brațul stelar (Tabelul ) Soluţie Pentru a determina condițiile pentru gradul maxim de descompunere, variabilele, natura influenței cărora este clară din ecuația de regresie, sunt luate egale: xr = + ; X = - Efectul concentrației de SOa în acidul fosforic este reprezentat în ecuație prin termeni liniari pozitivi și patratici negativi Scara optimă a lui xa, egală cu , , este determinată din condiția valorii extreme a lui r în xs Cu aceste valori ale factorilor xs, xs și x , ecuația de regresie va lua forma y \u d , - Sx] + , xx x Aducem această ecuație la forma canonică Coordonatele centrului suprafeței sunt: \u d - , - x + , x \u d , xX - ; x - , Oz - , Astfel, centrul suprafeței coincide cu centrul planului Polinom caracteristic Pk (X) = — , —X' , , —X Ph (X) \u d X + , X - , \u d\u d sau % Orez Căutați un extremum în prezența restricțiilor y^const -La Orez Hiperbole de ieșire egală Rădăcinile polinomului Xi sunt + , , Xr sunt - , Ecuație în formă canonică y - , \u d , XI - , X Suprafața de răspuns este un paraboloid hiperbolic În secțiuni ale suprafeței de răspuns prin plane y = const sunt hiperbole (Fig ) În centrul suprafeței - minimax Transformarea liniară este dată de sistemul: Xg \u d , + , x , X \u d - , , x Pentru a determina gradul maxim de expansiune, ieșim din minimax de-a lungul axei X\ (coeficientul formei canonice este pozitiv), echivalând Xb cu zero: Creșterea lui Y verificăm îndeplinirea condițiilor хі = În condițiile optime obținute (XI- + , ; X -+ ; xa-+ , ; l-+ , ; xb - ) și (es, , ; x - + ; -+ , ; x - - , ; x - - ) au fost efectuate experimente de control Gradul de descompunere s-a dovedit a fi egal cu , , respectiv , % Astfel, discrepanțele cu valorile calculate se află în limitele erorii experimentale (vy = V = ) Dacă procesul este descris prin mai multe ecuații de regresie, este necesar să se rezolve o problemă de compromis - să se determine valoarea extremă a unei funcții de răspuns sub restricțiile impuse de alte funcții de răspuns și limitele zonei de studiu (Fig ) Să se găsească extremul funcției y=f(хі, ,хк), care depinde de k variabile х^(/-= , ,k), legate la rândul lor prin relații = ) corespunde unei valori absolut inacceptabile a acestui răspuns; și d = l (D = ) este cea mai bună valoare a răspunsului, iar îmbunătățirea lui ulterioară este fie imposibilă, fie lipsită de interes Valorile intermediare ale dezirabilității și notele lor numerice corespunzătoare sunt date în tabel Tabelul Semnele de bază ale scalei de dezirabilitate Notă cantitativă pe scara dezirabilității Dezirabilitatea valorii de răspuns , , , + , , + , , + , , + , Foarte bun Amenda Satisfăcător Prost Foarte rău Această alegere a mărcilor numerice se explică prin comoditatea calculelor, deoarece = , ®!-— și = , ®- (V ) a ei Construit conform tabelului , scala d este o scară fără dimensiuni prin care orice răspuns poate fi convertit astfel încât să poată fi interpretat în termeni de utilitate sau dezirabilitate pentru orice aplicație particulară Cel mai simplu tip de transformare este acela în care există limite superioare și/sau inferioare de specificație, iar aceste limite sunt singurele criterii de calitate care nu pot fi modificate În afara acestor limite, valoarea d = Q,O, între ele valoarea rf==l Funcția de dezirabilitate parțială cu restricție unilaterală (Fig , a) are forma I , Y Y !/min Orez Cea mai simplă formă a unei anumite funcții de dezirabilitate iar este foarte bun , - bine , - , - , - , - satisfăcător rău - oh, oh- foarte J rău / -z- - ' Y; Tip Utah U Orez Funcția de dezirabilitate pentru constrângere bilaterală În mod similar, se obține o funcție de dezirabilitate parțială dacă specificația specifică o constrângere superioară Dacă există o restricție cu două laturi pentru o proprietate dată (Fig , ), atunci (O, U Umax d=J I • Umin ^Y ymin, o altă dependență exponențială este o formă mai convenabilă a transformării yvb (Fig ): d = exp [-exp(-y')] Orez Funcția de dezirabilitate pentru o proprietate mărginită de (V ) pe o parte În exprimare (V ) y' = bt + biy (V ) Coeficienții bo și b\ pot fi determinați prin setarea valorilor corespunzătoare de dezirabilitate d pentru două valori de proprietate y, de preferință în intervalul , + і>і = , , (V ) - &і = - , Soluția sistemului (V ) dă ii - , , io - , Astfel, funcția de dezirabilitate parțială are forma di = exp[—exp(— , - , і/і)] (v ) Funcțiile de dezirabilitate parțială rf și Numărul experienței Ui ^ l , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , LS , , , , , , , , І , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Pe baza prelucrării datelor experimentale, eliminării coeficienților nesemnificativi prin criteriul Student, s-au obținut următoarele ecuații de regresie: Wu! - , - , xx - - , x + , x + , x - , D -J- b b ^ Xg - - , xt x + , x x + , * h xv, y \u d , , l- ■ , * + , * - , * - , x , - , * * + - - , * * - , * * + , * * in! - , * ; y \u d , + , * s + , - , x + , *n + , * x + , * * - - , * * + , * * -;- * * s + , x * + , x x - , x x - - , * s * in - , x * b + , * * in + , + , - , * |; y \u d , - , * - , * - , * - , *v - , * Xg - - , + - , XJ - , x? - , x? Verificarea prin criteriul Fisher a arătat că ecuațiile rezultate descriu în mod adecvat experimentul Ecuațiile pentru - y fac posibilă calcularea valorilor indicatorilor de proces selectați pentru orice combinație a factorilor studiați din intervalul studiat Rezultatele unor astfel de calcule sunt prezentate în Fig - Pe fig arată dependenţa vitezei de cristalizare a hemihidratului (yi) de concentraţia ionului SOa în diferite condiţii de proces Dependența lui de concentrația de SO are o extremă Orez Dependența vitezei de cristalizare a sulfatului de calciu hemihidrazei de concentrația de SO ^ în soluție în diferite condiții - Orez Influența individuală a factorilor asupra ratei de cristalizare a hemihidratului: - SiFs'; - AI O ( % Fe O ); - Р О : - temperatura; - ResOs ( , % AIsOs) Valorile perechilor rămase corespund centrului planului - valorile tuturor variabilelor corespund centrului planului: PrOk - , %; SiFg' „ , %; AI O - , %; RegOz - % temperatura ° - - , % SiF -; ~ °C; - fără Fe O? şi , % ArOe; o - % P O ; valorile variabilelor rămase pentru curbele , , corespund prețurilor și py din plan caracter În acest caz, poziția și magnitudinea extremului (maximului) depind de alte condiții pentru procesul de cristalizare a hemihidratului (temperatura, concentrația P O , impuritățile SiFe~, AI O și P'erO ) Absența în ecuație a efectelor interacțiunii dintre conținutul de SO? și temperatura, precum și conținutul de SO? iar oxizii asociați duce la faptul că o creștere a temperaturii și a conținutului de R O modifică doar valoarea absolută a maximului Prezența efectului de interacțiune între concentrația de P O și SO? provoacă nu numai o modificare a valorii extremului, ci și o deplasare a maximului în regiunea valorilor SuJ inferioare (curba Pe fig Figura prezintă curbele care caracterizează influența individuală a factorilor (xr - ce) asupra vitezei de cristalizare a hemihidratului cu valorile variabilelor independente rămase corespunzătoare centrului planului Din figură rezultă că odată cu creșterea concentrației de acid fosforic, viteza de cristalizare a hemihidratului crește, iar odată cu creșterea concentrației de acid fluorosilicic, aceasta scade Odată cu creșterea concentrației unuia dintre oxizi hemihidratați în absența celuilalt în soluție, viteza de cristalizare a hemihidratului scade Din ecuația pentru yi rezultă însă că, odată cu prezența simultană a ambilor oxizi în soluție, efectul negativ al fiecăruia dintre ei scade și, la o anumită concentrație a unuia dintre ei, dispare O creștere a temperaturii duce la o creșterea ritmului procesului Pe fig Figura prezintă grafice ale dependenței duratei de filtrare a pulpei formate (y , s) de condițiile procesului de cristalizare, rezultă din Fig că o creștere a concentrației de acid fosforic are cel mai mare efect negativ pe procesul de filtrare (curba ), deoarece aceasta crește nu numai vâscozitatea soluției, ci și dimensiunea cristalelor de hemihidrat care cristalizează din această soluție scade O creștere a conținutului de AI O , ResO și SiF« crește vâscozitatea Orez Influența temperaturii și compoziției soluției asupra duratei filtrării pulpei; - temperatura; - SiF^; - FeaOa; - Algos; - PrOzsvob Valorile parametrilor rămași corespund centrului planului soluție osoasă, care mărește durata de filtrare (curbele , , ) Dimpotrivă, o creștere a temperaturii (curba ) reduce vâscozitatea și, în consecință, reduce timpul de filtrare Pe fig prezintă grafice ale interacțiunii dintre concentrațiile de A O și SiFg ~ în soluție: la concentrații de A O până la %, cu o creștere a conținutului de ion SiFȘ-, viteza de filtrare crește; la concentraţii de A O de – %, cu o creştere a conţinutului de SiF , viteza de filtrare scade Astfel, efectul pozitiv al prezenței în comun a ionilor AI + și SiFj în soluție în condițiile producerii acidului fosforic prin metoda hemihidratului este limitat la un interval mic de concentrație Pe fig prezintă grafice ale dependenței captării ionilor de fosfat de către precipitatul hemihidratului (noduri,%) de compoziția soluției și a temperaturii Dependenţa conţinutului de P O din hemihidrat de ionii SO în prezenţa impurităţilor R O şi SiFe - are un caracter extrem Poziția și mărimea extremului (minimului) sunt determinate de temperatura și concentrația de impurități (Fig ) Pe fig arată influența individuală a impurităților SO, (curba /), SIFg (curba ) și A O (curba ) asupra cantității de captare a P O, prin cristalizarea hemihidratului dintr-o soluție de acid fosforic care conține % P O și , % A O la o temperatură de °C În absența impurităților Fe C) , SiF și a unui conținut scăzut de A O ( , %), o creștere a concentrației de SOJ- în soluție duce la o „scădere lină a pierderi P O (curba /) În aceste condiții, o creștere a conținutului de SiF ~ în limitele studiate reduce brusc pierderea de P O (curba ), care, după toate probabilitățile, este asociată cu formarea de silicofluorură de calciu, care este insolubilă în soluții concentrate de acid fosforic , și apariția în fază lichidă a cantității corespunzătoare Orez Influența conținutului de AIRO și SiFe” asupra duratei filtrării pulpei U - - - - - S £',% Orez Influența temperaturii și compoziției soluției asupra cantității de captare a P O : - centrul planului; - % Р О ; - % A O , - , % A O ; - °С Valorile factorilor rămași corespund centrului planului calitatea ionilor SO?" Prezența impurităților în soluție afectează disociarea acidului fosforic În această soluție, ionii fosfat, aparent, sunt sub formă de compuși complecși cu impurități Gradul de disociere a unor astfel de compuși complecși depinde Orez Influența individuală a factorilor asupra mărimii cristalelor de hemihidrat- I - Si - PerOz- - temperatura; - Provsvob - - SOp; valorile parametrilor rămași corespund centrului planului Orez Influența individuală a factorilor asupra cantității de captare a P O : ~ SOJ - SiFg-,- - AIaOz,- - P'Aevob, valorile parametrilor constanți corespund valorilor minime ( % P O liber; , % AI O ; fără SOJ-, SiF|', AIUz) la °C; - temperatura; - Al ; valorile factorilor rămași în construcția curbelor , corespund centrului planului nu numai la temperatură, ci și la raportul componentelor din soluție O astfel de presupunere face posibilă explicarea prezenței în ecuația pentru uz a efectelor de interacțiune între temperatură și conținutul de SO ', SiF|\AbO , P O și între conținutul componentelor din soluția de SO - AbO ; Od ~ - P O ; SiFg - AI O ; AI O - P O ; AI O - Fe O , P O - Fe O Pe fig prezintă grafice ale mărimii cristalului în funcție de compoziție și temperatură În gama de factori studiate, se formează cristale hemihidratate în formă de ac, astfel încât lungimea cristalelor caracterizează dimensiunea acestora Valoarea medie a lungimii cristalului (noduri, μm) a fost determinată pentru fiecare experiment prin media lungimii a - de cristale Din fig Din fig rezultă că într-o soluție a cărei compoziție corespunde centrului planului, cu o creștere a concentrației de acid fosforic, dimensiunea cristalelor scade (curba ) \\ după toate probabilitățile, acest lucru se datorează faptului că o creștere a vitezei de cristalizare a semihidrei (vezi Fig , curba ) Dependența mărimii cristalului de concentrația de SO| este extremă (curba ) În intervalul de la la ~ , % SiF,? o creștere a concentrației de acid fluorosilicic duce la o creștere a solubilității CaSO • , H O, o scădere a vitezei de cristalizare a hemihidratului (vezi Fig , curba ) și o creștere a dimensiunii cristalelor Cu o creștere suplimentară a concentrației de SiFg, este posibilă cristalizarea micilor cristale de silicofluorura de calciu pe suprafața cristalelor hemihidrat și inhibarea creșterii acestora În mod similar, o creștere a concentrației de ErO în soluție duce la o scădere a vitezei de cristalizare a hemihidratului (vezi Fig b, curba ), în urma căreia dimensiunea cristalelor de hemihidrat crește ușor (curba ) Odată cu o creștere suplimentară a concentrației de FeO , o parte din acesta precipită, ceea ce duce la inhibarea procesului de difuzie a ionilor de Ca + și SO la suprafața cristalelor hemihidrat în creștere Dependența mărimii cristalului de temperatură este extremă: o creștere a temperaturii duce la o creștere a vitezei de difuzie, ceea ce facilitează creșterea cristalelor, dar o creștere a ratei de cristalizare a hemihidratului (vezi Fig , curba ) implică: precipitarea cristalelor mici (curba ) Ecuațiile de regresie rezultate pentru y, - yt au fost folosite pentru a rezolva problema optimizării procesului de cristalizare a sulfatului de calciu hemihidrat O analiză a sensibilității parametrice a procesului a arătat (Fig - ) că natura influenței factorilor controlați (concentrația SOJ, P O și temperatură) este semnificativ diferită După cum sa menționat deja (p ), una dintre cele mai de succes modalități de a rezolva problema optimizării proceselor cu un număr mare de răspunsuri, lipsite de dificultăți de calcul, este utilizarea așa-numitei funcție de dezirabilitate generalizată O propusă de Harrington ca criteriul de optimizare generalizat Pentru a construi o funcție de dezirabilitate generalizată D, este necesar să se convertească valorile de răspuns măsurate într-o scară de dezirabilitate adimensională d Construcția scalei dezirabilității stabilește relația dintre valoarea răspunsului y și valoarea corespunzătoare d (funcția privată de dezirabilitate) În cazul nostru, există restricții unilaterale asupra parametrilor de ieșire de forma y yt;n O formă convenabilă de conversie a y în d este relația exponențială (V ): coeficienții io și ii pot fi determinați prin setarea a două valori ale proprietății y la valorile de dezirabilitate corespunzătoare d, de preferință în intervalul , , atunci fie va fi necesar un număr mare de experimente, fie va fi necesar să se limiteze numărul de factori introduși în plan cu valoarea ( + ,) În plus, pentru unii factori calitativi, uneori nu este posibil să se stabilească mai mult de două niveluri În astfel de probleme sunt utile planurile complexe: experiment factorial k, combinat cu un pătrat latin de dimensiunea kX k Vă permit să introduceți mai mulți factori în planificare la niveluri = k și un număr destul de mare de factori cantitativi și calitativi la două niveluri Astfel de planuri pot fi construite numai pentru un experiment factorial k cu numărul de experimente egal cu pătratul întreg al numărului k, k - , , Tabelul -»(-!) și (+ ) хі(- ) хі(+ ) зсі(- ) X (- ) xs(- ) A B C D xs(+ ) B A D C X (+ ) xs(- ) D C B A xs(+ ) C D A B Pentru a combina experimentul factorial k cu pătratul latin, este convenabil să prezentați experimentul factorial k ca un tabel cu intrări k+], pe care se suprapune un pătrat latin de dimensiunea kX k, de exemplu, Tabel Atunci factorul introdus în planificare după schema pătratului latin este ortogonal la factori care definesc experimentul factorial complet Într-adevăr, toate nivelurile I = k ale acestui factor apar la fel de des în proiectare și fiecare nivel al acestuia apare cu orice nivel al factorilor inițiali de k de același număr de ori Planul original poate fi combinat cu un pătrat greco-latin de kX , sau chiar cu un pătrat hiper-grecolatin obținut prin suprapunerea ( k- ) pătrate latine ortogonale dacă există un set complet de pătrate latine ortogonale pentru un dat = k În acest caz, factorii introduși ( ^- ) sunt ortogonali cu cei k factori originali și sunt, de asemenea, ortogonali cu toate interacțiunile factorilor care definesc coloanele pătratului Planul va fi saturat dacă aceste interacțiuni sunt considerate nesemnificative și utilizate pentru a introduce factori suplimentari în plan la două niveluri De interes sunt cele mai diverse variante de desene ortogonale saturate obținute prin combinarea designului factorial de k cu un pătrat latin, două pătrate latine ortogonale, etc până la ( * -!) pătrate latine ortogonale Fiecare factor introdus în proiect la niveluri = k are ( fc-l) grade de libertate și este amestecat cu L- interacțiuni diferite a k factori ai experimentului factorial complet Dacă introducem m factori în plan (m nJ = ctOU]/-/N ^(nb), efecte de interacțiune libere de amestecare cu factorii principali Efectele factorilor introduși în plan la niveluri / = * se calculează separat pentru fiecare nivel Efectul factorului xj(J = n-t + , ,n) la nivelul q-lea (?^= , , , , - ) este egal cu suma răspunsurilor în toate experimentele în care factorul xj este setat la nivelul q, împărțit la numărul de apariții (/ = *) în proiectarea factorului xy la nivelul q, (* ••• */ ••• HP (v ) Dacă există motive pentru a presupune omogenitatea variațiilor în măsurarea răspunsului în toate experimentele, atunci testul r poate fi aplicat pentru a evalua semnificația diferenței dintre efectele acestor factori la diferite niveluri Dezavantajul acestui criteriu este că atunci când se evaluează semnificația diferenței dintre efectele acestor factori, de exemplu, x , la două niveluri / și /+ , nu sunt folosite toate informațiile, ci doar o parte a acesteia Testul de rang multiplu al lui Duncan face posibilă determinarea semnificației diferenței dintre efectele nivelurilor de factori introduse în plan la /> niveluri cu o mai mare fiabilitate, întrucât toate informațiile obținute în experiment sunt utilizate simultan Semnificația efectelor principale ale factorilor introduși în plan, atât la două, cât și /> niveluri, poate fi verificată folosind analiza multivariată a varianței și analiza factorială Pe baza rezultatelor analizei factoriale se poate face o ascensiune abruptă Pentru model liniar Ui ІЯ f = + Х І + • • • + \n-t)i + X(P-t + ) q + + Е і jq r — y, A^ = , > , — diferență semnificativă y —y — , = , ; y( >= , ; y( >= , ; y ( ) — y O) = , > , — diferență semnificativă y — y ( ) , > , — diferență semnificativă y ( ) — y ( ) , — y = , - pentru I> - , " - efectul este semnificativ == j = , > , — efect semnificativ , z , - efectul este invizibil N = = , > , — efectul este semnificativ , — efect semnificativ ) Efectul factorului x s-a dovedit a fi nesemnificativ; astfel, excesul de solvent nu afectează randamentul felavinei Creșterea timpului de extracție, creșterea numărului de extracții și creșterea temperaturii, i e restul de patru factori conduc la o creștere a randamentului felavinei Semnificația diferenței dintre efectele factorului xe introdus în proiect la patru niveluri a fost testată folosind testul de rang multiplu al lui Duncan Tabelul arată valorile producției medii pentru diferite niveluri ale factorului xe: Tip de solvent Niveluri de factori Randament mediu, % Metanol , % etanol , % alcool izopropilic % metanol , * Clasăm valorile medii de ieșire în ordine crescătoare; y( ) yC>) jXO / ) Eroarea de reproductibilitate a rezultatelor sy este de , % Eroarea mediei este " \u d V% p \u d K ' \u d - %" Scriem rangurile semnificative din tabelul Duncan pentru nivelul de semnificație p - , și numărul de grade de libertate Rp~/VOsp* R Ranguri , , , Cele mai puțin semnificative ranguri (LSR) înmulțite cu eroarea mediei sy sunt: p WIP X Sy , , , Să verificăm semnificația diferenței dintre mijloace; y'( > - = , - , = , > , - diferență semnificativă -y , - diferență semnificativă y ( > -= , - , = , - y ( ) = , - , = , > , - diferență semnificativă — , = , diferență semnificativă > y — y, > = , — , = , ' i=l Pentru un plan simplex, conform (V ), (V ) în timp ce pentru planuri precum k și k~r JOACĂ *bg n' Astfel, coeficienții ecuației de regresie obținuți prin planul simplex sunt determinați cu mai puțină acuratețe Este posibil să se construiască modele saturate cu elemente ± numai pentru numărul de factori egal cu a - , unde a este un număr întreg pozitiv De exemplu, pentru , , , etc factori Pentru utilizarea practică a matricei simplex (V ), valorile numerice ale elementelor sale sunt precalculate folosind formula (V ): , , , , , , - , , , , , , - , , , , , - , , , , x \u d - , , , (V ) - , , — , Proiectarea experimentului la scară adimensională pentru k factori constă din k coloane și k + rânduri ale matricei (V ) După implementarea simplexului original, rezultatele sunt analizate pentru a identifica cel mai rău punct Apoi, punctul cel mai rău este reflectat în raport cu centrul feței opuse a simplexului și astfel se găsesc condițiile pentru efectuarea unui nou experiment în locul celui exclus Condițiile pentru efectuarea experimentului în punctul reflectat pot fi găsite după cum urmează: x^+ ) = x este coordonata y a noului punct rezultat din reflexie; x^-J-a coordonată a centrului feței opuse: &-H 'L k (V ) Simplexul fc-dimensional original poate fi extins la unul (k + ])-dimensional prin introducerea unui singur punct nou O astfel de nevoie apare dacă, în prima etapă a studiului, dependența procesului studiat a fost luată în considerare doar de k factori, în timp ce depinde de al -lea factor Valoarea factorului (k+ ) nu s-a schimbat în experiment dintr-un motiv sau altul Atunci toate punctele simplexului fc-dimensional sunt de fapt puncte ale spațiului (&-|- )-dimensional care se află în hiperplanul xk+t = d, unde d este o valoare fixă a factorului (k+ ) în formă adimensională Din relațiile geometrice rezultă că pentru a construi un simplex de dimensiunea ) dintr-un simplex fc-dimensional, este necesar să se găsească centrul de greutate al punctelor simplexului fc-dimensional în (k + )-dimensional spațiu și trageți o perpendiculară prin acest punct pe hiperplanul care conține punctele simplex dimensionale Dacă pe această perpendiculară este trasat un segment de lungime hk+] (înălțimea unui simplex k+ -dimensional), atunci punctul rezultat împreună cu cele inițiale formează un simplex (k+ )-dimensional Coordonate noi de punct x , ), , x ) r(C) \u d i- ± - \u d ( r(e)= ± ± - , | = > r(e)=r - , + - = gM \u d - , + , - o\u e o Apoi coordonatele celui de-al nouălea punct sunt exprimate după cum urmează: d} ) = - , - , = , , z ( > = - , - , = , , z - - , - , = , , \u d - - \u d , z \u d - , - \u d , , z^ ) = - , - , = , , r = - , - , = , În mod similar, la reflectarea primului punct, s-au obținut coordonatele celui de-al zecelea punct În simplexul , , , , , , , , cel mai rău punct - Reflexia lui dă coordonatele punctului , reflexia punctului -/ dă coordonatele punctului În simplexul , , , , , , Și, cel mai rău punct este punctul Reflexia lui dă coordonatele punctului Ieșirea la punctul este mai mică decât la punctul Reflecția punctului va duce înapoi la punctul - Astfel, simplexul este buclă Definim ieșirea în centrul simplexului , , , , , , , Coordonatele centrului simplexului punctului ,- = , + - , + , + , + , + , = , ; ,(S) Z - , - + , - + , + , + , + , •— , ; ,(S) Z - , • + , + , + , + , , ; ^S> = - + - + + + + = ; (S) , - + - , + , + , + , + , z ' = , ; (S) , - + , - + , + , + , + , z ~ = , ; zP>= - , - + , - + , + , + , + , = , La punctul , experiența a fost implementată Valoarea obținută a densității optice ylS> - - , Astfel, cea mai bună valoare a criteriului de optimizare a fost obținută în centrul simplexului în experimente Metoda de ascensiune abruptă a necesitat experimente pentru a rezolva aceeași problemă Exemplul Am studiat reacția care se desfășoară după schema L + B + C» D într-o soluție apă-alcool Calitatea şi cantitatea produsului D(y) au fost influenţate de următorii factori: -zi — timpul de reacţie, h; gg este conținutul de alcool dintr-o soluție de apă-alcool, spun ei acțiuni; xs este concentrația substanței C, spun ei acțiuni; z este concentrația substanței B, mol acțiuni; zs este raportul molar al substanțelor B și A Principalul nivel și intervalele de variație a factorilor sunt prezentate mai jos Factori g, z z z z z° , , , , , Azy , , , , , Determinați condițiile pentru obținerea cantității maxime de produs (Umax) - Soluţie Să folosim metoda de planificare simplex Pentru k - extragem din matrice (vezi V ) o submatrice care contine cinci coloane si sase randuri Folosind formula de codificare (V ), obținem; zi - -° r - , z - , * , ' x*~ o ' Xs ~ , ' z - , g - , Xi ~ , ' Xb \u d , Atunci matricea simplexului original la scară naturală are forma (tabelul) Experienta numarul z * * * , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , După cum reiese din tabel, experiența este cea mai proastă Să înlocuim punctul cu imaginea sa în oglindă - punctul Găsim coordonatele noului punct folosind formulele (V ) și (V ) Să determinăm coordonatele punctului c — centrul feței formate din punctele , , , , : r(o = !^ ±^> = > ; r(O = - , + ^ + , = q g (O \u d - , + , + - , \u d Q g (O \u d - , + , + , Q , \u d - , - , \u d , : ) \u d - , - , \u d , ; ) \u d - , - , \u d , ; r^ ) = - , - , = , ; zg> \u d - - \u d Noul punct , împreună cu restul, formează un simplex , , , , , (tabel) Numărul experimentului Z] z z Y , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , După experimentul de la punctul , cel mai rău punct al complexului , , , , , s-a dovedit a fi punctul Reflexia lui relativ la fața , , , , oferă condițiile pentru următorul experiment etc După efectuarea celui de-al șaptelea experiment, se adaugă încă un factor zs - numărul de rotații ale agitatorului Până acum, ze a fost menținut la un nivel constant Zj - rpm Atunci coordonatele punctului în formă adimensională vor fi: „(*) g(l) r! s), + h xl > x xk+l g l+ - Aze = rpm este luată ca unitate de variație, iar zj — rpm ca nivel principal Apoi formula de codificare pentru ze ia forma ze- Xe== О xg>= Înălțimea simplexului cu șase dimensiuni se obține prin formula (V ) A- , Să determinăm valorile factorilor pentru experimentul Valorile primilor cinci factori reprezintă coordonatele centrului de greutate al simplexului cincidimensional , , , , , : (și , + - , + , p z = z = z| > = , ; z = z| ) = , ; z| ) = , ; z = + x^ ) = + ( x^> +fte) = rpm Experiența , împreună cu punctele - , formează deja un simplex cu șase dimensiuni , , , , , , (tabel) După implementarea celui de-al optulea experiment, este necesară analizarea rezultatelor și reluarea procesului de reflecție, luând în considerare șase factori Numărul experimentului r z * z Y , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Planuri Plackett-Berman saturate ortogonale Proiectele D-optimale ortogonale saturate pe două niveluri pot fi construite folosind replici fracționate din PFE pentru numărul de factori k = (u= ), k = ( V= ), fc= (/V= ) , fc- (N - ), etc Cu toate acestea, clasa modelelor saturate ortogonale poate fi extinsă considerabil Pluckett și Berman au dezvoltat o teorie matematică riguroasă pentru construirea și analiza planurilor ortogonale În special, s-a dovedit că într-un plan saturat, estimările celor mai mici pătrate ale efectelor au o acuratețe maximă pentru un număr dat de experimente N, care este aceeași pentru toate efectele dacă matricea de planificare are coloane ortogonale Pentru ca matricea să fie ortogonală, este necesar și suficient ca: ) fiecare factor să apară la fiecare dintre nivelurile sale de același număr de ori; ) fiecare doi factori cu orice combinație a nivelurilor lor întâlnite de același număr de ori; ) numărul de experimente a fost împărțit la pătratul numărului de niveluri, i e # = n/ , (V ) unde p este un număr întreg Cu această formulare a condițiilor de ortogonalitate, problema construirii unei matrice ortogonale (plan de experiment) se transformă într-o problemă pur combinatorie Dacă /V - nP, atunci numărul de factori ale căror efecte pot fi calculate Dat N este egal cu l = ("/" ])/(/ ]) (V ) sau o întreagă parte din ea Dacă numărul de niveluri pentru toți factorii este egal cu doi, atunci problema construirii unui plan optim se reduce la construirea unei matrice ortogonale formată din + și - , de mărimea L'XL', unde L' este un multiplu al patru, adică A = n Numărul maxim de factori care pot fi introduși în planificare este egal cu k—N- Pentru a construi planuri bogate pentru fc- , , și , folosim rândurile din Tabel Tabelul pentru k - , , și la N Combinații de caractere + + - + + + + - + + + + + + - + - + + + - + + + + + - + - + + + + + - + -+-+++ + + + + + - + + + + + - + + + - Când construiți un plan pentru = (vezi exemplul ), un rând din tabel este luat ca elemente ale primei coloane Obținem a doua coloană din prima, înlocuind primul element din ea cu ultimul și deplasând toate celelalte elemente în jos în consecință A treia coloană se obține prin înlocuirea primelor două elemente din prima coloană cu ultimele și deplasarea în jos a elementelor rămase etc Elementele ultimului rând sunt Planurile pentru fc = , și sunt construite într-un mod similar Pentru k ^ II, la construirea planificării matriceale se folosesc trei blocuri A, B și C, prezentate în tabel Aceste blocuri sunt scrise în ordinea permutației circulare: ABC TAXI B C A și li se adaugă din nou un șir, toate elementele fiind Tabelul A în cu + - + + + + + + -+ + + -+++++ + +-++++ + + - + + + + + + + - + + + - + + + + + + - + + + + + - + + + + + - + + + + + - + - + + - + - - + + ++-+-++-+ -++++-++- + -+-++-++ +-+++-+-+ ++ ++++ +++-+-++ +-++-+++-++-++ - -++ -++-++++-+ Pluckett și Berman au arătat cât de bogate modele de până la pot fi construite cu N divizibil cu (cu excepția lui N= ) Utilizarea planurilor Plackett-Berman face posibilă obținerea unor estimări separate ale efectelor liniare ale tuturor factorilor cu precizia maximă posibilă pentru un număr dat de experimente, care este aceeași pentru toate efectele Orice coeficient al ecuației de regresie liniară este determinat de formulă N U*L Ui /=l' k- Eroarea în determinarea bj în acest caz este egală cu sb / - sBocnp/l^^> unde £repro este eroarea de măsurare Deoarece matricea de planificare este ortogonală, această estimare a efectelor liniare coincide cu estimarea obținută prin metoda celor mai mici pătrate În plus, datorită ortogonalității matricei, estimările obținute ale efectelor liniare nu sunt amestecate între ele Raportul dintre bj și lvosp/t/N are o distribuție Student pentru ipoteza nulă, adică valoarea adevărată a lui B( = Acest raport poate fi folosit pentru a testa semnificația efectelor Pentru a testa semnificația diferenței dintre efectele, puteți folosi raportul uh - uh (v ) f = resc ^ / V "N având și o distribuție Student Exemplul S-a studiat posibilitatea obținerii de îngrășământ cu azot-fosfor-potasiu prin înlocuirea parțială a potasiului cu amoniac în timp ce se neutralizează extractul de acid azotic Procesul de neutralizare poate fi caracterizat prin reacția globală: Са (NO ) + ЗНЭРО + NH + ЗК СО = ЗСаНРО + NH NO + + KNO + CaCO În studiul neutralizării succesive a extractului cu amoniac și potasiu, a fost de un interes deosebit elucidarea gradului de retrogradare a formelor asimilabile de oxid de fosfor (V) Prin urmare, indicatorul procesului (y) a fost gradul de digestibilitate al compușilor fosforici formați (procentul cantității de forme solubile în apă și solubile în lămâie de fosfor față de cantitatea totală de fosfor din produsele de reacție) Au fost aleși ca factori independenți următorii factori: zi este temperatura de amoniație ( °C); z este durata amoniației ( – min); za este norma amoniacului ( ^ % din norma stoichiometrică); z", zs, ze, z? - continutul de impuritati din extractul original, respectiv , % Mg (NOs) ; +■ , % Fe(NOa) ; - , % AI(MO )e; - - , % HaSiFe; ze este temperatura în timpul interacțiunii componentelor extractului amoniat cu o soluție de carbonat de potasiu ( - - °C); za este durata interacțiunii cu carbonatul de potasiu ( h - min); zio este rata carbonatului de potasiu ( - % din rata stoechiometrică) Conținutul de PrOb ( , %) și CaO ( %) din extract a rămas constant Soluţie Primele coloane ale planului Plackett-Berman au fost folosite ca plan experimental (Tabelul ) Tabelul Matricea de planificare și rezultatele experimentale Numărul experimentului X, x * X * ■* x -* X$ X\o Y + + - + - - - + + + - + + + - + - - - + + + , + - + + - + - - - + + , + + - + + - + - - - + , + + + - + + - + - - - , + + + + - + + - + - - , + - + + + - + + - + - , + - - + + + - + + - + , + - - + + + - + + - , + + - - - + + + - + + , + + - - - + + + - + , + - - - - - - - - - , b , , - , - , , , , , - , - , , 'j , , , , , , , , , , , Valorile medii ale gradului de asimilare y au fost determinate prin două experimente paralele Varianta de reproductibilitate este ^repro = , Numărul de grade de libertate /VOsp " Valoarea tabelară a criteriului Studentului ta, oz ( ) - , Astfel, toți coeficienții ecuației de regresie s-au dovedit a fi semnificativi, iar ecuația are forma y \u d , - , * j - , * - , * + , * + , * + + , * in + , *, - , * in - , *, + , * (V ) Deoarece planul experimental a fost nesaturat (N - , k = ), există un grad de libertate pentru a verifica adecvarea ecuației la experiment Dispersia adecvării este egală cu L ( - Ui) o \u d iad —n = , Valoarea raportului F F = i ad/NV SPr este , Valoarea tabelară a criteriului Fisher Fo (l, ) = , Astfel, ecuația este adecvată experimentului Ecuația de regresie rezultată ne permite să determinăm condițiile care asigură producerea de îngrășăminte pe bază de extract de acid azotic din fosfați, care conțin tot fosforul într-o formă asimilabilă În absența experimentelor paralele, efectele așa-numiților factori imaginari pot fi folosite pentru a estima pvosp Factorii imaginari sunt introduși dacă designul este nesaturat, adică k A+ , bk+ , bN i~ efectele factorilor imaginari, atunci cantitatea = [b +i + ^+ + • • • + bN-\ servește ca estimare a erorii pătratice în determinarea efectelor factorilor principali și ^ref= N ( bk+l + bk+ + • • • + bN~l ) X k i este o estimare a erorii pătrate în măsurarea răspunsului, presupunând că nu există efecte de interacțiune Planurile Plackett-Berman sunt în unele cazuri planuri mai economice decât replici fracționate pentru o anumită dimensiune a spațiului factorilor Experimente de screening Metoda de echilibrare aleatorie Pentru a reduce numărul de experimente, valorile unor factori sunt adesea stabilizate în timpul procesului de cercetare fără motive suficiente La hotărâre problema de optimizare, aceasta duce la determinarea doar a extremelor locale ale procesului Pentru sarcini multifactoriale, experimentele de screening sunt efectuate în prima etapă a studiului Întrucât intensitatea influenței unui factor este legată de gama modificării acestuia, mulți factori care sunt suspectați a fi semnificativi pe baza unor informații a priori se pot dovedi a fi nesemnificativi Prin urmare, experimentele de screening sunt eficiente nu numai în studiul noilor procese, ci și ca primă etapă în studiul proceselor multifactoriale cu suficiente informații a priori, dacă numărul de factori este prea mare pentru a planifica imediat un experiment care să vizeze găsirea optimă condiţiile procesului Pentru a filtra factorii cantitativi și calitativi cu un număr de niveluri egal cu două, se pot folosi replici fracționale dintr-un experiment factorial cu un grad suficient de ridicat de fracționalitate, precum și planuri Plackett-Berman ortogonale saturate Aceste planuri fac posibilă obținerea unor estimări separate ale efectelor liniare ale tuturor factorilor cu precizia maximă posibilă pentru un număr dat de experimente, care este aceeași pentru toate efectele Acesta din urmă este deosebit de valoros în etapa de screening, deoarece nu se știe care efecte vor fi semnificative Dezavantajele acestor planuri includ cerința pentru absența efectelor semnificative de interacțiune Pentru a izola efectele semnificative - interacțiuni liniare și perechi - Sutterswyth a propus o metodă de echilibrare aleatorie In aceea metoda, se propune ca planul de experiment să fie suprasaturat - numărul de experimente N din matricea de planificare este mai mic decât numărul de experimente luate în considerare efecte, adică la începutul studiului, numărul de grade de libertate / X X + ^ X X + • • • + & £ - , k xk-l xk + ®play> unde sprob este eroarea experimentală Presupunând că unele Xj denotă interacțiuni de pereche, modelul divizat poate fi scris într-o formă liniară: Y = b<>+Y,b}x}+'Y C}x + u = ba+'^b}x} + a, /=i /=i /=i unde I-t este numărul de efecte semnificative; m este numărul de efecte legate de câmpul de zgomot; T „ = c> x' , ; Gv \u d ^ - \u d \u e Astfel, ambii coeficienți sunt semnificativi Apoi, efectele factorilor semnificativi sunt eliminate pentru a identifica efectele mai slabe Pentru a face acest lucru, valorile bj sunt scăzute din toate valorile lui y pentru care factorul Xj a fost la nivelul + Deci, la eliminarea efectului factorului, valoarea - , a fost scăzută din y , y , y , uv, yI , y , y , V , iar la eliminarea efectului factorului x , valoarea x , a fost scăzut din y , y , Uv- Yi-Ub-Tio-Tiz> Tis- Rezultatele obţinute ale valorilor lui y' sunt date în tabel Apoi se repetă întreaga prelucrare a datelor experimentale, începând cu construirea unui grafic de dispersie pentru valorile y corectate Criteriul pentru sfârșitul screening-ului este raportul de dispersie F = ®av/ ®aosp ' unde *, este varianța rezultatelor în raport cu media aritmetică a acestor rezultate În exemplul luat în considerare, întregul ciclu a trebuit să fie repetat de șapte ori; efectele selectate și estimările lor sunt date în Tabelul Tabelul Tabelul rezumativ al rezultatelor experimentului de screening Numărul etapei Efecte identificate Estimarea coeficienților de regresie Dispersia c sr Numărul etapei Efecte identificate Estimarea coeficienților de regresie Dispersia sr I x, , , V x| - , , X * * ~HJ II Х Х , , VI , , xbx III - , , VII x - , - , , IV XeXio Dispersia de reproductibilitate F = , / , = , mai mic decât tabelul F ; ( , ) = , Prin urmare, diferența dintre varianțe și nu este semnificativă Prelucrarea manuală a rezultatelor prin metoda echilibrului aleator este extrem de laborioasă Se propune un algoritm de procesare a rezultatelor unui echilibru aleator pe un computer digital, așa-numita „strategie de ramificare” A fost dezvoltat un algoritm pentru a evidenția cele mai mari efecte din diagramele de dispersie Această etapă nu introduce nimic nou în comparație cu prelucrarea manuală Pentru cuantificarea efectelor identificate, se utilizează analiza de regresie convențională Până la douăzeci de coeficienți de regresie pot fi evaluați simultan În această etapă, se face deja o îmbunătățire semnificativă Dacă în loc de trei, douăzeci de efecte sunt evaluate simultan, atunci varianța reziduală scade brusc și, prin urmare, sensibilitatea metodei crește Selectarea efectelor semnificative se realizează în două etape În primul rând, efectele care diferă de zero cu mai puțin de Zu sunt eliminate; în caz contrar, eliminarea ulterioară strategia de ramificare ar fi o sarcină greoaie Efectele rămase sunt apoi eliminate cu o strategie de ramificare; valoarea varianţei reziduale este utilizată ca criteriu Semnificația strategiei de ramificare este că efectele sunt eliminate secvenţial în toate combinaţiile posibile câte unul câte unul, câte două și așa mai departe Acest proces secundar de screening continuă până când este selectat grupul de efecte care asigură variația reziduală minimă Dacă luăm în considerare corelarea coeficienților de regresie (datorită neortogonalității matricei de planificare), atunci este clar că strategia de ramificare are avantaje neîndoielnice față de evaluarea separată a coeficienților de regresie folosind criteriul / Evaluarea separată evidențiază adesea efecte false Întregul proces de cernere se repetă de mai multe ori până când variația reziduală este redusă la valoarea dorită Planificarea unui experiment în determinarea constantelor ecuațiilor cineticii formale În prezent, se dezvoltă intens o nouă direcție în aplicarea metodelor statistice pentru studierea mecanismului și determinarea constantelor cinetice ale reacțiilor chimice complexe Să luăm în considerare cele mai simple tehnici bazate pe utilizarea ideilor și metodelor de planificare a experimentelor extreme pentru a determina constantele ecuațiilor cinetice formale Cea mai răspândită metodă de prelucrare a datelor cinetice, care constă în liniarizarea dependențelor cinetice folosind transformări speciale De exemplu, viteza de reacție t p (V ) i=i /=i unde Aj și Bj sunt materii prime și respectiv produși de reacție; a și bj sunt coeficienți stoichiometrici; poate fi descris prin ecuația formală r = /(P • (V > /=i unde K este constanta vitezei de reacție; p, este ordinea reacției față de substanța i-a; [А, ] este concentrația substanței /-a Dacă dependența constantei de viteză de temperatură respectă legea Arrhenius, atunci avem A \u d Ko exP I -\u e \RT) (v ) unde Ko este factorul pre-exponențial; t este energia de activare; L este constanta gazului; Г este temperatura absolută Aritmarea loi a ecuației (V ) în vederea (V ) dă pG = pK - + p p[Ai] (V ) Notând variabilele astfel: L În l / r \u d Z , În [L ] - Z , Pi \u d a , ІП ІАп] - ^m + i, obținem o ecuație liniară L U = „o t+ + Ș аІгІ- i= (V ) (v ) Pentru a determina coeficienții ecuației (V ) prin metoda de planificare a experimentelor, se pot folosi planuri liniare ortogonale cu numărul de experimente m - Exemplul Cinetica procesului de copolimerizare a a, p, p-trifluorostiren cu acid metacrilic a fost studiată la grade scăzute de conversie a monomerilor iniţiali într-un polimer Ecuația cineticii formale a procesului a fost luată în considerare sub următoarea formă: UZ = K[Cl + C' ]"*[/Z]"*> (V ) unde [Сі] și [Сі] sunt concentrațiile monomerilor inițiali, mol/l; [//]—concentrații ale inițiatorului (peroxid de benzoil), mol/l; i, și P sunt ordinea reacției în ceea ce privește concentrația totală de monomeri și inițiator; E este energia de activare, kJ/mol; K este constanta de viteză totală a reacției Ecuația (V ibb) face posibilă determinarea ordinii de reacție în raport cu concentrația totală de monomeri și inițiator, energia de activare și constanta vitezei de reacție Constanta vitezei de reacție K depinde de temperatură conform ecuației (V ) În experiment, concentrația totală a monomerilor inițiali a variat în intervalul , -t - - , mol/l, concentrația inițiatorului - în intervalul , si = Ă [A] + k [B] [C], ag unde k/(i - , , , ) sunt constantele de viteză ale reacțiilor pseudomonomoleculare; [A], [B], [C], [D] sunt concentrațiile de reactivi în faza gazoasă Astfel, este necesar să se determine constantele vitezei a nouă reacții din datele experimentale și să se arate că sistemul (V ) descrie în mod adecvat cinetica reacției Soluţie Să notăm concentrațiile inițiale ale substanțelor A, B, C zi, zi, respectiv za Considerând ratele de acumulare a substanțelor ca răspunsuri y și concentrațiile inițiale ca factori independenți, obținem din (V ) Ui = b u+blux +b ux + b ux , u= , , , - (V ) Variabilele хі, х%, хз sunt legate de concentrațiile inițiale prin transformarea liniară obișnuită *і = (zj - ?)/Dy Termenul liber din ecuația (V ) este legat de termenii liniari: zY z, zY ba = b, + bg + b - °“ c Дгі „Az s“ Az (v ) Cinetica dehidrogenării butenei a fost studiată la diferite temperaturi Pentru fiecare temperatură a fost compilată propria matrice de planificare Să luăm în considerare un experiment cinetic la – K Aria de studiu este dată în tabel (z este concentrația de butene, fracții molare): , , , , , , Pentru fiecare funcție de răspuns (V ), trebuie găsiți doar trei coeficienți Prin urmare, atunci când planificați, este suficient să luați din PFE Aici, totuși, este utilizat PFE de tip (Tabelul ), deoarece acest lucru permite să se determine mai precis constantele vitezei de reacție și să se analizeze mai detaliat abaterile valorile prezise de ecuațiile din datele experimentale Pentru a determina varianța de reproductibilitate, unele experimente au fost duplicate Tabelul Matrice de planificare pentru studierea cineticii dehidrogenării buteiovului Numărul experienței -* x •* > Oz - - - , , , , - - + - , , - , , - + -I - , - , , , + - - - , , , , + + + - , - , - , , + + -I - , - , , , + - + - , , - , , - + + - , - , - , , Viteza de acumulare a substanțelor y (fracții mol /s) a fost determinată prin formulele: [A]t - [A], [Bjx - [B], Y = -, yi = , t z (v ) [С]х- [С]o [D]ț-[D] Uz \u d - Ui \u d • t t Aici, indicele „ ” se referă la concentrațiile inițiale de butene, iar indicele m se referă la concentrațiile corespunzătoare timpului t Astfel, segmentele inițiale ale curbelor cinetice au fost aproximate printr-o linie dreaptă În acest caz, estimările calculate ale vitezei de acumulare a reactivilor vor fi părtinitoare, iar părtinirile depind de viteza de reacție Cu cât reactantul este mai activ din punct de vedere chimic, cu atât polarizarea va avea constantele de viteză corespunzătoare Tabelul prezintă valorile constantelor și erorile acestora, obținute din datele experimentului planificat și din metoda estimărilor neliniare (INR) Constantele ratei, s~' / ' kі Ă k* kъ Zs k? Zs kts Planificat Conform INR , , , , , , , , , s- (c) , , , , , , , , , Compararea constantelor de viteză cu erorile lor arată că un număr de constante nu ies în evidență pe fundalul zgomotului Pentru a reduce erorile constantelor, este necesară creșterea intervalelor de variație Estimările constantelor obţinute au fost rafinate prin metoda estimărilor neliniare (INE) Conform acestei metode, constantele vitezei de reacție trebuie alese în așa fel încât suma abaterilor pătrate să fie minimă (V ) Concentrațiile se obțin prin integrarea sistemului (V ) de la t - la T - m în condiții inițiale (vezi tabelele de la p ) Însumarea a fost efectuată pentru toate experimentele, iar termenii au fost introduși cu mase egale, deoarece s-a dovedit că erorile în reproductibilitatea concentrațiilor tuturor substanțelor sunt omogene Constantele determinate conform planului au fost folosite ca aproximare inițială Apoi, conform criteriului Fisher, a fost verificată adecvarea modelului matematic (V ) la experiment: F (^) min S°CT“ , - ' unde φ( Noi medii y, - = y * Calculul efectelor Calculul intervalelor de încredere (IC) Efect acid#! - y (ѵz + st - y? - uz) = - , Efect de încărcare | (uz + y - y - y) \u d - , Efectul de interacțiune D \u d I (uz + uz -ud - Uz) - , Efectul modificării mediei M - y (uz + uz + beat + ys - yi) = , CI pentru mijloace noi - \u d s y \u d + , CI pentru efecte '/'r =± , V'n D CI pentru efectul modificării medie = , y p D Concluzii: nu s-au constatat efecte semnificative Eroare a priori a experienței zy = , Soluție: repetați experimentul În blocul „Calculul efectelor” sunt luate în considerare efectele factorilor și interacțiunile acestora Efectul unui factor este egal cu creșterea valorii parametrului de optimizare (răspuns) ca urmare a trecerii factorului de la nivelul inferior la cel superior În cazul a doi factori, efectul lui x efect x „Uz + Ui - Ug - Ui = n Uz+ Ub—Ug—Ui P (V ) (V ) efect de interacțiune Blt = y? + uz (V ) Astfel, efectele factorilor sunt egale cu dublul coeficienților ecuației de regresie atunci când se utilizează modele liniare ortogonale, cum ar fi k sau k~p (vezi formula V ) Rezultatele calculelor B}, B și B] sunt introduse în rândurile , și din blocul „Calculul efectelor” În a -a linie a acestui bloc se înregistrează valoarea efectului modificării parametrului mediu de optimizare (c) Acest efect este calculat ca diferența dintre media generală pentru toate rezultatele experimentului și valoarea lui y} obținută în centrul planului (la nivelul principal): efectul de schimbare medie Y r + Uz + Ui + Ub + Ui - Ug + Uz + Ui + Uz - tyi „, , I* \u d -~ Ui \u d - (v ) Efectul modificării mediei caracterizează suprafața de răspuns în zona de experimentare Dacă suprafața are un extremum, iar experimentele sunt efectuate departe de acesta, efectul modificării mediei va fi semnificativ, iar dacă q este semnificativ mai mic decât zero, suprafața are un maxim, dacă q este semnificativ mai mare decât zero, are un minim Dacă extremul este plat sau există un punct de șa, efectul modificării mediei în regiunea extremului va fi aproape de zero În blocul „Calculul intervalelor de încredere”, intervalele de încredere (IC) sunt calculate pentru medie, efectele și efectele modificării mediei Pentru a construi intervale de încredere, se utilizează distribuția lui Student Pentru efectele factorilor și efectele de interacțiune "v CI = ± t ■; (v ) Da N - Da p pentru efectul de modificare medie Când L' = + = pentru probabilitatea de încredere B= , valoarea t « și, prin urmare, din (V ) și (V ) avem su CI=± —, (V ) Ambalaj CI^ = ± , — (V ) Sf Efectele sunt considerate semnificative statistic dacă sunt egale sau mai mari decât valoarea absolută a intervalelor de încredere În exemplul luat în considerare (vezi Tabelul ), eroarea experimentală a priori sy = , a fost utilizată pentru a calcula intervalele de încredere Compararea dimensiunilor efectului cu limitele intervalelor lor de încredere indică despre absența efectelor semnificative și necesitatea de a trece la cel de-al doilea ciclu al acestei faze, adică să repeți din nou toate experimentele planului Masa de lucru se umple după al doilea ciclu În al doilea ciclu, începem să primim informații despre eroarea experimentului din datele aceleiași faze Secvența de completare a foii de lucru pentru al doilea ciclu este următoarea Datele din rândurile și din blocul „Calculul mediilor” din tabelul primului ciclu (vezi Tabelul ) sunt transferate în rândurile și ale aceluiași bloc din tabelul celui de-al doilea ciclu (Tabelul ) ) Tabelul EVOP, faza , ciclul l = Calculul mediilor Calculul erorii de experiență Suma numerelor anterioare ale punctelor de experiment Suma erorilor anterioare , , , , , CICLURI DE CICLURI = Media anterioare Eroarea medie a anterioare total cicluri , , , , , total cicluri = Observații noi , , , , , Interval = , Diferențe ( ) - ( ) - , , , , - , Nouă valoare de eroare Sy = span „ = , Sume noi ( ) + ( ) , , , , , Sumă nouă de erori n = , Mijloace noi Eroare medie nouă p , , , , , / i '-eu!' Calculul efectelor Calculul încrederii intervale (CI) Efectul acidității B\ \u d / (nod + nod ~ Oz - nod) - - , Efect de încărcare Vz \u d / (kn + kn - kn - U ) \u d - , Efectul de interacțiune Vіz \u d / (uz + uz - st - uz) \u d - , Efectul modificării mediei М = '/s (уі + Uz + U + uz - уі) - - , CI pentru mijloace noi CI pentru efecte ^-= , CI pentru efectul modificării I medie 'sv - , - , /G/T Concluzii: nu s-au constatat efecte semnificative Eroare anterioară sy - Soluție: repetați ciclul Rezultatele experimentelor obținute în al doilea ciclu sunt înregistrate în rândul În rândul , diferențele rândurilor și sunt trecute cu semnele lor, cele mai mari și cele mai mici diferențe sunt subliniate Acestea sunt diferențele între rezultatele medii ale ciclului anterior și rezultatele noi În rândul al blocului „Calculul mediilor”, notați noile sume ale rezultatelor observației pentru fiecare punct al planului, adică sumele rândurilor și În rândul , notați noile medii pentru fiecare punct al planului ca coeficienti din împărțirea noilor sume (linia ) la numărul de cicluri la un timp dat, adică, în cazul nostru, coeficientul din rândul împărțit date prin n - După formulele (V ) - (V ) calculați efectele factorilor, interacțiunile acestora și efectul modificării mediei Valorile obținute sunt introduse în blocul „Calculul efectelor” În blocul „Calcul de eroare din experiment” pentru al doilea ciclu, rândurile și nu sunt completate; în rândul se scrie intervalul de valori Ay, definit ca diferența dintre cele mai mari și cele mai mici valori ale diferențelor Ay subliniate în rândul din blocul „Calculul mediilor” În cazul nostru, intervalul este Au \u d , - (- , ) - , În rândul , notați valoarea noii erori a experimentului, care este considerată prin metoda intervalului conform formulei (V ) Valoarea lui fNn este găsită din tabel În cazul nostru, L \u d și \u d , / \u d , Valoarea erorii = , - , = , Până la sfârșitul celui de-al doilea ciclu, încă nu există estimări ale erorii experimentale din ciclul anterior, astfel încât rândurile și din blocul „Calcul de eroare experimentală” sunt rescrise cu valoarea noii erori (linia ) Conform formulelor (V ) - (V ), se determină intervale de încredere pentru noua medie, efecte și efectul modificării mediei, care se înregistrează, respectiv, în rândurile , și din blocul „Calculul intervalelor de încredere” Comparați dimensiunea efectelor cu intervalele de încredere corespunzătoare și concluzionați că efectele sunt semnificative din punct de vedere statistic, apoi decideți asupra pasului următor În cazul nostru, toate efectele luate în considerare în al doilea ciclu s-au dovedit a fi nesemnificative din punct de vedere statistic și s-a decis să se repete toate experimentele planului încă o dată, adică să se treacă la al treilea ciclu Începând cu al treilea ciclu, toate rândurile foii de lucru sunt completate (Tabelul ) Toate calculele se efectuează în același mod ca pentru ciclurile anterioare, doar în blocul „Calcul eroare experiment” din rândurile și sunt rescrise rezultatele de la rândurile și ale aceluiași bloc al ciclului anterior În rândul se introduc rezultatele obținute prin însumarea liniilor și , rândul se obține prin împărțirea liniei la n - Evaluarea rezultatelor celui de-al treilea ciclu indică o diferență semnificativă față de zero a efectului acidității În această situație, se ia una dintre cele două soluții posibile Continuați faza actuală de planificare pentru a evidenția mai bine efectele semnificative Începeți o nouă fază de planificare Începând o nouă fază, ca centru al noului plan, puteți alege punctul planului din faza anterioară cu cea mai bună valoare a funcției obiectiv În exemplul luat în considerare (Tabelul ), s-a decis începerea unei noi etape cu centrul planului la punctul (Fig ) În noua fază de planificare, dacă este necesar, puteți schimba Tabelul EVOP, faza , ciclul l = Calculul mediilor Calculul erorii de experiență Numărul de puncte experimentale Suma precedentului Suma erorilor precedentului cicluri , , , , total CICLURI = , Media anterioare Eroarea medie a anterioare total cicluri , , , , , total cicluri sy = , Observații noi , , , , , Interval - , Diferențe ( ) - ( ) , - , , , Noi valori de eroare Sj,-span Dg „- Sume noi ( ) + ( ) , , , , , Sumă nouă de erori J „= , Mijloace noi Eroare medie nouă - p , , , , , b=-T?T- , Calculul efectelor Calculul încrederii intervale (CI) Efectul acidității Bi \u d (uz + U - y - Uz) \u d - , Efect de încărcare Bg \u d -y (uz + ' - y - Y ") - , Efectul de interacțiune V \ g \u d -Ș- (yi + uz - U - uz) - , Efectul modificării mediei •O \u d (yi + uz + y + uz - yi) \u d , CI pentru mijloace noi —- , (V ) „*■ ? - ? • Tabelul Matricea de planificare Numărul experienței X X X X X X X X X X X X Y > - I - + + + - Y + + - + - - - bloc Oz + - + - + - - St - + + - - + - J V - - + + - + Y? II + + + + + + + Te bloc + - - - - + + X - + - - + - + D' Experimentele blocului II constituie o altă semireplică ~ cu un contrast definitoriu =x^x x Sistemul de notare este următorul: I}' ->■ + ? • B "-? + Rez (V ) ÎN? -*■ Z + S' Combinarea experimentelor blocurilor I și II duce la un experiment factorial complet , care face posibilă obținerea de estimări nemixte pentru efectele de interacțiune liniare și de pereche: Яі = ~ (SJ + SV) - ?i ^+s") ■* fi = ~ ("g + g) ₽ (V ) ^ = ~ [B^ Bj) [iu, \u d y (S "-S H B =~ (Bj' -B, ) -> ? Să luăm în considerare procedura de completare a meselor de lucru într-o problemă cu trei factori folosind exemplul unui experiment industrial pentru a minimiza pierderea de azot legat S-a decis începerea fazei a doua a experimentului cu centrul la punctul (vezi Fig ) Să introducem al treilea factor, temperatura (x ), în considerare Nivelul principal, intervalele de variație și limitele zonei de studiu sunt date în tabel Tabelul Factori Aciditate, sarcină G/l NHg, m /h Temperatura, °C Nivelul principal Intervalul variază , Xj =+ , , În urma a trei cicluri ale fazei a doua, s-au obținut următoarele valori de răspuns \u b\u b, g / l (Tabelul ) Tabelul Blocul I II Experimentează numerele de puncte Ciclul , , , , , , , , , , Ciclul , , , , , , , , , , Ciclul , , , , , - , , , , Calculele pentru fiecare ciclu sunt efectuate în două foi de lucru, separat pentru fiecare bloc În tabel prezintă tabelul de lucru al blocului I al celui de-al treilea ciclu al celei de-a doua faze a unui experiment industrial pentru a minimiza pierderea de azot legat în producția de azotat de amoniu și în tabel - blocul II al aceluiași ciclu al aceleiași faze În blocul „Calculul efectelor” (Tabelele , ), estimările efectelor sunt calculate în conformitate cu (V ) și (V ) În blocul „Calcul final al efectelor” din fișa de lucru pentru blocul II, estimările pentru efecte sunt mediate în conformitate cu (V ) Aici, efectul modificării mediei în funcție de datele ciclului complet este calculat prin mediarea efectelor corespunzătoare pentru blocurile I și II Calculul erorii experimentale se efectuează în funcție de datele fiecărui bloc, în timp ce intervalul este determinat ca suma dintre cele mai mari și cele mai mici Tabel EVOP, faza , ciclul l - , bloc I Calculul mediilor Eroarea de calcul a experienței Experimentează numerele de puncte Suma anterioare Suma erorilor anterioare cicluri ale blocului I , , , , , cicluri (toate blocurile) - , Media anterioare Eroarea medie a anterioare total cicluri ale blocului I , , , , , total cicluri (toate blocurile) sy — , Observații noi , , , , , Interval - , ; - , Diferențe ( ) - ( ) , - , , , - , Noua valoare de eroare s y - interval fN,n - , Noi sume ale blocului I Noua sumă de erori (toate ( ) -( ) , , , , , blocuri) s - , Mijloace noi Eroare medie nouă bloc I / l p , , , , , s? n s , Calculul efectelor Efectul acidității c} \u d ~ (uz + yi - y - y) - - , Efect de încărcare ÎN? \u d -y (uz + yi, - y - uz) - - , Efectul temperaturii Vz \u d (yi + uz - yi - y\u e) - - , Efectul schimbării blocului mijlociu I M \u d -y (U + Uz + U (V ) /I pentru efecte CI \u d , "y, (V ) /l pentru efectul de modificare medie CI = , —— „y Vl (V ) - Tabel EVOP, faza , ciclul l - , bloc Calculul mediilor Eroarea de calcul a experienței Experimentează numerele de puncte Suma celor anterioare cicluri , , , , , Media precedentului Cicluri NV , , , , Observații noi , , , , , Diferențele ( ) - ( ) -D C) , " , Sume noi ( ) + ( ) , , , , , Medii noi , , , , , Suma erorilor ciclurilor anterioare (toate blocurile) s “ Eroare medie a ciclurilor anterioare (toate blocurile) sy - , Interval - , Noua valoare a erorii sy este intervalul Drn = , Sumă nouă de erori (toate blocurile) s - , Noua eroare medie Calculul efectelor Calculul intervalelor de încredere (IC) Efectul acidității # U= (uv + - nod - th) - - , Efect de încărcare Y P \u d ~ y (SW + YYU „U?” „U > „+ , Efectul temperaturii Vz ^ \u d (uz + ys - uv - th) - , Efectul modificării mediei Mg \u d -y- (y? + uv + uv + uy - uv) - , CI pentru mijloace noi CI pentru efecte , \—F \u d sy- , 'u tu u' CI pentru modificarea efectului mediu , ~~y=sy = , V p Calculul final al efectelor B' = -Ț (v} + I, ) - , ; Vz \u d "(B ' + B ^> ~ " , ; Ve \u d -y- (yz + - „ , ; Efectul modificării mediei Y \u d C-GyzP-dz) - - , Bіz \u d "(Y -Bz) - + , \u d - (B' ~ Bb-" , m \u d ~ ( * + g) ) - , Concluzii: efectul aciditatii este semnificativ Soluție: începe o nouă fază În exemplul luat în considerare, după al treilea ciclu al celei de-a doua faze (Tabelul ), efectul Wu s-a dovedit a fi semnificativ statistic ( , > , ), așa că s-a decis trecerea la o nouă fază Astfel, metoda de planificare evolutivă face posibilă efectuarea unui experiment activ în condiții industriale fără a încălca regimul de producție, în timp ce nu este nevoie de un echipamentelor și personalului de cercetare, toate calculele necesare pot fi efectuate cu ușurință de către personalul de întreținere al unității de producție Capitolul prezintă metode de planificare a unui experiment care vizează găsirea unui extremum în condițiile în care mecanismul procesului este necunoscut În aceste cazuri, se utilizează un model de proces polinomial Alegerea planului este determinată de formularea problemei de cercetare Fiind suficient de departe de extrem, cercetătorul pune la cale experimente pentru a se apropia de condițiile optime Pentru a rezolva această problemă, se folosesc planuri liniare ortogonale Un model liniar este utilizat pentru a determina gradientul în metoda suprafeței de răspuns la urcare abruptă Modelele simplex pot fi, de asemenea, folosite pentru a se deplasa către un extremum Pentru a descrie zona apropiată de extremum, „zonă aproape staționară”, puteți folosi planurile compoziționale de ordinul În multe studii experimentale, împreună cu factorii cantitativi, este, de asemenea, necesar să se modifice factorii calitativi În astfel de situații, planurile complexe se dovedesc a fi foarte eficiente - experimentul factorial A, combinat cu pătratele latine Pentru sarcini multifactoriale, experimentele de screening sunt efectuate în prima etapă a studiului Pentru a rezolva această problemă, pot fi folosite replici fracționale dintr-un experiment factorial cu un grad ridicat de fracționalitate, planuri Plackett-Berman și metoda echilibrului aleatoriu Pentru a căuta optimul în condiții industriale, se folosește metoda de planificare a experimentelor evolutive, cu ajutorul căreia se rezolvă problema izolării unui semnal slab util pe fundalul zgomotului Exerciții Ce proprietăți are matricea de planificare a unui experiment factorial complet *? Construiți matrice de planificare pentru k - și k - Construiți desene ortogonale liniare pentru k - folosind replici fracționale din experimentul factorial complet și desene simplex Ce proprietăți au aceste planuri? Experimentul factorial complet , prezentat în tabel: a fost folosit pentru a studia dependența raportului dintre formele apoase și totale de PaOs (g %) de temperatura procesului de amonizare (zi, °C) și conținutul de apă în faza alcoolică (z , %) în producerea de fosfat monocalcic prin descompunerea acidă a fosfaților prin extracție lichidă Fiecare experiment a fost repetat de ori * a) Determinaţi ecuaţiile de regresie la scări adimensionale şi naturale; b) evaluează semnificația coeficienților; c) verificați adecvarea ecuației de regresie la experiment În timpul maturării superfosfatului obținut din apatită, s-a studiat dependența coeficientului de descompunere (y, %) de timpul de păstrare (zi, zile) la diferite conținuturi de HaSO (Z , părți de masă) Zona de studiu este prezentată în tabel Z\ Z min zj /shah ZJ A fost utilizat un plan ortogonal de ordinul al -lea (Tabelul ), s-au obținut următoarele valori ale coeficientului de expansiune: Numărul experienței y, Numărul experienței Y, % Varianta de reproductibilitate s Bocnp este de %,/reproductibilitatea „ Determinați: a) ecuația de regresie de ordinul doi; b) semnificația coeficienților și adecvarea ecuației de regresie la experiment; c) tipul suprafeţei; d) condiţii optime care asigură descompunerea aproape completă (y - %); e) analizați sensibilitatea factorului de descompunere la modificările timpului de păstrare și ale conținutului de HjSO Conform condiţiilor din exemplul pentru a studia cinetica copolimerizării: a) se face O-optimal simplex; b) folosind rezultatele experimentului (Tabelul ) și proprietățile planului simplex, începeți să mergeți spre creșterea vitezei de copolimerizare CAPITOLUL VI PLANIFICAREA EXPERIMENTULUI ÎN CÂND SE STUDIAZĂ COMPOZIȚIA DIAGRAMEI - PROPRIETĂȚI Metoda rețelelor simplex Când se studiază proprietățile unui amestec care depind numai de rapoartele componentelor, spațiul factorilor este un simplex regulat (q - )-dimensional Pentru sisteme, relația eu *i=i este concentrația componentei; g este numărul de componente Pentru sistemele cu două componente, simplexul - conținutul direct al componentelor este determinat de raportul segmentelor Pentru q - simplexul regulat este un triunghi echilateral Fiecare punct al triunghiului corespunde unei compoziții specifice a sistemului ternar și invers, fiecare compoziție este reprezentată de un punct specific Compoziţia poate fi exprimată în moli, fracţii sau procente în greutate sau volum Vârfurile triunghiului corespund substanțelor pure, laturile sistemelor binare Coborând înălțimea de la fiecare vârf al triunghiului, împărțind fiecare dintre ele în zece segmente de dimensiuni egale și trasând linii drepte paralele cu laturile triunghiului prin diviziunile rezultate, obținem o grilă triunghiulară Apropierea de la fiecare latură dată de vârful opus corespunde unei creșteri proporționale a conținutului componentei corespunzătoare, prin urmare, liniile drepte paralele cu această latură, atunci când se deplasează succesiv de la o dreaptă la alta, reflectă o creștere a celei de-a treia componente de L % CU,% Orez Triunghiul de concentrare Gibbs Orez Triunghiul de concentrare Roseboom % Pe fig , la înălțimile corespunzătoare ale triunghiului, este indicat conținutul fiecăreia dintre componente (în punctul = X - XjX - x x , (VI ) xg \u d x - x, x - x x - Să substituim (VI ) și (VI ) în (VI ) și să prezentăm termeni similari: Y - (&o + + &ii) xi + (^o + ^ + bm) x + (b" + b + baz) xa + + (bі - bn - bm) xgxg + (bіz - bc - baz) xtx + (b - b - b ) x x (VI ) Denota Рі = feo + bІ + bii> Vij = bij - bn - bjj (VI ) Apoi obținem un polinom redus de gradul doi în trei variabile: Y == ? * + ? * + ? * + ? * * + ? S* * * ? S* * - (VI ) Astfel, numărul de coeficienți a scăzut de la la Polinomul redus de gradul doi în q variabile J= ?iXf+ Vіііхіхі (VI ) l *i* (*i - XJ) + l + " * * (" „*h)” + + « * * (* - *h)g + ? * * * + ? *l*f * + P *l* *| (VI • ) pentru componenta ^ Y \u d + ?i)хіхі + Wi (*і - + i Numărul de coeficienți necunoscuți din (VI ) este același ca și din (VI ) Polinoamele de mai sus au primit o aplicație mai largă și pe viitor le vom folosi doar pe acestea Astfel, minimul numărul de puncte experimentale pentru determinarea coeficienților unui polinom de grad n în q variabile este SDN I (Tabelul ) Tabelul Număr componente Gradul polinom Număr componente Gradul polinom (incomplet) (incomplet) Planuri Scheffe simplex-zăbrele În prezent, planurile simplex-zăbrele propuse de Scheffe au primit cea mai mare utilizare Aceste planuri oferă o răspândire uniformă a punctelor experimentale pe simplexul (q - )-dimensional Punctele experimentale reprezintă rețeaua {q, u} pe simplex, unde q este numărul de componente ale amestecului; și este gradul polinomului Planurile simplex cu zăbrele sunt planuri saturate Pentru fiecare componentă, există (u + ) niveluri egal distanțate xz - , /n, /n, , și sunt luate toate combinațiile posibile cu astfel de valori ale concentrațiilor componentelor Deci, de exemplu, pentru o rețea pătratică {q, }, care oferă o aproximare a suprafeței de răspuns prin polinoame de gradul doi (u = ), trebuie utilizate următoarele niveluri ale fiecăruia dintre factori: , / și , pentru un cubic (u - ) —O, / , / și etc Unele rețele { , u} sunt prezentate în Fig , și { , și} - în fig Aceste planuri sunt parțial compoziționale O rețea cubică incompletă { , *}, de exemplu, poate fi obținută din { , } prin adăugarea unui singur punct în centrul simplexului, rețeaua { , } prin adăugarea de puncte la rețeaua , } După ce au scris coordonatele punctelor rețelei simplex, obținem matricea de planificare Să construim matrici de planificare pentru rețelele { , ), { , } și { , } Indicii proprietăților amestecului indică conținutul relativ al fiecărui component din amestec De exemplu, amestecul (Tabelul ) Tabelul Matricea de planificare pentru rețeaua { , } N Xi X xs bks N Xi X xs bks Yi / / Y Y Uz / / Z^іz Uz / / T constă numai din componenta xt, proprietatea acestui amestec se notează y, amestecul este format din '/sx, iar Vrx^, proprietatea se notează y, Matricea de planificare pentru rețeaua simplex { , } este dată în tabel În tabel prezintă matricea de planificare pentru construirea unui polinom de gradul al patrulea într-un sistem cu trei componente Orez { , n}-rețele pentru polinoame: a, de ordinul doi; b - ordinul III incomplet; c - ordinul al treilea; r-al patrulea ordin Orez { , nTrete pentru polinoame; a - ordinul doi; b - ordinul III incomplet; în ~ al treilea ordin; d - ordinul al patrulea Tabelul N Xi X X Y N Xi X X y Ti / / U U / / U Tz / / J' / / U / / L'izz / / U / / / U' I tabelul Matricea de planificare pentru , }-latice Xi X He Y Xi X X Y / / U U / / J' Tz / / V' / / U / / J/ / / Uіz / / / U / / Y / / / J/ / / UіІі / / / U / / U Coeficienții polinoamelor reduse se obțin folosind proprietatea de saturație a planului Pentru a obține coeficienții polinomului U - * + + Z * + * * + * * + Z * -'% vom substitui secvenţial coordonatele tuturor celor şase puncte ale matricei de planificare (Tabelul ) în ecuaţie Apoi, înlocuind coordonatele primului punct (xt - , x " , x ® ), obținem Ui ~ Zi- (VI ) Respectiv Zz = Uz și Zz = Uz ■ (VI ) Când înlocuim coordonatele celui de-al patrulea punct în ecuație, obținem Vig \u d Zi Vz + Zg Vg + Zig / ZM = U" - tyi ~ Zyj- (VI ) Cu substituție succesivă într-un polinom de ordinul trei pentru un amestec cu trei componente Y \u d * + Zz * + Zz * + * * + * * + Zr * * + * * (* - * z) + + * * (* ~ *b) + * * (* - *b) + * * * coordonatele punctelor matricei de planificare (vezi Tabelul ) avem: Zi - Ui > Zg - Ug • Zz - Uz "(VI ■ ) Ung - Vі / + Ug vz + g/g + */ ? (VI ) U = Ui vz + Vg / + / - -, • (VI ■ ) Însumând ecuațiile (VI ) și (V ), obținem Zi \u d (Yng + Ui - Ui - Ug) - (VI ) Scăzând ecuația (VL ) din (VI ), obținem pentru : \u d \u d '/" (Suts - Zuі - yi -f- y,) (VI ]) Înlocuind polinomul redus pentru coordonatele punctelor și , obținem: ? \u d * / "(Y + Yrzz - Yr - Uz) • (VI ) za - (Zuggs ^ Vrzz Uz + Uz) • Și în mod similar, după înlocuirea coordonatele punctelor și ?із - alt (Uiz + Uіzz - Uі - Uz) • (vi ) іz \u d (Zutsz - Zu - Ui +uz) Înlocuind coordonatele ultimului punct, ținând cont de relațiile (VI ) - (VI ), obținem Uіrz ~ Vі vz + Yr '/s + Uz vz + (Uig + Uіrg - Ui - Uz) + ~ - / ?іrz > P = y - r'/ (Uz + yigr + Ujz + Ui + Ugz + Ugzz) + + "/ (Ui + Ug + Uz) - (VI- ) Coeficienții polinomului de ordinul trei pentru amestecul q-component Y \u d y Mi + (* "-x /) + l ?i=Ui etc , (VI ) ? \u d y - yi - yi etc , (VI ) ? \u d y - (y Yi + Uz) + (yi + Ug + Uz) • (VI - ) Pentru amestecul de componente g Y= Р (VI ) ?i = UIT d (VI ) = ^ - /, - / , (VI ) і \u d / s (-Ui + yi - / + Ug), Este \u d '/s (- Ui + usz - ț / i + y ), (VI • ) І z \u d / s (- Ug + uggz - / + Uz) • іg \u d ' / z (- Ui + Vipg - buig + J / - Ug) • (VI • ) іz \u d / z (- Ui + ushz - buіz + uіzzz - Uz), (VI ) ®gz \u d / s (- Ug + ygz - ygz + ygzzz - Uz) ”(VI ) ?urz \u d (Zi/i s - Uіzzz - Uіgzz) + /s (Оуі - Uz - Uz) I® (Uіr + Uіz) - - / s (boupig + tufiș - Z / - ZUіzzz - Uggz - Uzzz)\u e (VI • ) ? \u d ( x / - Yngz - Vygzz) + / s (bYg - Yi -Uz) - b (Yig + Uіz) - - / s ( yi + i / - t / J - Z / - Y - Ui s) • (VI • ) Pirsz \u d ( / - Ure - Viggs) +v / Ui -Ug) - I® (Uіz + Uzz) - - / (Sy + ^ - Zi/u - Zi/ gm - Ush Vizz ) • (VI • ) Pentru amestec //-component Y= ^,Xi+ S ?MX|X/+ іux/x/(хі - xy) + l хі*Дх( — Xy) + + l ) (VI ) hj \u d / s (- Ui + tyiiij - §Uіi + UiSh - Ui), (VI ) Ріі/Л \u d (Zushchk - uschk - yijkk) + » / (bui - yj - yk) - (уц + уік) - - v / z (tyiiij + - Zuzz - Zuikkk - yjjjk - yjkkk), (VI ) ?ijjk = ( УіЯй - Уііік - Уііккі + /з (tyj - Уі - Ук) - uts + tjjk) - - v / z fiyijjj + tyjjjk - tyiuj - tyjkkk - ushk - ytkkk) '(VI ) Pijkk = (Зуі]кк - уіі]к - уц}к) + / ( yk - Уі - Уі) - (yik + yjk) - /z (tyikkk T~ ^Uikkk ^Tsіііk ^UL)k - Y iii] - Y ijjj) > (VI ) ?ijkl - tyijkl - + Uііkk + Uіkki + Uіkki + Uіііі + Uіііі + Uіkіі) + For/z (Uіііі + Uіііk + Uіііі + + Uіііі + Uіііk + Uііііі + Uіkkk + Vjkkk + Ukkі + Uіііі + UNI + Ukііі)■ (VI ) După determinarea coeficienților ecuației de regresie, este necesar să se efectueze o analiză statistică a rezultatelor obținute: se verifică caracterul adecvat al ecuației și se construiesc intervale de încredere pentru valorile răspunsului prezise de ecuația de regresie Atunci când se realizează un experiment pe modele simplex-latice, nu există grade de libertate pentru a verifica caracterul adecvat al ecuației, deoarece aceste modele sunt saturate Pentru a verifica adecvarea, experimentele sunt efectuate la așa-numitele puncte de control suplimentare Numărul punctelor de control și coordonatele acestora sunt legate de formularea problemei și de caracteristicile experimentului În același timp, încearcă să prevadă posibilitatea utilizării punctelor de control pentru a îmbunătăți modelul în caz de inadecvare Precizia predicției răspunsului nu este aceeași în diferite puncte ale simplexului Varianta valorii răspunsului prezis poate fi determinată din legea acumulării erorilor Să arătăm acest lucru prin exemplul unui polinom de ordinul doi pentru un amestec cu trei componente În acest caz, presupunem că x, - sunt determinate fără erori, varianța de reproductibilitate Sy este aceeași în toate punctele planului, iar valorile răspunsului sunt rezultatul medierii l și pi a experimentelor paralele la nivelul corespunzător punctele simplexului Atunci varianțele lui y și y sunt egale; = (VI ) ȘI (VI ) cu Să înlocuim în polinomul redus Y \u d Mі + Ma + i'z-^z + + ?іzХі*z + ?az* *z coeficienți prin expresia lor în termeni de răspunsuri: ₽і \u d Ui , \u d tyij - ^Ui - ty} - Apoi primim Y \u d Y X + Yr x + Uz*e + ( Tu - Ti - a) * іXa + ( ^іz - Ui - Ue) * і * z + + ( f/ - Ya - ys) xaxs = yt (Xj - xjXa - rx^) + y (xa - xjXa - xax ) + + z (*z - *i*z - * * ) + n*i*a + із*і*е + аз*а*з• (VI ) Folosind condiția хі+хг+хз = , transformăm coeficienții la Xj - X]Xa - XJX = xj - xj (xa + Xg) = Xj - xj (I - xj = = x ( x - ), etc (VI ) obține Y = *i ( xj - ) y~i + *a ( xj - ) + *z ( xz - ) + xjxa i + XJX + xx z (VI ) Prin introducerea notaţiei la = xt ( xr - ), ar; = xrxy, (VI ) ținând cont de relațiile (VI ) și (VI ), obținem o expresie pentru dispersie : În mod similar, se obțin expresii pentru polinoame de ordinul trei, al treilea și al patrulea incomplet Unde eu bі=~~хі ( xf- Xj+l)- ^]x}, /= bi} = xtxj ( xj + xy - ), bijk = XiXjXh (VI ) (VI ) (vi ) (VI ) Pentru un polinom de ordinul al treilea NE PIL ' cijk nijk (VI ) Unde q \u d Va * i (Zxi - ) (Zxi - ), Cii; = '/a*i*;( xr - ), ssh = '/ *i*/( * - ctjk = xtxjxk (vi ) (VI ) (VI ) (VI ) Pentru un polinom de ordinul al patrulea Unde l (VI ) pentru un polinom incomplet de ordinul trei Кі , , , , , , , , , > , , , , , , , , , j/ , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , V , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , > , , , , , , , , , Continuare Numărul experimentului Denumirea răspunsului ■^ex Du₽ VP -Iex >p d>p t₽ tn U , , , , , , , , , ^ , , , , , , , , , și , , , , , , , , , ^ , , , , , , , , , U , , , , , , , , , U , , , , , , , , , Coordonatele punctelor de control sunt alese astfel încât să se poată construi un polinom de ordinul al patrulea dacă ecuațiile de regresie (VI ) și (VI ) sunt inadecvate (vezi Fig , d) Pentru fiecare punct de control, a fost compilat un raport t: t Ay Kn sy vm Eroare de reproductibilitate la determinarea reactivității cocsului sȘ - , , la determinarea porozității - sJJ - , Numărul de grade de libertate D- Numărul de experimente paralele în fiecare punct este n - Pentru fiecare punct de testare, a cantitatea „= l , se pot folosi Lambrakis palans, rețelele Scheffe simplex obișnuite, dar nu includ componente pure în aceste rețele, ci în schimb se pot realiza experimente la q puncte cu coordonate X =x = Orez Planul Lambrakis De exemplu, la construirea unui polinom de gradul doi într-un sistem cu patru componente, patru puncte cu coordonatele хі = х? \u d \u d x \u d X \u d (a se vedea Fig , a) înlocuiți cu patru puncte cu coordonatele xi \u d xi \u d xs \u d X \u d '/s (Fig ) Astfel planul Lambrakis în loc de patru experimente la vârfurile tetraedrului, include patru experimente la centrele triunghiurilor care formează acest tetraedru (хігз, Хі , хіз și Х ) și șase experimente la centrele fețelor tetraedrului (хір, хіз, Х , Х , Х și Х ) Planificarea centroidului simplex Planurile de centroid simplex ale lui Scheffe conțin ’- puncte, dintre care q sunt componente pure, Q sunt amestecuri cu două componente, Q sunt amestecuri cu trei componente etc și o observație este pentru fiecare amestec de componente ^ Coordonatele punctelor din planurile centroide simplex ( , , , ), ( / , / , , , ), ,( /q, / > T 'A 'A J Uz I 'A 'A 'A > R 'A 'A 'A > '/ 'A yi 'A 'A 'A U 'А > 'А 'А 'А >гз 'A 'A > '/ '/ / '/ > 'A 'A >^ Folosind proprietatea de saturație a planului, substituind succesiv coordonatele punctelor experimentale - - în polinom (VI ), determinăm coeficienții polinomului: = h = yt> ₽ = g/s; fa=y* ?m ~ ^Ui - yr ?іz \u d "/іz - Ui - uz ?I \u d J / - /! — і/ ₽ \u d / gz - u - i / ₽ = і/ — і/ — yt ₽ \u d "/z - i / s - і / Різз = і/ з - (і/ + j/ S + j/гз) + (yi + yt -} - y ), ₽і = і/ - (і/ + y + і/ ) + ( /i + Ug + Uz) > P = {/! - ( /! + yXi + Um) + ( / Ch " / + Uz), ?gz = / - ( /gz + Yu, + I/sz) + ( /" + Uz + Uz) (VI ) (VI ) (VI ) ₽ М - у - (і/іае + i/і + / + / ) + (і/ + yl + s/ - j/ a - + / + I/ ) - (I/! + y + Uz + Ui) (VI ) În mod similar, pentru polinomul (VI ) al amestecului //-component, avem: $i = Ui, (VI ) = *Уіі - уі - y}= [ ut - (y{ + y))], (VI ) \u d j / i> ft - (yi + uіk + yjk) + (yi - y} + uk) - = [ uik - (y{) + Uіk + Uіk) + (Ui + Ui + Uk)] (VI ) РіЛаі - uі)kt - (uі)k + ui)t + wikt + y)kt) - (uts + l/ife + + Uіt + UѢ + U)t + Ukt) - (Ui + U] + Ub + Ut) - [& uі]kt - (ui]k - + Uі)t + Uikt) + (Uі) + Uіk + Uіt + Uіk + Uіt + Ukt) ~ — (Ui + Ui + Uk + Ut)]- (VI ) În cazul general, formula pentru coeficienții ecuației de regresie obținute din designul centroidului simplex este ₽" ="• (VI ) unde r este numărul de indici ai coeficientului B*, ; SS, - suma rezultatelor experimentelor tuturor amestecurilor de componente r, luate în proporții egale ( /i) De exemplu, pentru coeficientul yuk avem r- (/, j, k) și trei sume: Yi + yj + yk-SS! — pentru Uіі + Uіk + Uіk = SSj — pentru Uііk \u d SS - Pentru / \u d i / /*= (VI ) (VI ) (VI ) Prin urmare, PiA = ((- l)®" l"" SSt + (- )®"* "-i SS + (-i)*~® "-" SS ] = = [( /i + Ui + Uk) - (yi + Uik + Uik) + yiih] (VI ) Verificarea adecvării ecuației de regresie obținută prin planul centroid simplex și construirea intervalelor de încredere pentru valorile proprietăților prezise de ecuație se realizează în aceleași moduri ca și în metoda rețelei simplex Exemplul S-a studiat influenţa compoziţiei asupra activităţii (yi) şi rezistenţei (yi) a unui catalizator de platină pe un suport metalic neporos la °C Masa totală a componentelor de la experiență la experiență a fost menținută constantă După ce a acceptat pe unitate, putem scrie că L x, - , unde x, - - componenta Pt / AIROz - zdrobit- catalizator de reformare uzat; xg și xe - componente - oxizi anorganici ai metalelor din grupele II și III ale sistemului periodic de elemente ale lui D I Mendeleev Soluţie A fost aplicat un design centroid simplex pentru - Matricea de planificare și rezultatele experimentelor sunt prezentate în tabel Numărul experimentului h • Xz y, % y , % , , , , , , , , , , , , , , , , Conform formulelor (VI ) și (VI ), coeficienții ecuațiilor de regresie pentru activitate L Ui - Xi + * + xs + xjX - Ex^ + * * + + , хіХ Х (VI ) și rezistența catalizatorului y \u d xі x xd - xjX + xiXg x Xg -f- XiX Xg (VI ) Eroare de reproductibilitate la măsurarea activității catalizatorului jy, - , , la măsurarea rezistenței sy - , Adecvarea ecuațiilor de regresie (VI ) și (V ) a fost verificată prin testul t Student la punctele de control și Yu (tabel) Număr de experiență * * xs U, U U , , , , , , , Pentru toate punctele de control, valorile testului r pentru un nivel de semnificație de p - , sunt mai mici decât valoarea tabelară Pe fig prezintă linii de valoare egală a activității catalizatorului yi și a rezistenței uz, construite conform ecuațiilor (VI ) și (VI ) Cea mai mare activitate a catalizatorului corespunde regiunii în care valorile componentei xt > , O putere de % este destul de satisfăcătoare De cel mai mare interes sunt punctele situate la intersecția liniei de producție egală yi - % cu linia de producție egală yi - % Experimentul (vezi tabelul și Fig ), efectuat în zona specificată, a dat o concordanță bună (în cadrul erorii experimentale) între rezultatele calculate și cele experimentale Planificarea unui experiment în studiul secțiunilor locale de diagrame Când se studiază diagramele compoziției-proprietăți ale amestecurilor de componente ^, devine adesea necesar să se investigheze dependența proprietății de compoziție nu în întregul interval de modificări ale concentrației componentelor în secțiunea locală a diagramei: O p) + • • • + Xg ( p - r ) , r* = ' + b ( ) - °) + * (r ) - °) + •••+*"( " ~ ")- (VI ) = xP + x gr + • • ■ + xP r'” = xp p + x rp + x r + xp r + • • • + x = xrr + xp p + xp ( > + + x( ) r(,)= xP gr+xp gr+xprz+ + x^ gr = o (VI ) x * + x r = O x(x) ( ) > x( = , , gP - , g - z} ) = , O bgrP' - z( ) - zj ) = , (VI ) lz^> -O?' * - z* ) = , , g = , z =l , (VI ) z = - , r® = Folosind soluțiile obținute în sistemul de ecuații (VI ), obținem formulele pentru legătura dintre coordonatele naturale xy și sistemul de coordonate Zy: gh = - , x - x r \u d , x , (VI ) Zg - Z | Z - Xg Înlocuind (VI ) în (VI ), obținem ecuația de regresie în coordonatele inițiale: L y \u d , + , x - , x + , x xs - , / + , x' - - , fd - , x x | + , x| x* + , l^ - , x^ + + , x^ - , x x| - , x^ + , xs (VI ) Pentru ușurință în utilizare, au fost trasate linii de temperaturi egale conform ecuației de regresie (Fig ) K CO Orez , linii de temperatură egală Zona studiată este un poliedru În prezența restricțiilor privind modificarea concentrațiilor componentelor, regiunea studiată formează în general un anumit poliedru La construirea unui plan de experimente, este necesar să se distribuie cumva punctele experimentale în funcție de rezultatul obținut din condiție și b x xs — — , x x -]- , x xs — , x x — — xsx (VI ) Deoarece dependența proprietății de compoziție este descrisă în mod adecvat de ecuația de regresie de ordinul doi, s-a dovedit a fi posibilă determinarea condițiilor optime prin aplicarea metodei programării neliniare Condițional Wies care oferă luminozitatea maximă a strălucirii au fost determinate în conformitate cu restricțiile (VL ): Vmax = , cu * = , , x = , , x = , , x = , Pe măsură ce numărul componentelor amestecului crește, numărul de puncte experimentale din planul McLean și Anderson crește rapid Pentru a reduce numărul de experimente, se pot exclude unele dintre centrele fețelor sau astfel de puncte, după excluderea cărora cele rămase se dovedesc a fi distribuite mai mult sau mai puțin uniform pe regiunea studiată D-Planuri optime Dintre diferitele criterii binecunoscute pentru optimitatea planurilor, cele mai importante sunt cerințele optimității D și G D-optimal este un proiect care minimizează volumul elipsoidului de împrăștiere al estimărilor coeficienților ecuației de regresie Proprietatea G-optimalitate oferă cea mai mică variație maximă a valorilor de răspuns prezise în zona de studiu Modelele simplex-redare au proprietățile optimității D și G numai atunci când construiesc polinoame de ordinul doi și al treilea incomplet Modelele Scheffe de ordin superior nu sunt D-optimale Rețeaua simplex D-optimă pentru un polinom de ordinul trei a fost mai târziu construită de Kiefer Dacă luăm în considerare un set de planuri cu coordonatele punctelor *i= , *> = xfc = , Xf = - Xj = b, xx = , b xіx> + W ("i - *>) + Kі" Kі l, pe această hartă de contur (pentru comparație, vezi Fig , b pentru o rețea simplex simplex) Cu toate acestea, este mai dificil să se pregătească amestecuri cu un conținut de componente de , și , decât amestecuri cu un raport de componente % și 'A Un dezavantaj semnificativ al planurilor considerate D-optimale de ordinul al treilea și al patrulea este lipsa compoziționalității Pe fig arată locația punctelor experimentale în planurile D-optime pentru sisteme cu trei componente Dacă este necesar, puteți trece de la un plan de ordinul doi sau un al treilea incomplet (Fig , a, b) la un plan de ordinul al patrulea (Fig , d), păstrând în același timp proprietatea D-optimalitate În acest caz, punctul din centru poate fi folosit ca control Când treceți de la un plan de ordinul al treilea (Fig , c) la un plan de ordinul al patrulea, trebuie fie să refaceți experimentele, fie să lucrați conform unui plan care nu are proprietatea D-optimalitate Xj(H O) % x [(NH ) HPO J a Orez Zona de studiu a vâscozității în sistem (NH p$zo> - ' > (VI ) ( °) q( ) = -■ (VI ) unde este factorul corespunzător din ecuația (V ); - corespunzatoare coeficientul se obţine în ecuaţia (VI ) Înlocuind (VI ), (VI ) și (VI ) în (VI ), p( ) p( ) + ■ ■ '~ T-(VI• ) Ecuația finală pentru regresia vâscozității soluțiilor pe temperatură și compoziție are forma L y \u d ( , - , T) zx - ( , - , T) z - ( , - , T) z - - (- , - , T) zxz - (- , - , T) zxz - ( , - , ') z za - - ( , - , T) zxz (zx - z ) - ( , - , T) zxz (zx - z ) - - ( , - , T) z z (z - z ) - ( , - , T) zxz z (VI ) Să obținem ecuațiile de legătură dintre pseudocomponentele zj și variabila naturală xj Sistemele de ecuații (VI ) pentru problema luată în considerare au forma , g - g * > - , g (p \u d Оz} » - , z' ) - , zf > = Ogr -Ogr - g = Soluțiile sistemelor (VI ) au forma gr = , , grP=O, r} > = , r = , , (VI ) gr = , gr = Înlocuind (V ) în sistemul de ecuații (VI ), obținem rx = , ( - x - x ), r = , x , (VI ) z \u d - zx - z \u d , xg - , * - , Exemplul S-a determinat compoziţia optimă a unui solvent multicomponent utilizat în procesul de curăţare a drojdiei de hidrocarburi Principalul indicator al epurării drojdiei este conținutul de hidrocarburi din biomasă după extracție (y) Pe baza considerentelor tehnologice și tehnice și economice, planificarea experimentului a fost realizată pe o secțiune locală a triunghiului de concentrare (Fig ) În zona de studiu, conținutul în amestec (%): acetonă x, + , zj ) + lzj ) = , z{'> - zJ > - , z} ) = O zj ) + z + zj > = (IV ) , ^ + , r = , g + z^>+l, z + z^ > + z = Ca urmare a rezolvării sistemului (VI ) avem z} '= , ; zj )= , ; r = , ; r* = , ; z = - , Folosind soluțiile găsite, obținem formulele pentru legătura dintre coordonatele naturale x și sistemul de coordonate z: zt \u d - , + , x + , x , z \u d , - , x - , x , (VI ) r \u d - yy - r \u d - , + , x + , x Verificarea adecvării ecuaţiei de regresie obţinută a fost efectuată după criteriul Studentului la patru puncte de control Rezultatele testului de adecvare sunt prezentate în tabel X, X -*z Y Y / * tab , , , , , , , , , , , , , , , , , , , , , , , , ' , , , , , , , Ecuația (VI ) este adecvată pentru a experimenta la un nivel de semnificație p = , Calitatea produsului rezultat este considerată satisfăcătoare dacă conținutul de hidrocarburi reziduale din biomasă este mai mic de , % Pentru a determina compozițiile solventului care asigură această cerință, conform ecuației (VI ), se construiesc linii de valori egale (Fig ) Compoziții de solvenți care îndeplinesc cerința y n ■ Au fost propuse principiile utilizate în alegerea proiectelor adecvate mai devreme de Box și Draper Draper și Lawrence și-au făcut planuri sisteme componente pentru grade de polinoame u, = , n = , u = și u = Pentru comoditatea construirii planurilor, Draper și Lawrence introduc un nou sistem de coordonate În sistemele cu trei componente, în planul triunghiului de concentrare (xp x , Xg), noul sistem de coordonate dinat a coincis cu triunghiul, dar această axă (Fig (xi, x , xz) și dreptunghiulară (z,, Zj) ) are următoarea relație Orez Sistem de coordonare nat pentru planurile Draper-Lawrence pentru trei și patru este ales în așa fel încât originea coorului să fie centrul de greutate al triunghiului, unul dintre vârfuri să fie pe axa z , iar celelalte două să fie simetrice față de - ) Între coordonatele triunghiulare r, = - Kz - I - * i - x + x (VI ) ȘI * - Zz ( j -Z Pz"-|-/p), = V ( z^ - z m), (VI ) xs - ( za V tn), unde /u~ este lungimea laturii triunghiului de concentrare Când se studiază întreaga diagramă, r = , când se studiază secțiunile locale ale diagramei, m ); ) vârfurile dreptunghiului (c, d), (-c, -d),'(c, - d), (-c, d) După implementarea unuia sau altui plan Draper-Lawrence pentru sisteme cu trei componente, polinoamele sunt construite pentru două variabile independente z și z de ordinul întâi („i = cu u = ) A Y \u d bd * * b z (VI ) sau de ordinul doi (n] = cu u = ) Y = ba + ftjZj + b z + b z z -zf -^ • (VI • ) Tabelul Set de puncte Număr de experimente în centru „ Număr total de experimente N Parametri ( , ) p = , , g = , ( , ) p- , , g - , ( , ) р- , , g = , ( , ) р = , , g „ , ( , ) p \u d , , a - , ( , ) p = , , b- , ( , ) p \u d , , c \u d V , Im - t / * ( , , ) p-\ \u d , , g \u d , p - , ( , , ) p \u d , , gt - , g, - , * Valoarea lui d este aleasă în mod arbitrar Pentru a construi un polinom de gradul I în raport cu sistemele cu trei componente ( (VI ) x \u d g / i (- Z + Z + Z + m)\u e X = /t (- Zl - Z - Z + m), Punctele planului pentru sisteme cu patru componente sunt selectate (în coordonatele zb z , z ) din următoarele seturi: ) vârfuri ale unui tetraedru asemănător unuia de concentrație cu coordonate de vârf (a, a, - a), (a, - a, a), (- a, a, a), (- a, - a, - a); ) vârfurile tetraedrului ( , b, b), (b, - b, - b), (b, b, - b), (-b, - b, b); ) puncte de pe axe (±/i, , ), ( , ztft, ), ( , , ztft); ) vârfuri de tetraedre cu coordonate (-r, -s, -/), (-r, s, i), (r, -s, /) (r, s, -t); (-i, r, s), (t, -r, s), (t, r, -s), (-t, -r, -s); (—s, t, r) (s, -t, r), (s, t, -r), (-s, -t, -r) Tabelul Set de puncte Număr de experimente în centru n Număr total de experimente N Parametri ( , ) a = , , b = , ( , ) a = , , b = , ( , ) a = , , b = , ( , ) a = , , b = , (h ) a \u d , , b \u d , ( , ) a = , , h = , ( , ) a = , , h = , ( , ) a = , , h = , ( ) g = , , s = , t = , , ( ) /•= , , t = , t = , , După implementarea unuia sau altuia plan Draper-Lawrence pentru sisteme cu patru componente, se construiesc polinoame pentru variabilele independente sin zp z și z de ordinul întâi (u, = cu u = ) L Y - '^O + & Z + & Z + b Z (VI ) sau de ordinul doi (", = pentru n = ) Y = ba + &izi + ^ Z + b z b yy + b z z + b z z * + ^llzf +^ z| H”^ z| • (VI ) Parametrii (în fracțiuni de m) ai unor planuri Draper-Lawrence care conțin cel mult puncte (pentru q = , u = , q = ) sunt dați în Tabel Tabelul Set de puncte Număr de experimente în centru Po Număr total de experimente N Parametri ( , , ) a \ b az r s t ( , , ) , , , , , ( , , ) , , , , , ( , , ) , , , , , ( , , ) , , , , , ( , , ) a , b , /• , , t , ( , , ) , , , , , ( , , ) , , , , , ( , , ) ' , , , , , ( , , ) a , h , r , s , t , ( , , ) , , , , , ( , , ) , , , , , (G, , , ) a\ , az , b , h , ( , , , ) , , , , ( , , , ) , , , , ( , , , ) , , , , ( , , , ) , , , , ( , , , ) , , , , ( , , , ) , , , , ( , , , ) ay , az , , valorile parametrilor r, s, t și N sunt determinate din sistemul de ecuații r + s + / \u d V "i / , (VI ) primul = , adică pentru sisteme cu trei componente, puncte cu coordonate ( z „ z ) ar trebui luate, iar pentru sistemele cu patru componente - cu coordonate ( Z,, z > Oz^ Parametrul Ѳ depinde de eroarea aleatorie și de coeficienții polinomului și este aproape de unitate dacă eroarea aleatorie nu domină pentru sistemele cu trei componente și pentru sistemele cu patru componente sisteme Să transformăm, de exemplu, pentru a minimiza eroarea totală, planul ( , , ) pentru # = , u, = , u = , dat în Tabelul Coordonatele primului punct, ținând cont parametrul Ѳ= , , se obțin astfel: x \u d \u d i / (- ! - Ѳr Kz + "i) \u d / s (- , - , UT-G ) \u d , , x * \u d / s (+ - vgvT ) / n) \u d і / z (+ - , - , - , - , / z - | - ) \u d \u d , , x \u d і / (gbg / + t) \u d Vz ( - , - , / z + ) \u d , Planificați complet ( , , ) pentru # = , u = , u = , minimizând eroarea totală, este dat în tabel Tabelul Matricea de planificare ( , , ) pentru W , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Tabelul Cuantile ale distribuției normale r - r r , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Tabelul Cuantile de distribuție t a lui Student Numărul de grade de libertate f Niveluri de semnificație p , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Continuare Numărul de grade de libertate f Niveluri de semnificație p , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , OO , , , , , , , Tabelul Cuantile de distribuție Pearson X j p Numărul de grade de libertate f Niveluri de semnificație p , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , p s , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , - , , , , , , , , , , , , , , , , , , , , , , , , , , , I , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Continuare Număr de grade de libertate / niveluri de semnificație p , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Tabelul Cuantile ale distribuției Fisher Fy-p pentru p - , fi OO , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , Continuare L ' OO , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , s, s , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , OO , , , , , , , , , Tabelul Cuantile ale distribuției Cochran G*\p pentru p == , n f OO I * Toate quantilele G]~p sunt mai mici decât unu, prin urmare, în tabel arată numai săptămâna zecimală după punctul zecimal ^ înainte de care, atunci când utilizați tabelul, trebuie să puneți zero „-’ De exemplu, pentru n - , /" avem , - , Tabelul Ranguri semnificative ale testului de ranguri multiple al lui Duncan la p , „Oh r , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , ] , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , oo , , , , , , , , , , , , , , , , Tabelul Matricea echilibrului aleatoriu Numărul experimentului * l * l L * x L Liu LI Li - + - + + - - - - - - - - - + - - - - - - - - - + + - + - - - - - - - + - - - + - - - - - - + - - + - - - - - - - - - + - + - - - - - - + + - - - - - - - - - - + + - - - - - - - - - - + + - - - - - - - - - + + - - - - - - - - + + - + - - - - - - - - - - - - - - - - - - + + + + - - - - - - - - + - - - - - - - - - + - - - - - - - - + - - - - - - - - - Continuare Numărul experimentului *, A -Xz * * * Xj *ѳ * -Hyu L- * + + + + - + - + + + - + + - + - + - + + - - - + - + - - + + -I -I - - + - + - + + + + + - - + + - + + + - - - + + + ' + + + - - - + - + - + + + + - + + + - - + + + - - + - - - -I - - -I Tabelul Matricea echilibrului aleatoriu Numărul experimentului * * L- * L x *a * Liu L * * * * L, * * * ) - + + + -J + + + - + - - - - - I - - + - - + + - + + - - + + - + - - - - - + + + + - + + + + - + - - - - - - + + + + + - + - + - + + + - - - - - - - - - + + - + + + - - - + - - - - - + - + - + + - - - + - - + - - - - - + + - + + + - - + + - - - - - + + - - + + - - - - + - + - - - - - - - - + - - - - + -t + - + - - - - - - + - - - - - - - - - - - - - - - - - - + - - - + + + + + + + - + - - - - - + + + - - - + - + + - + + - - - - - - - + - + + - + - + + + - - - - - + + + - - - + - + - - - - - - - - - + + + + + + - + + - - - - - - - - + + - + - + + + + + + - - - - - - - - + + - + - + - - + - + - - - - - - + - - + - + - - - + + + + + - - - - - - + - + - + - - - - + + - - - - - - - + + + + - - + - +r - - + - - - - - - - - - - - - - + - - - - - + - - - - - - + - + - - + - + - - - - - + - - - - + + - - + - - - + - - - - - - - - - + + - - - + + + - + - + + + - - - - - - - + - + - + + - + + - - + - - - - - + - - - + - + + + - + + + - - - - - - + - - + + - - + + + - - - - - - - + + + - + - + - - - - - + - - - - - - - - - + + + - - + + + + + - - - - - - - + + + + - - + + + - - + - - - - - - - + + + + - - - + - - + - - - - - - - + + - + + + + - + + - + + - - - - - - CUPRINS dezirabilitate ( ) Planuri complexe Factorizați experimentul \ combinat cu pătratul latin ( ) Metoda de planificare secvenţială simplex ( ) Planuri Plackett-Berman saturate ortogonale ( ) Experimente de screening Metoda echilibrului aleatoriu ( ) Planificarea unui experiment în determinarea constantelor ecuațiilor cineticii formale ( ) Planificarea unui experiment într-un mediu de producție ( ) Exerciții Capitolul VI Planificarea unui experiment în studierea diagramelor compoziție-proprietate Metoda rețelelor simplex ( ) Design Scheffe Simplex Lattice Planificarea centroidului simplex Planificarea unui experiment în studiul secțiunilor locale de diagrame ( ) D-Planuri optime ( ) Planuri cu minimizarea părtinirii sistematice ( ) Planificarea unui experiment atunci când se studiază dependența unei proprietăți de raporturile componentelor ( ) Exerciții Cuvânt înainte Introducere Literatură Aplicație Prima parte METODE DE ANALIZĂ STATISTICĂ A EXPERIMENTULUI Capitolul I Principalele caracteristici ale variabilelor aleatoare Variabile aleatorii Axiomele teoriei probabilităților Legile de distribuție ( ) Caracteristici numerice ( ) Proprietăți ale așteptării și dispersiei matematice ( ) Distribuție uniformă ( ), Distribuție normală ( ), Sisteme de variabile aleatoare ( ) Stochastic conexiune ( ) Exerciții Capitolul II Definirea parametrilor funcției de distribuție Populația generală și eșantionarea aleatorie ( ) Metoda maximei probabilități ( ) Estimarea așteptării și varianței matematice ( ) , Clasificarea erorilor de măsurare ( ) Legea adunării erorilor ( ) Erori de măsurători indirecte ( ) Determinarea varianței din măsurătorile curente ( ) Intervalele de încredere și probabilitatea de încredere ( ) Testarea ipotezelor statistice ( ) Yu Estimarea așteptării matematice a unei variabile aleatoare distribuite normal ( ) Estimarea dispersiei unei variabile aleatoare distribuite normal ( ) Comparația a două dispersii ( ) Comparația mai multor varianțe ( ) Comparația a două medii ( ) Comparația mai multor medii ( ) Verificarea omogenității rezultatelor măsurătorilor ( ) Comparație între distribuția eșantionului și distribuția populației ( ) Criteriul de bunăstare a potrivirii < ( ) Testul Wilcoxon ( ), Testarea ipotezei de normalitate pe un set de eșantioane mici ( ) Exerciții Capitolul II Analiza varianței Problema analizei varianței ( ) Analiza unidirecțională a varianței ( ) Analiza bidirecțională a varianței ( ) Planificarea unui experiment în analiza varianței latină și hiper-greacă-latină quadra- tu ( ) Cuburi latine (І ) Exerciții H Capitolul IV Metode de analiză de corelare și regresie Coeficient de corelație eșantion ( ), Coeficienți de corelație parțială ( ) Regresia aproximativă Metoda celor mai mici pătrate ( ) Regresia liniară de la un parametru ( )' Regresia parabolică ( ) Polinoame Cebyshev ( ) Regresia transcendentală ( ) Estimarea etanșeității unei conexiuni neliniare ( ) Metoda corelării multiple ( ) Yu Analiza regresiei sub formă de matrice ( ) Metoda grupului noua prezentare a argumentelor (MGUA) ( ) Metoda componentelor principale ( ) Exerciții Partea a doua METODE DE PLANIFICARE EXPERIMENTALĂ Capitolul V Experiment factorial complet ( ) Replici fracționate ( ) Optimizare prin metoda ascensiunii abrupte de-a lungul suprafeței de răspuns ( ) Descrierea zonei apropiate de extremum Planuri compoziționale ale lui Box - Wilson ( ) Planuri ortogonale de ordinul doi ( ) Planuri rotative de ordinul doi Box-Hunter ( ) - Criterii pentru optimitatea planurilor ( ) Investigarea suprafeţei de răspuns Rezolvarea problemei de optimizare ( ) Funcția Svetlana Lazarevna Akhnazarova, Viktor Viaceslavovici Kafarov Metode optimizarea experimentului în tehnologia chimică Cap editat de S F Kondrashkov Editor T S Kostyan ml redactor S M Erokhin Artistul V M Borovkov Editor de artă T M Skvortsova Redactor tehnic A K Nesterova Corector S K Zavyalova IB nr Ed N? Khim - Predat setului / / Semnat pentru a apăsa / / Format X '/iv Bum birou carte-jurnal Font literar Imprimare offset Volumul , arb cuptor l , arb kr -ott , ed l Tiraj de exemplare Zach Pret rub k Editura „Școala Superioară”, , Moscova, GSP- , str Neglinnaya, / Fabrica de poligraf Yaroslavl a Soyuzpolygraphprom sub Comitetul de Stat URSS pentru editare, tipar și comerț cu cărți , Iaroslavl, st Libertate, 